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Abstract 


A  critical  review  of  the  breakup  of  liquids  is  given,  with  emphasis  oh 
liquid  drops.  Evidence  is  presented  that  the  breakup  is  a  rate  process 
arid  hence  must  be  considered  from  a  time-dependent  approach  rather 
than  from  purely  dimensionless  parameters,  such  as  Weber  or  Bond 
numbers. 

A  theoretical  model  of  the  aerodynamic  breakup  of  liquid  drops  is  de¬ 
veloped,  postulating  that  the  dual  breakup  mechanisms  (bag  and  stripping! 
of  liquid  drops  result  from  pressure  and  frictional  drag  on  the  drop 
The  derivation  of  a  quantitative  expression  is  presented,  predicting 
aerodynamic  breakup  times  of  liquid  drops  as  a  function  of  the  gas -Mow 
properties  arid  the  physical  properties  of  the  drop.  A  quantitative  ex¬ 
pression  is  also  developed  that  predicts  the  mass-mean  diameter  of 
droplets  produced  by  aerodynamic  breakup  of  the  drop  as  a  function  of 
the  gas -flow  properties  arid  physical  properties  of  the  drop 

Using  shock-tube  techniques,  a  comprehensive,  experimental  parametric 
study  was  conducted  oh  the  breakup  times  of  drops  of  water,  mercury, 
bis,  arid  three  silicone  fluids  with  different  viscosities  The  results 
ire  compared  to  the  available  theories,  and  are  shown  to  be  in  very  good 
agreement  with  the  preceding  theory 

Ah  experimental  study  of  the  size  distribution  of  droplets  produced  bv 
the  breakup  of  bis  droplets  was  conducted  for  various  gas-fiow  velocities 
and  initial  drop  sizes.  The  results  are  shown  to  be  in  good  agreement 
with  the  developed  theory. 
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PRODUCTION 


Knowledge  of  the  dispersion  (atomization,  disseminatm1..  or  breakup) 
c  naracteristics  of  liquids  is  of  importance  in  a  wide  \.iriot'-  of  applica¬ 
tions,  including 

•  Liquid-fuel  atomization  m  rockets,  internal-combustion  engines, 
and  gas  turbines 

•  Insecticidal  and  agricultural  spraying 

e  Spray  drying 

•  Meteorological  studies 

•  Production  of  aerosols  for  therapeutic  work  and  chemical/biological 

warfare 

The  mechanism  of  dissemination  of  liquids  is  net  well  understood  There 
is  no  general  theory  available  with  which  to  estimate  the  particle  -  size 
distribution  of  a  liquid  as  a  function  of  the  physical  properties  of  the  liquid 
01  to  evaluate  the  atomizing  technique  employed  and  its  operating  conditions 
Nevertheless,,  many  features  of  the  dissemination  process  have  been  eluci¬ 
dated  from  considerable  experimental  work  and  some  limited  theoretical 
work  performed  over  the  years  (References  1  through  5 >. 

The  factors  that  influence  the  breakup  of  liquids  are,  to  some  extent,  a 
function  of  the  techniques  used  for  performing  the  breakup,  e  g.  ,.  explosive 
dissemination,  nozzle  dissemination,  and  electrical  dispersion  Generally, 
factors  that  have  been  found  to  influence  the  breakup  of  liquids  include  the 
surface  tension,  viscosity  and  density  of  the  liquid,  the  density  and 
viscosity  of  the  gas  (or  other  fluid)  through  which  the  liquid  is  moving,  the 
relative  velocity  between  the  liquid  and  gas,  the  wavelength  of  any  surface 
disturbances  resulting  in  instability,  and  the  hydrodynamic  properties  of 
the  gas  flow  (laminar,  turbulent,  sub-  or  supersonic,  etc.  ) 

The  physical  properties  of  the  liquid  and  of  its  environment  define  part 
of  the  p-eceding  factors  in  a  given  case,  whereas  the  techniques  employed 
for  dissemination  define  the  remainder  of  the  factors  For  example,  it  is 
well  known  that  the  efficiency  of  atomization  and  dispersion  of  a  liquid 
through  an  injection  valve  depends  both  01  the  properti'-s  of  the  liquid 
(for  a  given  shape  of  the  orifice)  and  on  tne  shape  of  the  orifice  (for  a 
given  iiquidl 
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Th.'e  is  evidence  that  atomization  of  a  liquid  occurs  in  several  stages, 
and  that  there  is  little  difference  whether  the  breakup  of  the  liquid  is 
brought  about  by  the  emergence  of  liquid  into  still  air  at  h;oh  speed  from 
a  nozzle,  or  by  the  interaction  of  a  stream  of  liquid  and  .  fa st -flowing 
gas  stream.  The  various  stag  is  occur  in  the  following  order- 

a.  Stretching  of  the  liquid  into  films,  sheets,  streams,  or  jets  as  a 
result  of  accelerating  the  liquid  by  some  prescribed  means. 


b.  Initiation  of  small  disturbances  at  the  surface  of  the  liquid,  in 
the  form  of  local  ripples,  protuberances,  or  waves. 


c.  Formation  of  ligaments  as  a  result  of  air  (fluid)  pressure  and 
shearing  forces. 


d.  Collapse  of  ligaments  into  drops  as  a  result  of  surface  tension 


e.  Further  breakup  of  these  drops  in  movement  through  the  air  or 
other  fluid. 


This  report  is  concerned  with  a  theoretical  and  experimental  investigation 
of  the  last  of  these  stages,  i.e.  ,  the  aerodynamic  breakup  ,f  liquid 
drops.  The  studies  have  emphasized  the  mechanism  and  rate  of  the  breakup 
process,  with  some  attention  directed  to  the  resulting  particle  size  of 
the  dispersed  aerosol.  Also  considered  was  le  role  of  the  physical  prop¬ 
erties  of  a  liquid  drop  on  its  dispersion  characteristics. 


2.  SUMMARY 

A  critical  review  of  literature  on  the  aerodynamic  breakup  of  liquids  is 
given,  with  emphasis  on  liquid  drops.  It  is  shown  that  the  breakup  of 
any  geometry  liquid  may  be  considered  to  occur  in  two  main  phases: 

(1)  the  shearing  of  sheets,  films, or  jets  from  the  liquid  by  the  gas-flow 
forces,  and  (2)  the  rapid  breakup  of  these  liquid  streams  by  instabilities, 
which  grow  exponentially  with  time.  The  breakup  of  liquids  is  shown  to  be 
a  rate  process  and  hence  must  be  considered  f;oin  a  time -dependent 
approach  rather  than  from  purely  dimensionless  parameters,  such  as 
Weber  o’-  Bond  numbers 

A  theoretical  model  of  the  aerodynamic  breakup  of  liquid  chops  is  de¬ 
veloped,  postulating  that  the  dual  breakup  mechanisms  (bag  and  stripping) 
observed  for  liquid  drops  result  from  pressure  and  frictional  d’-ag  on  the 
drop.  The  frictional  drag  is  usually  greater  than  the  pre-  sure  drag; 
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however,  the  original  deformation  of  the  drop  produced  by  the  aerod\nair,;c 
flow  produces  a  sharper  curvature  at  the  outer  edges  of  the  drop  than 
in  the  middle  of  the  drop.  This  curvature  produces  a  higher  surface- 
tension  pressure  at  the  edges  than  in  the  middle  of  the  drop.  The  surface- 
tension  pressure  resists  liquid  deformation  by  the  gas  pi  essure,  and 
hence  drops  undergo  bag  breakup  for  velocities  near  the  critical  velocity 
required  for  breakup.  For  higher  velocities,  however,  the  friction  stress 
at  the  drop  edges  becomes  larger  than  the  surface-tension  pressure;  the 
drop  then  undergoes  stripping  breakup,  s;nce  the  inction  stress  is 
greater  than  the  pressure  stress.  For  certain  conditions,  the  drop  ex¬ 
hibits  the  characteristics  of  both  bag  and  stripping  breaknn. 


Quantitative  expressions  are  derived  to  predict  the  breakup  times  lor 
both  the  bag  and  stripping  breakup  of  liquid  drnp<?  with  given  physical 
properties  (size,  surface  tension,  and  viscosity)  subjected  to  aerodynamic 
flow  of  a  given  velocity  .  These  expressions  are  based  on  how  quickly  the 
middle  of  the  drop  can  be  sheared  away  from  the  edges  of  the  drop  (bag 
breakup),  and  how  quickly  the  edges  of  the  drop  can  be  sheared  away  from 
the  middle  of  the  drop  (stripping  breakup).  For  most  practical  situations, 
however,  one  expression  is  sufficient  to  describe  the  breakup  times  as 
accurately  as  they  can  be  measured.  The  breakup  time  t  is  shown  to  be 
given  by 

_  d _ 

|{A2  +  BP)1/2  -  aJ 


where 


P  = 


,2  _ 


2<r 

d 


d  =  the  original  drop  diameter 

q  -  viscosity  of  drop 

<r  -  surface  tension  of  drop 

Oj  =  density  of  drop 

=  density  of  gas  flow 
u  =  velocity  of  gas  flow 
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For  low-viscosity  liquids  and  negligible  surface-tension  forces  (either 
low  surface  tension  or  high  velocity), the  preceding  equation  simplifies 
to 


For  very -high-viscosity  liquids  and  negligible  surface-tension  forces, 
the  equation  becomes 


t 


32  q 


2 

u 


It  is  shown  that  the  proposed  breakup  model  can  probably  be  modified  to 
include  the  erosion  of  solid  particles  in  gas  streams. 


A  theoretical  model  is  developed  for  the  breakup  of  the  liquid  film  (or 
sheet)  sheared  from  the  drop  by  the  gas  flow;  this  breakup  produces  the 
droplets  that  are  formed  by  the  breakup  of  the  drop.  The  liquid  pieces 
produced  by  breakup  of  the  film  (and  which  subsequently  form  into  drops 
by  surface-tension  forces)  are  considered  to  have  a  thickness  that  is  de¬ 
termined  by  the  shearing  of  the  film  from  the  drop,  a  width  determined  by 
divergence  (expansion)  of  the  sheet  as  it  leaves  the  drop,  and  a  length 
determined  by  the  growth  of  instabilities  produced  in  the  lengthwise  direc¬ 
tion  (parallel  with  the  gas  stream)  of  the  sheet  by  the  gas  stream.  The 
mass -mean  diameter,  D,  of  droplets  produced  by  the  breakup  of  the 
original  drop  is  given  by 


D 


T  136  q  <r3/Z  d  1,2  1 

I  ,  —\ 

L  Pa  Pi  u  J 


1/3 


(The  velocity  and  initial  drop-diameter  dependence  predicted  in  this  equa¬ 
tion  has  previously  been  found  experimentally  by  Weiss  and  Worsham.) 

A  comprehensive,  experimental  parametric  study  has  been  conducted  on 
the  breakup  times  of  drops  of  water,  mercury,  bis.  and  three  silicone 
fluids  with  widely  different  viscosities,  employing  various  flow  velocities 
and  drop  sizes.  The  liquids  were  chosen  to  include  a  w.de  variation  oi 
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physical  properties,  which  were  investigated  using  shock-tube  and  high¬ 
speed  photographic  techniques  Descriptions  of  the  experimental  apparatus 
and  test  procedure  .'.re  presented  in  detail,  and  the  o:  .ginal  experimental 
data  are  appended.  Experimental  breakup  times  are  compared  to  avail¬ 
able  theories,  and  are  shown  to  be  in  very  good  agreement  with  the  de¬ 
veloped  theory.  The  results  are  also  in  agreement  with  the  Gordon 
theory,  if  the  surface-tension  pressure  used  by  Gordon  is  decreased  by 
a  factor  of  two  to  four.  Experimental  photographs  of  the  breakup  process 
demonstrate  that  the  bag  and  stripping  breakup  mechanisms  are  separated 
by  a  transition  region  that  includes  characteristics  of  both  modes  of 
breakup  behavior  as  predicted  Various  characteristics  of  the  breakup 
process  are  discussed, including  the  change  of  drag  coeinoent  with  time 


An  experimental  study  was  performed  to  sample  and  assess  droplet-size 
distributions  resulfng  from  the  breakup  of  single  drops  of  bis  subjected 
to  various  aerodynamic  flow  velocities  and  exhibiting  both  mechanisms  ol 
breakup.  The  mass-mean  diameter  of  the  resultant  droplets  was  found  to 
be  in  good  agreement  with  the  derived  equation. 


It  is  concluded  that  the  present  breakup  theory  is  probably  adequate  to 
describe  the  breakup  times  of  liquid  drops  of  any  liquid,  the  derived  mean- 
drop-size  expression  is  probably  adequate,  but  needs  further  checking 
using  other  liquids.  Further  studies  are  also  required  to  elucidate  theo¬ 
retically  the  entire  drop-size  distribution  function  Associated  conclusions 
and  recommendations  are  also  discussed 


3.  DISPERSION  BACKGROUND 


3.  1  DISPERSION  OF  BULK  LIQUID 

Breakup  of  a  moving  jet  of  liquid  was  first  treated  theoretically  by  Rayleigh 
(Reference  6)  and  later  extended  by  Weber  (Reference  7,.  in  these 
treatments  small  perturbations  or  disturbances  acquired  in  the  slip 
streams  at  the  nozzle  ;or  base  of  the  liquid,  are  shown  to  be  carried  for¬ 
ward  with  ine  jet  at  the  same  time  growing  exponentially  until  violent 
instability  occurs  followed  by  breakdown  of  the  liquid  coiumn.  At  low 
jet  velocities,  the  air  does  not  appreciably  affect  the  shape  of  the  jet,  anu 
rotationaily  symmetric  ^varicose;  disturbances  are  formed  on  the  cylindri¬ 
cal  jet  The  disturbances  are  assumed  to  grow  in  conformance  with  the 
equation 


S 


Sj  exp  < q* ,  cos  <2  it  x / \ • 


(1) 


whe  r  •» 

S  =  amplitude  of  the  disturbance  at  time  t  and  distance  x 
from  the  nozzle 

Sj  =  average  initial  amplitude  of  the  disturbance  imparted  to 
the  jet  stream  at  its  base  (the  nozzle) 

q  =  rate  of  growth  of  the  disturbance  of  wavelength  X 

Breakup  of  the  jet  into  droplets  occurs  where  the  jef  is  nicked  off  by  the 
growth  of  the  varicosity.  The  wavelength  lh«.l  dominates  the  lit  tkup  is 
the  optimum  wavelength,  X  Qpt,  for  which  the  rate  of  growth  is  a  maxi¬ 
mum,  i.  e.  ,  qm-  The  breakup  criteria  is:  Sj  exp  =  dj?2,  where  d. 

is  the  jet  diameter.  Breakup  time,  t^,  is  then  given  by 


tb  =  d/qm>  In  (dj/2Sj) 


(2) 


The  optimum  wavelength  according  to  Weber  is 


where 


q  -  viscosity  of  liquid  jet 
cr  =  surface  tension  of  liqiid  jet 
p  2  =  density  of  the  liquid  jet 

Assorting  the  average  jet  diameter  is  equal  to  the  diameter  of  the  jet 
orifice,  the  average  diameter,  d,  of  the  spherical  drops  initially  produced 
by  breakup  of  the  cylindrical  jet  is 


!-5dj2  Xopt 


1/3 


The  growth  rate  of  the  optimum  disturbance  is  given  by 


H) 


'/qm  =  ( pj/tr)0'  ^  dj  1,5  +  3  ti«Mt 


(5) 


-6- 


The  breakup  distance  is  -  u  where  is  the  jet  velocity  Thus, 
increases  with  increase  in  jet  velocity  white  t^  is  independent  of 
velocity 

As  the  velocity  of  the  jet  is  further  increased,  the  breakup  i  to  drops 
becomes  partially  influenced  by  interactions  of  the  jet  with  the  air.  The 
air  velocity  increases  over  the  wave  crests  and  decreases  over  the 
troughs;,  simultaneously,  the  pressure  decreases  over  the  crests  and 
increases  over  the  troughs  The  net  result  is  that  the  optimum  wavelength 
for  breakup  decreases,  and  the  growth  rate  of  the  optimum  wavelength 
disturbance  increases,  leading *o  a  decrease  in  breakup  time  breakup 
distance,  and  resulting  drop  size  with  an  increase  m  jet  velocity.  As 
the  jet  velocity  is  still  further  increased,  the  jet  becomes  sinuous  iwavy* 
in  nature  as  a  resu't  of  increa  ;ed  surface  interactions  with  the  air.  the 
faster  the  jet  velocity  the  more  the  wave  motion  is  intensified.  The 
wave  motion  leads  to  rapid  breakup  of  the  jet  into  drops  However,  the 
treatment  of  this  case  by  Weber  is  semi  - empe ncal  in  nature,  and  does 
not  lead  to  simple  equations  by  which  the  breakup  characteristics  of  the 
jet  may  be  readily  predicted 

The  theoretical  results  of  Rayleigh  and  Weber  have  been  verified  by 
various  experimental  studies  (References  1  through  5i  As  shown  by 
the  study  of  Haenlein  (Refererce  8<  four  characteristic  breakup  forms 
were  evidenced  as  the  jet  velocity  was  increased  viz.  drop  formation 
without  air  influence,  drop  formation  with  air  influence  wave  formation 
leading  to  drop  formation,  and  immediate  complete  disintegration  of  the 
jet  into  drops.  The  latter  situation,  which  was  not  discussed  by  Rayleigh 
and  Weber  is  generally  the  situation  that  occurs  m  practical  applications 

Castleman  (Reference  9'  considered  the  breakup  of  a  high-velocity  jet  to 
occ^r  through  the  formation  of  ligaments  or  threads  which  are  drawn 
from  the  main  jet  mass  and  collapse  because  of  their  instability  into  a 
number  of  drops.  The  ligaments  arise  through  air  friction  on  the  jet 
an  effect  that  causes  surface  distu.bances  to  grow  according  to  the 
Rayleigh-Weber  theory  At  higher  airspeeds,  finer  ligaments  are  formed 
and  break  up  to  form  smaller  dropr,  higher  surface  tension  would  cause 
qu-cker  collapse  of  the  ligaments  before  ihey  are  drawn  too  finely  and 
would  result  in  larger  drops  These  trends  are  verified  by  experiment 

Fogler  and  Kleinschmidt  -Reference  10.  consider  the  liquid  to  be  drawn  out 
into  thin  flat  sheets  or  films  rather  than  ligaments  or  threads  The 
waves  m  the  film  builc  up  rap-dly  and  ciuse  a  whipping  of  the  surface  so 
that  it  curls  back  on  itself  to  form  a  hol.ow  tube  Such  a  tube  is  unstable 
and  breaks  off  immediately  into  a  series  of  hoMuw  spheres,  which  form 
droplets 
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Schweitzer  (Reference  11)  stressed  the  importance  of  turbulent  flow  on 
the  breakup  of  high-velocity  jets.  The  turbulence  gives  the  liquid  a  radial 
component  of  velocity,  helping  to  overcome  the  surfaco-tensic-.  forces, 
and  providing  surface  disturbances  for  interactions  with  the  air.  Schweitzer 
considered  turbulence  to  be  the  dominant  factor  in  the  breakup  of  very -high- 
velocity  jets,  with  air  fi  iction  enhancing  the  process,  he  concluded  that 
viscosity  is  the  most  important  liquid  property  that  influences  the  breakup 
of  liquid  jets.  The  turbulent  breakup  mechanism  would  help  explain  the 
fourth  characteristic  breakup  form  observed  by  Haenlein. 

The  views  of  Schweitzer  regarding  the  effects  of  turbulence  on  the  breakup 
of  high-velocity  jets  are  also  supported  oy  the  work  of  various  Oii.e. 
workers  (References  1  through  4),  including  Merrington  and  Richardson 
(Reference  12)  who  found  that  the  mass -mean  droplet  diameter,  d,  depended 
primarily  on  the  jet  velocity  and  the  liquid  viscosity;  no  influence  of  nozzle 
size  or  shape  was  detected.  Their  experimental  data  at  high  velocities 
satisfied  the  empencal  equation 


dUj  =  500  (q/pj)°-2 


(6) 


At  lower  jet  velocities,  the  drop  size  reached  a  limiting  value,  as  expected 
for  varicose  breakup  (refer  to  Equation  4). 

Ohnesorge  (References  1  and  3)  concluded  that  the  breakup  of  liquid  jets 
occurs  quite  generally  ir  three  different  stages,  with  the  particular  stage 
depending  on  the  value  of  the  Reynolds  number,  Re,  of  the  jet.  These  stages 
consist  of  varicose  breakup,  sinuous  breakup,  ar.d  direct  atomization  from 
the  nozzle;  the  stages  are  separated,  consecutively,  in  a  plot  of  the 
Ohnesorge  number,  Z,  versus  Re,  using  the  lines  given  by 


z  = 

i . 025  - 

5.  08  (in-4)  Re 

(?a) 

z  = 

1.430  - 

1.  42  (10-4)  Re 

(7b) 

z  =  r)/  {(T  pj  dj)*^2  (7c) 

Re  =  u,  u  d  /  n 
1  j  J 

These  relationships  show  that  the  jet  flow  chant---.,  'rom  a  laminar  to  a 
turbulent  behavior  as  the  velocity  is  increased. 
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Dityakin  and  Yagodkin  (Reference  13)  concluded  that  fluctuations  of  the 
jet-flow  velocity  and  the  density  of  the  ambient  medi"m  lead  to  a  decrease 
in  the  average  droplet  size  produced  by  jet  breakup,  a-:  wall  as  to  a  dis¬ 
tribution  of  droplet  sizes. 

Dunne  and  Cassen  (Reference  14)  investigated  the  instabi'ity  of  a  liquid 
jet  with  a  velocity  discontinuity  (shock)  superimposed  on  the  jet.  They 
found  that  the  perturbation  grows  in  time,  and  is  propagated  through  the 
jet  as  a  thin-disc,  rotationally  symmetric  wave  whose  velocity  is  the  mean 
of  the  instantaneous  particle  velocities  immediately  in  front  of  and  behind 
the  discontinuity.  Eisenklam  and  Hooper  (Reference  15)  suggested  that 
the  breakup  of  very -high-velocity  jets  could  be  partially  produced  by  the 
internal  pressure  of  the  liquid,  and  could  periodically  alternate  with  the 
turbulent  breakup  mechanism. 

Middleman  and  Gavis  (Reference  16)  found  that  capillary  jets  of  liquids 
expand  or  contract  upon  ejection,  depending  upon  conditions  of  ejection 
and  fluid  properties.  Straubel  (Reference  17)  showed  that  an  electrical 
potential  placed  across  an  injector  effected  fine  atomization  of  a  previously 
coarse  stream.  Miesse  (Reference  18)  found  that  a  transverse  ambient - 
pressure  oscillation  brought  about  rapid  mixing  of  a  pair  of  adjacent  streams. 


Weiss  and  Worsham  (Reference  19)  studied  the  drop  sizes  obtained  on 
injecting  a  cylindrical  liquid  jet  of  a  synthetic  wax  into  a  large,  hot  air- 
stream  sustained  at  high  velocity.  The  relative  velocity,  ur,  between  the 
liquid  jet  and  the  airstream  was  found  to  be  the  variable  having  the  largest 
effect  on  particle  size.  The  mass-median  diameter  of  the  drop  size 
was  found  to  be  approximated  by 


v.^  =  injection  velocity  of  the  liquid 

D.  ,  =  diameter  of  the  jet 

inj  J 

Pao  =  density  of  the  air  at  300°F  and  1 -atm  pressure 

n  -  air  density 

a 
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The  inv .  rse  velocity  exponent  1.  33  compares  to  a  value  ol  1.  0  found 
earner  by  Nukiyama  and  Tanasawa  (Reference  20)  for  small  pneumatic 
atomizers,  and  to  the  value  of  1.  11  for  a  metal  alloy  and  1.  68  for  a  wax 
found  by  Marshall  (Reference  21)  for  large  venturi  atomizers. 

Weiss  and  Worsham  concluded  that:  "The  atomization  of  liquids  by 
large,  high-velocity  airstreams  occurs  by  direct  action  of  the  airstreams 
on  the  exposed  liquid  surface.  Therefore,  the  relative  velocity  between 
the  liquid  and  the  airstream  is  of  primary  importance.  Physical  prop¬ 
erties  ol  the  fluids  do  affect  spray  fineness,  but  their  net  influence  is 
less  critical.  The  exact  way  in  which  the  liquid  is  introduced  to  the  air, 
i.  e.  ,  the  geometry  and  operation  of  the  injecto*-,  is  of  least  importance, 
particularly  at  very  high  air  velocities  and  for  customary  range  of 
variables.  " 

The  preceding  conclusions  appear  to  be  generally  applicable  under  all 
c  itions  in  which  high-velocity  gases  are  involved,  with  the  exception 
of  conditions  where  high-vapor-pressure  liquids  are  concerned  (in 
which  case  vaporization  may  also  be  an  important  parameter). 

The  breakup  of  a  tangentially  moving  liquid  sheet,  or  a  conical  film 
such  as  may  be  produced  by  a  swirl  type  of  injector,  has  been  considered 
by  various  investigators  (References  22,  23,  and  24)  to  occur  in  much  the 
same  manner  as  that  of  a  plain  jet.  .small  disturbances  in  the  sheet  develop 
tc  ‘orm  waves,  which  rapidly  grow  and  break  up  the  sheet;  in  addition, 
holes  often  appear  in  the  sheet.  Squire  (Reference  23)  deduced  the  optimum 
wavelength  to  be  given  by 


k  opt 


4  Tra/pa 


2 

u 


(9) 


where 

na  =  air  density 

u  =  relative  tangential  velocity  cf  the  air 
The  maximum  growth  rate  is 


" 

r 

r 

2  nu/\>pt 

1  2  o-/p,  u2  b 

L 

. 

(10) 
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where  b  is  the  film  thickness.  The  average  drop  size  produced  by  the 
breakup  is 


d 


1/3 


(11) 


Equation  9  may  also  be  used  for  estimating  drop-size  values  in  the  sinuous 
breakup  of  liquid  jets.  It  is  of  significance  that  Equation  1  1  predicts  the 
drop  size  to  vary  inversely  as  the  4/3  power  of  the  relative  velocity,  as 
found  by  Weiss  and  Worsham  (Equation  8). 

Borodin  and  Dityakin  (Reference  25)  showed  that  the  possibility  of  several 
predominant  wavelengths  leading  to  breakup  could  lead  to  a  distribution 
of  drop  sizes.  Ycik  et  al  (Reference  24)  assumed  that  each  wave-formed 
ripple  on  a  liquid  sheet  behaves  as  a  circular  jet  of  radius  (b\/u)l^  and 
obtained  for  the  drop  size 


d  =  2(bXopt)l/2  (12) 

Equations  11  and  12  have  both  been  shown  to  be  in  fair  agreement  with 
experimental  data. 

Taylor  (Reference  26)  considered  the  acceleration  of  two  fluids  in  contact, 
and  showed  that  the  boundary  between  high-  and  low -density  fluids  will  be 
unstable  when  the  acceleration  (or  force)  is  directed  from  the  low-density 
fluid  to  the  high-density  fluid.  This  destabilizing  influence  is  known  as 
Ta^lor_Instabilit^.  Acceleration  from  the  denser  towards  the  lighter 
fluid  exerts  a  stabilizing  tendency.  The  relative  tangential  velocity  on  a 
flu.d  is  always  a  destabilizing  influence  (known  as  Helmholt^Instabilit^). 
Richtmyer  (Reference  27)  has  considered  Taylor  instability  in  the  shock- 
acceleration  of  compressible  fluids. 


3.  2  DISPERSION  OF  LIQUID  DROPS 

Lane  (Reference  28)  conducted  one  of  the  first  comprehensive  studies  ol 
the  aerodynamic  breakup  of  liquid  drops  in  an  air  stream.  He  subjected 
a  water  drop  to  an  essentially  steady  (gradually  increasing)  air  stream, 
and  found  that  it  becomes  increasingly  flattened;  at  a  critical  velocity  of 
air,  the  drop  was  blown  out  'n  a  concave  manner  into  the  form  of  a  hollow 
bag  attached  to  a  roughly  circular  rim.  Bursting  of  this  bag  produced  a 
shower  of  very  fine  droplets,  and  the  rim,  which  contained  at  least  70%  of 
the  mass  of  the  original  spherical  drop  Voke  up  later  into  larger  drops. 
This  type  of  breakup  is  known  as  bag  breakup. 
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Lar.e  pointed  out  that  the  phenomenon  has  special  interest  in  connection 
with  the  atomization  of  sprays  in  that  in  both  cases  the  mass  fraction  of 
very  small  particles  produced  is  small,  and  that  these  small  pa.  tides 
result  from  the  breakup  of  a  stretched  film.  He  found  that  the  cm. cal 
velocity,  uc,  required  to  break  the  drops  was  given  by:  (uc  -  v)^  d  =  612, 
where  v  is  the  velocity  of  the  drop  at  breakup  and  d  is  the  diameter  of 
the  orig.nal  drop.  A  relationship  of  the  form 


Cd  pa  (u-v)2/2  =  4<r/d 

was  expected  on  theoretical  grounds,  where  Cd  is  the  drag  coefficient  of 
a  sphere.  This  work  led  to  a  value  of  the  constant  of  about  1200  rather 
than  612.  However,  since  the  drop  does  not  remain  spherical  but  approxi¬ 
mates  a  lens  shape  before  it  bursts,  the  drag  coefficient  should  be  closer 
to  that  given  for  a  circular  disc  {about  twice  the  value  for  a  sphere), 
yielding  a  value  of  about  600  for  the  constant. 

Equations  of  the  form  of  Equation  13  have  also  been  used  by  other  early 
investigators  of  droplet  breakup  (References  1  through  4).  The  interpre¬ 
tation  of  this  equation  is  that  the  relative  flow  velocity  is  brought  to  rest 
in  front  of  the  drop  at  the  stagnation  point  and  is  converted  to  a  pressure 
This  pressure  pushes  the  interior  of  the  drop  more  than  the  sides  of  the 
dro-  (as  a  result  of  pressure  distribution),  so  as  to  form  a  bag,  which 
then  breaks.  Viscosity  was  found  to  influence  the  breakup  only  when  it 
was  very  great,  and  then  only  tended  to  retard  the  breakup. 

The  breakup  of  drops  subjected  to  abrupt,  fast  (transient)  air  blasts  was 
also  studied  by  Lane.  In  this  case,  the  drop  deformed  m  the  opposite 
direction  to  that  of  bag  breakup,  and  formed  a  convex  surface  to  the  flow 
of  air,  the  diameter  being  about  twice  that  of  the  original  spherical  drop. 
The  edges  of  the  saucer  shape  were  first  drawn  out  into  a  thin  sheet,  then 
into  thin  filaments  that  broke  to  form  drops.  Tins  type  of  breakup  is 
known  as  shear  breakup  or  stripping  breakup. 

it  was  found  that  the  velocities  required  for  breakup  in  the  fast  blasts  were 
lower  than  in  the  steady  stream,  the  uivci6..,ice  .ncreasmg  for  small  drops 
This  effect  was  explained  by  some  work  of  Taylor,  who  deduced  that 
(uc  -v)  steady  r  (tic  '  v)  transient-  A  measurement  of  the  drop  sizes 

produced  by  the  breakup  of  a  drop  showed  the  resultant  mean-mass  drop 
size  decreased  with  an  increase  in  velocity. 


<H) 
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Calculations  by  Taylor  (Reference  28)  of  the  boundary-layer  thickness 
(and  hence  resultant  drop  size)  of  the  drop  undergoing  stripping  were 
roughly  in  agreement  with  experimental  data  obtained  rt  low  air  velocities, 
but  were  too  large  at  the  higher  velocities. 

Hinze  (Reference  29)  theoretically  investigated  the  breakup  of  liquid  drops 
in  a  gas  stream  by  an  examination  of  the  slight  deformation  of  the  globules, 
using  linearized  hydrodynamic  equations.  He  found  the  critical  condition 
for  breakup  of  a  nonviscous  liquid  suddently  exposed  to  a  constant-velocity 
gas  stream  to  be 


6/R  =  0.  1?  We 


41 


where 

6  =  deformation  of  the  drop  in  the  radial  direction 
R  =  initial  drop  radius 

In  Equation  14,  the  Weber  number,  We,  is  defined  by 

We  =  pa  u2  R  U  (15) 

where 

p  =  gas  (air)  density 

u  =  relative  velocity  between  the  drop  and  gas  stream 
a  -  surface  tension  of  the  liquid 


For  liquids  with  large  viscosities,  Hinze  found 


6/R  -  0.095  We  (16) 

Hinze  thus  associated  the  breakup  of  liquids  with  a  critical  value  of  the 
Weber  number.  For  the  cases  of  the  breakup  of  a  liquid  in  a  continuously 
increasingly  gas  flow  (e.g.  ,  a  falling  drop),  the  constants  in  Equations  14 
and  16  should  both  be  0.  095.  The  correlation  of  these  equations  with 
limited  experimental  data  of  Merringtcr.  and  Richardson  (Reference  12) 
suggested  that  the  critical  value  of  6/R,  i.e.  ,  (8/R)c  ..  is  generally  of  the  order 
of  unity  for  breakup;  however,  highly  viscous  liquids  required  larger  v-tlues 
(e.  g.  ,  about  2  in  one  case). 
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The  breakup  time,,  tb,  of  a  nonviscous  liquid  was  computed  to  be 


tb  =  (1.  16  R/u)  ^(p,/pa)  (6/R)cr 
and  that  of  a  highly  viscous  liquid  was  found  to  be 

tfa  =  (10qi/pau2)  (6/R)Cr 


(17) 


(18) 


It  should  be  noted  that  Equations  17  and  18  do  not  include  the  case  of  medium 
viscosity  effects;  in  addition,  the  values  of  (6/R)cr  required  vary  with  the 
viscosity. 

An  experimental  study  of  the  aerodynamic  breakup  of  liquid  drops  was 
conducted  by  Hanson  et  al  (References  30  and  31).  These  investigators 
found,  as  opposed  to  the  findings  of  Lane,  that  it  is  possible  to  produce 
bag  breakup  under  suddenly  applied  (transient)  flow  conditions  if  the  air 
velocity  is  only  slightly  greater  than  the  critical  velocity.  When  the  air 
velocity  is  considerably  above  the  critical  value,  shear  breakup  always 
occurs.  They  also  found  that  with  some  drops  undergoing  bag  breakup, 
the  bag  develops  a  reentrant  portion  near  its  middle;  this  so-called 
"stamen"  increases  in  length  with  time,  and  in  some  cases  completely 
in’  erts  the  bag  before  breakup  occurs. 

Also  observed  was  that  small  dark  spots  surrounded  by  concentric  rings 
appeared  in  the  bags  (it  was  not  known  whether  these  were  the  beginnings 
of  rupture).  A  general  conclusion  was  made  that  the  passage  of  a  normal 
air  shock  over  a  liquid  drop  does  not  of  itself  cause  breakup,  but  rather 
that  it  is  the  relative  velocity  between  the  drop  and  the  shocked  air  which, 
if  sustained  for  a  sufficient  time,  will  cause  the  drop  to  break  up.  The 
breakup  of  distilled  water,  methyl  alcohol,  and  three  silicone  oils  with 
viscosities  of  10,  50,  and  100  centistokes,  over  the  100-  to  700-micron 
range  of  drop  diameters  was  found  to  ob&y  the  following  emperical  equations:- 


2  , 

uc  d  - 

6.  21  ( 1 06> 

(19a) 

uc2d  = 

2.  71  (106) 

(19b) 

9.  69  (105) 

(29c) 

162d  = 

6.84  (i05) 

(lVd) 

l-36d  = 

2.  32  (105) 

(19c) 
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where  d  is  the  drop  diameter  in  microns,  u  is  the  gas  velocity  in  ft/sec, 
and  the  subscript  c  denotes  critical  quantity. 

The  dependence  of  critical  velocity  on  surface  tension  for  constant  drop 
size  was  found  *i  go  approximately  as  the  surface  tension  to  the  ore-third 
power.  This  effect  contrasted  with  the  one-half  power  variation  found 
by  Lane.  The  studies  showed  a  negligible  effect  of  viscosity  in  the  cri¬ 
tical  breakup  velocity  for  viscosities  of  about  10  centistokes  or  less, 
above  this  value,  however,  the  effects  became  significant.  Increasing 
viscosity  increased  the  critical  velocity  required  to  break  up  a  drop  ol 
given  diameter.  The  trend  became  more  pronounced  as  the  drop  diameto- 
was  decreased. 


Priem  (Reference  32)  studied  the  breakup  of  liquid  drops  by  shock  waves 
and  concluded  that  the  drops  were  broken  up  by  the  high  gas  velocity  behind 
the  shock  front.  A  similar  conclusion  was  reached  by  Sato  (Reference  33) 
on  the  basis  of  photographic  evidence. 


Magarvey  and  Taylor  (Reference  34)  studied  the  free-fall  breakup  of  large 
liquid  drops  and  concluded  that  the  regular  breakup  is  not  triggered  by 
internal  vibrations.  As  shown  in  Reference  35,  the  time  tn  (period) 
required  for  a  complete  oscillation  of  a  liquid  drop  fer  the  nth  spherical 
function  (harmonic)  is  given  by 


*n  = 


3  it  m  /  a  n  (n-1)  (n+2) 


1/2 


where  m  is  the  mass  of  the  drop. 


(20a) 


The  case  of  n  =  1  results  in  no  motion,  i.  e.  ,  the  shape  of  the  drop  is 
spherical;  the  case  n=2  is  that  of  the  slowest  vibration  and  represents 
deformation  into  what  is  close  to  an  ellipsoidal  shape.  In  this  case 


t 


n 


3  n  m  /  8it 


1/2 


(20b) 


The  cases  of  n  =  3,  4  .  .  .  ,  correspond  to  more  rapid  partial  vibrations  to¬ 
ward  the  ellipsoidal  shape;  these  vibrations  are  anharrnonic  because  their 
periods  do  not  have  a  direct  relation  to  the  period  of  the  ellipsoidal  de¬ 
formation. 


Engel  (Reference  3o)  photographically  studied  the  breakup  of  w:.ter  drops 
behind  an  air  shock.  She  found  the  drops,  which  unde -went  shear  breakup, 
to  first  flatten  on  the  leeward  side  subsequent  to  the  passage  of  the  shock 
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front.  This  affect  was  followed  by  a  radial  flow  of  the  drop  to  give  a 
ring  shape,  which  was  then  bent  in  the  wind  direction  by  the  airflow.  The 
deformation  was  accompanied  by  streamers  of  mist  torn  off  by  'ha  moving 
boundary  layer.  Flattening  of  the  drop  perpendicular  to  the  dnecls*,-j  of 
the  airflow  is  the  result  of  the  pressure  distribution  around  the  drop. 

According  to  Burgers  (Reference  36), the  preceding  effect  is  initially 
described  by 


<W  =  Pt2  /  3  R  P 


(21) 


where  douj  is  the  outward  displacement  at  the  equator  of  the  drop 
(perpendicular  to  the  direction  of  the  airflow  at  time  t),  p  is  the  density 
of  the  liqu  d  drop  of  radius  R,  and  p  is  the  pressure  difference  between 
either  the  windward  or  the  leeward-face  stagnation  point  and  points  around 
the  equator  of  the  drop. 

After  the  initial  rapid  flattening,  the  deformation  remained  approximately 
invariant  with  time,  and  finally  increased  again  (but  at  a  slower  rate  than 
the  initial  increase).  In  the  final  stages  in  the  breakup  of  the  drop,  a 
distinct  corrugation  appeared  in  the  windward  face,  an  effect  whicn  maybe 
caused  by  surface  waves  produced  by  the  wind  that  blows  out  radially 
around  the  stagnation  point.  The  crests  of  these  waves  may  break  or  be 
blown  off,  and  breakup  of  the  remaining  portion  of  the  drop  into  separate 
sections  can  then  occur. 

It  is  possible  that  a  hole  forms  in  the  center  of  the  liquid,  forming  a  ring 
that  breaks  into  droplets  It  was  found  that  large  drops  lag  behind  small 
drops  in  developing  various  stages  of  fragmentation  (this  behavior  decreased 
as  the  flow  velocity  was  increased).  A  change  in  velocity  was  found  to  be 
more  effective  on  the  rate  of  breakup  than  a  change  in  droj  diameter. 

Rabin,  Schallenmuller,  and  Lawhead  (References  37  and  38)  investigated 
the  aerodynamic  breakup  of  burning  and  non-burning  liquid  drops  They 
found,  as  did  previous  investigators,  that  the  flow  field  following  the  rela¬ 
tively  weak  shock  causes  the  droplet  breakup,  rather  than  any  impulsive 
action  of  the  front  itself  The  critical  velocity  required  to  break  up  a 
burning  drop  was  found  to  be  slightly  lower  than  for  a  non-burning  drop  (an 
effect  attributed  to  the  lower  surface  tension  of  a  burning  drep).  It  may 
also  be  noted,  however,  that  the  boundary  layer  around  the  horning  drop 
will  also  differ  frcrri  that  of  a  non-burning  drop  Both  bag  and  shear 
breakup  behaviors  were  observed. 
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An  analysis  of  the  data  revealed  that  the  breakup  mechanism  (bag  vs  shear) 
could  not  be  explained  on  the  basis  of  the  flow  duration  as  related  to  the 
drop  vibrational  period  (Equation  20b).  Data  of  Lane  <*.id  of  Hanson  et  at. 
(Reference  30  and  31)  also  were  not  explainable  on  this  basis.  At  velo¬ 
cities  greatly  in  excess  of  the  critical  value,  only  shear  breakup  was  ob¬ 
served.  at  elevated  gas  pressures,  only  shear  breakup  was  observed,  and 
the  critical  velocity  decreased.  The  critical  velocity  that  auld  just  cause 
breakup  for  short -duration  flow  was  found  to  agree  with  the  equation 


We  /  Re  1/2  "S'  0.  45 


(22) 


where  the  Reynolds  number  is  based  on  the  sphere  diameter,  and  the  Weber 
number  is  based  on  the  sphere  radius.  The  drag  coefficient  for  burning 
and  non-burning  drops  was  found  to  be  approximately  unity. 

Gordon  (Reference  39)  has  derived  an  expression  from  Newton's  second 
law  for  the  aerodynamic  breakup  of  liquid  drops  undergoing  bag  breakup 
(it  appears,  however,  that  the  results  are  also  applicable  to  shear  breakup, 
as  shown  later  in  this  report).  An  assumption  was  made  that  the  stagnation 
pressure  at  the  front  of  the  drop  extrudes  a  cylindrical  plug  from  the  drop, 
the  extrusion  is  retarded  by  surface  tension,  viscosity, and  inertial 
forces.  The  acceleration  of  the  plug  is  given  by 


dv  P  1 

dt  pjd  pjd 


l-,2  pau2  -  (8a/d)  -  (16r)  v/d) 


(23a) 


where  d  is  the  drop  diameter,  P  is  the  resultant  pressure  on  the  face  of 
the  drop  (composed  of  the  stagnation  pressure  minus  the  surface  tension 
pressure  minus  the  viscous  pressure),  u  is  the  relative  velocity  between 
the  drop  with  surface  tension,  a.  a-id  viscosity,  q  ,  and  the  air  stream  w-lh 
density, pa,  and  pj.  is  the  liquid  density.  Integrating  Equation  23  gives  the 
velocity  of  the  plug  as  a  function  of  time 


where 


v  =  (A/B) 


1  -  exp  (-Bt) 


> 


nau“  «<r  16q 


(c3b) 


(2L) 
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Integral. ig  a  second  time  gives  the  plug  displacement  as  a  function  of  time, 
and  considering  breakup  to  occur  when  the  displacement  is  equal  to  the 
drop  diameter  gives 

2  <16t,)2 

Pjd2  (pau2  -  16cr/d) 

where  t  is  the  breakup  time.  This 
Thus 

t  =  (2d/u)(Pl/pa)1/2  ;  for  n<  <du(paPl)1/2  (25a) 

<r<  <du2pa 

t  =  32q  /pau2  ;  for  q> >  du  (papj)  ^2 

<r<  <  du2pa  ( 25b) 


16q  t 
Ptd2 


1  +  exp  I 


-16rj  t  * 


(24) 


equation  simplifies  for  special  cases. 


Equations  25a  and  b  may  be  compared  with  Equations  17  and  18  as  derived 
by  Hinze.  Except  for  a  constant  factor,  the  equations  are  identical.  A 
useful  approximation  to  Equation  24  is 


t 


2dPll/2  32 n 

. . . .  . . .  i-  .  ..  .  . . 

(pauZ  -  16<r/d)1^2  pau2  -  16<x7d 


(26) 


Gordon  states  that  this  approximation  is  never  too  small,  and  is  at  most 
37%  too  large.  (Experimental  data  were  not  available  for  a  rigorous  test 
of  the  Gordon  theory  at  the  time  of  its  publication. )  His  theory  is  com¬ 
pared  to  results  of  the  present  experiments  in  a  later  section  of  this  report. 

Dodd  (Reference  40)  considered  the  bag  breakup  of  liquid  drops  and  argued 
that  the  critical  velocity  required  for  breakup  is  the  velocity  that  makes 
the  radius  of  curvature  of  the  deformed  drop  a  minimum.-  He  concluded 
that  the  form  of  the  condition  for  breakup  is  the  same  as  that  found  by 
Lane,  i.e. ,  Equation  13. 

Morrell  considered  the  breakup  of  liquid  drops  and  of  jets  by  transverse 
shocks  in  a  series  of  papers  (Reference  41,  42,  and  43).  His  results  with 
jets  are  directly  related  to  drops,  since  breakup  was  consider  :d  under 
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conditions  where  the  breakup  was  due  entirely  to  the  transverse  shock  and 
not  to  the  jet  instability  per  se  Originally  (Reference  41;  he  utilized 
the  method  of  Hinze  (Reference  29;  for  drops  and  derived  an  analogous 
breakup-time  expression  for  a  cylindrical  jet.  He  observed  experimen¬ 
tally  that  the  breakup  occurred  by  a  shear  type  of  mechanism  (the  breakup 
time  wa3  found  to  depend  on  d,/u).  However,  large  values  of  6/R  were 
required  to  make  the  data  agree  with  the  theory  Morrell  pointed  out  that 
his  data,  as  well  as  those  of  Rabin  and  Lawhead  tReferenee  37)  and  of 
Lane  (Reference  28),  showed  that  the  minimum  or  critical  Weber  number 
for  breakup  should  increase  as  the  initial  drop  tor  jet}  diameter  increases, 
whereas  the  data  reported  by  Hanson  et  al. (References  30  and  31)  and 
by  Volynskii  (Reference  44)  show  the  opposite  The  reasons  for  this 
discrepancy  are  not  apparent 

In  a  later  paper  (Reference  42),  Morrell  derived  an  expression  for  the 
breakup  of  drops  as  a  result  of  stripping  It  was  assumed  that  the  liquid 
issuing  from  the  periphery  of  the  drop  forms  a  sheet,  which  breaks  when 
the  frictional  force  on  the  sheet  is  equal  to  the  tensile  force  of  the  liquid. 
The  time-dependence  of  the  breakup  was  not  considered,  however.  The 
case  in  which  the  dynamic  pressure  causing  breakup  of  the  drop  decays 
exponentially,  l  e  , 


iU2  =  faiu,2 


exp  i 


-t/L 


<2 


was  also  considered  for  both  bag  and  shear  breakup  In  this  equation, 
the  subscript  i  indicates  initial  values,  and  ta  is  the  action  time  (time 
constant)  of  the  decay.  Morreii  assumed  that  when  the  action  time  is 
greater  than  the  natural  period  tn  of  the  drop  vEquation  20b)  or  jet.  the 
liquid  mass  will  oscillate  and  break  by  a  deformation  (bag,  mechanism. 
For  ta<tn,  the  mass  will  not  oscillate  and  will  break  by  a  stripping 
mechanism.  A  reasonable  correlation  of  the  data  was  obtained  using  this 
assumption.  However,  Rabin  and  Lawhead  tReferences  37  and  38/  had 
earlier  been  unable  to  corre!a»e  either  their  data  or  those  of  Hanson  et  al. 
by  means  of  a  similar  assumption. 

It  may  also  be  argued  that  the  drop  can  not  initially  know  either  the  pres¬ 
sure  duration  cr  its  decay  characteristics,  and  hence  cannot  know  what 
frontal  shape  to  assume  tconvex  or  concave)  to  help  establish  the  breakup 
mechanism 
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In  his  later  paper  (Reference  43),  Morrell  states  that  the  flow  duration 
(action  time)  merely  affects  the  extent  of  breakup.  Deformation  breakup 
appeared  to  become  increasingly  important  as  the  velocity  w-s  increased. 
The  deformation  6  was  determined  (photographically)  to  be  given  by 


6/R  =  2+4  We/Re1/2 


(28) 


where 

We  =  pau2R/<r 
Re  =  PauR/n 

R  =  the  jet  or  drop  radius 

Substitution  of  this  expression  in  the  equation  originally  derived  for  break¬ 
up  (Reference  41)  gave  a  complicated  expression  for  the  breakup  time  as 
a  function  of  the  Reynolds  and  Weber  number.  This  expression  always 
predicted  values  of  the  breakup  time  that  were  low  when  compared  with 
the  available  experimental  data.  The  similarities  between  Equations  28 
and  22  are  noteworthy,  and  the  possible  use  of  Equation  28  in  Equations  17 
and  18  derived  by  Hin2e  is  also  of  interest. 

.  lorrell  also  derived  a  rate  expression  for  the  stripping  breakup  of  jets 
in  his  later  paper,  assuming  the  liquid  to  be  stripped  off  in  sheets  per¬ 
pendicular  to  the  jet.  The  volumetric  removal  rate  of  liquid  from  the  drop 
was  taken  as 


dV 

(29) 


where  V  is  the  volume  of  liquid  per  unit  length  of  jet,  dbj  is  the  thickness 
of  the  liquid  sheet  (boundary  layer),  ujis  the  average  velocity  in  the 
boundary  layer,  and  the  right-hand  side  of  Equation  29  is  to  be  evaluated 
for  the  length  of  the  sheet  at  breakup..  This  length,  6,  was  assumed  to 
be  given  by  Equation  28. 

Integrating  Equation  29  gives  for  the  breakup  time 


t  =  ttR2  /  (2dbluj) 


(30) 
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A  modified  form  of  the  general  velocity  profile  suggested  by  Sandborn  and 
Kline  (Reference  45)  was  used  for  the  boundary  layer,  rather  than  the 
exponential  profile  used  earlier  by  Taylor  (Reference  26).  The  boundary 
layer  thickness  was  found  to  be  given  by 


dbl  =  (vj  x/ub)1,/2 
and  the  average  boundary  layer  velocity  by 


(31a) 


Uj  =  10  bu 


(31b) 


where 

b  =  0.04,  lpa/p!)2/3  (va/v!)1/3 


(31c) 


and  v  is  the  kinematic  viscosity,  q  /p. 

Substitution  of  Equation  31  in  Equation  30  gives  for  the  breakup  time 


(32) 


A  comparison  of  the  experimental  data  with  Equation  32  showed  satisfactory 
agreement,  except  near  the  threshold  (large  values  of  t).  It  was  in  this 
region  however,  that  the  breakup  was  found  to  occur  by  the  deformation 
mechanism.  Equation  32  may  be  expected  to  also  hold  for  liquid  drops, 
except  for  a  small  numerical  con-tant  (which  accounts  for  the  change  in 
geometry).:  The  equation  is  compared  with  experimental  data  in  a  later 
section  of  this  report.- 

A  theory  of  the  breakup  of  liquids  in  high-velocity  gas  streams  has  re¬ 
cently  been  presented  by  Mayer  (Reference  46),  who  considered  the  wind 
(air3tream)  to  induce  disturbances  (waves  or  ripples)  in  the  fluid.  Waves 
of  very  small  wavelengths  decay  because  of  viscous  dissipation,  and  very 
long  wavelengths  develop  slowly  because  of  inertial  :-i fects.  The  other 
waves  grow  at  an  exponential  rate,  which  is  determined  by  the  fluid 
properties. 
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It  was  pc  '.tulated  that  when  a  wave  of  a  certain  wavelength  has  grown  to 
an  amplitude  comparable  to  its  wavelength,  the  crest  of  the  wave  is  shed 
as  a  ligament  from  which  droplets  (whose  diameters  are  proportional  to 
the  wavelength)  are  also  formed. 

The  derived  mean  droplet  diameter  for  the  primary  breakup  of  the  liquid 
is 

d  =  71,2  B  |h  jta/pj)172  /  paua2|  2/3  {3 

where  B  =  F/B^^.  The  parameter  F  i3  associated  with  the  crest  configu¬ 
ration  at  the  instant  of  erosion,  ar.d  the  sheltering  parameter,  B,  is 
associated  with  the  portion  of  the  wave  crest  exposed  to  the  driving  effect 
of  the  wind.  The  value  of  B  is  conceptually  of  the  order  of  unity,  and  for 
one  specific  case  it  was  calculated  to  be  about  0.  3.,  Equation  33  predicts 
that  the  average  drop  size  should  vary  inversely  as  the  4/3  power  of  the 
relative  velocity  between  the  liquid  and  the  gas  stream  (precisely  the 
velocity  effect  on  drop  size  found  experimentally  by  Weiss  and  Worsham 
in  Reference  19). 

When  a  liquid  drop  mo/es  through  a  fluid  (or  vice  versa),  a  shearing  stress 
(skin  friction)  is  exerted  on  the  surface  of  the  drop,  an  effect  that  may  set 
up  an  internal  circulation  within  the  drop.  In  addition,  the  drop  will  be 
deformed  as  discussed  previously.  These  effects  significantly  influence  both 
the  drag  and  the  stability  of  the  drop,  and  were  studied  by  Bond  and  Newton 
(Reference  47),  These  investigators  showed  that  the  relevant  parameter 
is  the  Bond  number,  Bo.  which  may  be  interpreted  as  being  approximately 
the  ratio  of  the  hydrodynamic  head  (accelerational  or  gravitational)  pressure 
to  the  surface-tension  pressure  of  a  liquid. 

For  a  liquid  drop  accelerated  through  a  caseous  medium,  the  Bond  number 
may  be  written  as 


<Pt  *  Pa)  "  dv 

Bo  =  -  —  (34) 

<r  at 

where  dv/dt  is  the  acceleration  of  the  droplet  relative  to  its  gas  environ¬ 
ment,  Experiments  have  shown  (e.  g.  ,  Reference  12)  that  when  the  Bond 
number  of  an  accelerated  drop  exceeds  a  certain  number  (about  8  to  12, 
although  the  exact  value  also  appears  to  be  a  function  of  the  viscosity  of 
the  liquid),  the  drop  becomes  unstable  and  will  probably  break  up.  The 
time  required  for  breakup  to  occur,  however,  is  not  specified 
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It  may  also  be  recalled  from  the  previous  discussions  that  a  critical 
value  of  the  Weber  number  is  aisc  usually  required  for  droplet  instability 
and  breakup.  The  Weber  number  (Equation  15)  may  b^  interpreted  as 
being  approximately  the  ratio  of  the  drag  pressure  to  the  surface-tension 
pressure  of  a  liquid  drop. 


4.  THEORY  OF  THE  MECHANISM  AND  RATE  OF  BREAKUP  OF 
LIQUID  DROPS 


4.  1  INTRODUCTION 

A  common  assumption  often  made  in  discussing  the  breakup  of  liquids  is 
that  when  the  maximum  force  tending  to  disrupt  the  liquid  exceeds  the 
surface-tension  force  that  tends  to  hold  it  together,  the  liquid  will  burst, 
i.e..  Equation  13  is  often  assumed  to  define  the  breakup  criteria.  This 
assumption  was  used  in  the  classical  work  of  Taylor  (Reference  48),  and 
has  also  been  used  in  many  contemporary  treatments  of  the  breakup 
process  <e.  g. ,  Reference  49).  The  assumption,  however,  is  true  only 
for  small  rates  of  stress  loading,  and  is  not  valid  for  high  (e  g.,  shock) 
loading  rates,  since  the  flow  or  breakup  of  a  liquid  has  been  shown  experi¬ 
mentally  to  be  a  rate  process  Moreover,  in  any  system  in  which  the 
stress  tending  to  break  the  liquid  undergoes  a  change  in  a  time  less  than 
the  breakup  time  of  the  liquid,  it  will  be  expected  that  the  assumption  will 
also  be  erroneous. 

The  preceding  will  be  true  in  many  situations  involving  shock-loading  of  a 
liquid.  An  examination  of  the  various  theoretical  treatments  of  liquid 
breakup  summarized  in  the  previous  section  shows  that  only  hydrodynamic 
and  mechanic  approaches  have  been  pursued  Since  the  available  evidence 
suggests  that  the  breakup  of  a  liquid  is  a  rate  process,  the  breakup  process 
should  perhaps  be  amenable  to  kinetic  approaches  heretofore  not  used 
This  consideration  does  not  indica'e  that  the  hydrodynamics  and 
mechanics  of  the  problem  should  be  ignored  they  should  be  incorporated 
in  the  proper  kinetic  expression  of  the  problem  The  proper  kinetic 
model  should  be  capable  of  predicting  loth  the  breakup  times  of  the  liquid 
drop  and  the  initial  sizes  of  the  resultant  drops  (perhaps  without  recourse 
to  assumptions  previously  necessary  for  solution  of  theproblemi.  in 
ar.y  case,,  the  expressions  may  allow  a  deeper  elucidation  of  the  roles  of 
both  the  physical  properties  of  a  liquid  and  the  hydrodynamical  properties 
of  the  flow  to  the  breakup  cheracterisf.es  of  liquids 
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4.  2  LIQUID  FLOW  AND  BREAKUP  AS  RATE  PROCESSES 


Evring  (References  50  and  51)  was  one  of  the  first  to  conside*  flow  as  a 
rate  process.  He  reasoned  that  for  molecular  flow  to  take  place  in  a  con¬ 
densed  material  (solid  or  liquid),  it  is  necessary  that  a  suitable  "hole" 
or  site  be  available  for  a  molecule  to  "jump"  or  flow  into.  Production  of 
the  holes  was  pictured  as  being  brought  about  by  thermal  fluctuations  of 
the  molecular  energy,  i.  e.  ,  the  molecular  energy  distribution  continuously 
provides  a  fraction  of  the  molecules  with  enough  energy  to  push  their 
neighboring  molecules  aside  so  as  to  form  holes  into  which  other  molecules 
will  immediately  jump. 

An  alternative  picture  was  that  in  order  for  a  molecule  to  jump  to  a 
neighboring  position,  it  has  to  acquire  by  chance  a  sufficiently  high  kinetic 
energy  (i.  e. ,  it  must  be  activated),  which  enables  it  to  leave  the  potential 
sphere  of  its  immediate  neighbor  to  a  new  equilibrium  position.  In 
either  case,  an  (activation)  energy  of  sufficient  magnitude  is  required  to 
both  form  the  holes  or  site  of  the  new  position  and  move  the  molecule  into 
the  hole.  Eyring  argued  that  the  energy  required  to  form  a  hole  ot  mole¬ 
cular  size  in  a  liquid  is  equal  to  the  energy  oi  vaporization  per  molecule 
of  the  liquid.  He  pointed  out,  however,  that  a  distribution  of  hole  sizes 
will  always  be  available  as  a  result  of  thermal  motion,  and  that  it  may 
not  necessarily  require  a  hole  of  molecular  size  for  molecular  flow  to 
occur. 

The  absolute  specific  rate  constant,  kr,  for  any  thermal  rate  process  was 
derived  by  Eyring  (References  51  and  52)  in  terms  of  molecular  constants, 
and  is  given  by 

kx  F*  kT  *  1 

kr  -  K—  — •  exp  (-  AE*  /RT)  =  K—  exp  {- AF  /RT)  sec 
r  h  Fj  °  h 

where 

K  -  transmission  coefficient 

k  =  Boltzmann's  constant 

h  =  Planck's  constant 

T  =  temperature  (°K) 

F*  =  appropriate  partition  function  for  the  activated  molecular 
state 

Fj  =  appropriate  partition  function  for  the  nitial  molecular  state 


(35) 
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AE*  *  difference  in  the  zero-point  energies  of  the  initial  and  the 
activated  state,  i.e. ,  the  activation  energy  of  the  process 
at  the  absolute  zero  of  temperature  (0°K) 

R  =  gas  constant 

AF*  =  Gibbs'  free  energy  of  activation  when  the  initial  and  ac¬ 
tivated  molecules  are  all  in  their  standard  states 


Equation  35  gives  the  specific  rate  for  any  molecular  thermal  rate  process. 
The  reciprocal  of  the  specific  rate  constant,  i.e. ,  l/kr,  is  the  time  it 
takes  for  the  specified  molecular  process  to  occur., 

Eyring  showed  that  to  properly  consider  a  flow  process  in  a  fluid  as  a 
rate  process,  the  specific  rate  constant  (Equation  35)  must  be  interpreted 
as  the  frequency  (specific  rate)  with  which  molecules  jump  in  various 
directions  as  a  result  of  thermal  motion  in  the  fluid  when  there  is  no 
applied  stress.  The  influence  of  an  applied  shear  stress  on  the  specific 
rate  was  then  incorporated  by  multiplying  Equation  35  by  an  appropriate 
expression  involving  the  shear  stress.  The  resulting  specific  rate  is 


ks  =  2kr  sinh  ■ . —  Px..  sec 

2kT 


(36a) 


also 


1  L/Na 


(36b) 


where 

ks  =  absolute  specific  rate  constant  unde  application  of  the  stress 

PXy  =  applied  shear  stress 

L^Lj  =  effective  area  of  a  molecular  "flow  unit"  on  which  the  force 
is  acting,  i.  e.  ,  the  area  of  the  flow  unit  in  the  shear  plane 

L  =  average  jump  distance  of  a  molecular  flow  unit  during  the  flow, 
i.  e.  ,  the  distance  between  two  equilibrium  positions  in  the 
direction  of  motion 

Na  =  average  number  of  molecular -flow-unit  bombs  per  unit  area 
in  the  shear  plane  of  the  flow  unit 

kr  =  thermal  rate  constant  (Equati...:  35) 
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The  influence  of  the  applied  stress  is  to  lower  the  activation  energy 
barrier,  and  hence  the  resistance  to  flow  in  the  direction  in  which 
the  stress  is  acting. 

Liquid  flow  is  the  result  of  shear  strain,  and  the  net  velocity  of  flow 
in  the  shear  (axial)  direction  is 


v  =  Lks 


(37) 


The  rate  of  shear  of  the  liquid  *s  given  by 


dv  /dr  *  ds/dt  =  \  /Lj  *  Lkg/Lj 


(38) 


where  Lj  is  the  distance  between  the  shear  planes,  and  r  is  the  radial 
distance.  Thus,  Equation  38  indicates  that  the  shear  rate  is  equal  to 
the  net  number  of  molecular  jumps/sec  made  by  the  flow  unit  in  the 
shear  direction,  ks,  multiplied  by  the  average  distance  traversed  per 
jump,  L,  divided  by  distance  between  the  fluid  planes  undergoing  the 
shearing,  Lj. 

4.  3  BREAKUP  MECHANISM  AND  BREAKUP  TIME  OF 
LIQUID  DROPS 

When  a  liquid  drop  is  subjected  to  a  nonumform  pressure,  a  shearing 
of  the  liquid  will  occur  if  the  pressure  is  sufficient  to  overcome  the 
surface  tension  of  the  drop.  The  geometry  of  the  liquid  undergoing  the 
shearing  will  depend  on  the  (generally  time -dependent)  pressure  distri¬ 
bution  on  the  drop.  When  the  drop  is  subjected  to  an  aerodynamic  fluid 
(gas)  flow,  the  total  drag  force  exerted  on  the  drop  is  composed  of 
two  components,  viz. ,  the  pressure  drag  due  to  the  pressure  distribu¬ 
tion  over  the  surface  of  the  drop,  and  the  friction  drag  due  to  viscous 
shear  at  the  surface  of  the  drop  (Reference  53).  It  is  natural  to  believe 
that  these  two  individual  stresses  are  responsible  tor  the  two  (extreme) 
breakup  behaviors  of  liquid  drops,  and  in  particular  that  the  pressure 
drag  produces  bag  breakup,  while  the  friction  drag  is  responsible  for 
shear  or  stripping  breakup..  The  preceding  statement  appears  self- 
evident;  however,  the  breakup  mechanisms  of  liquid  drops  have  never 
been  expla-  .ed  in  this  manner  before.  In  the  following  discussion,  this 
concept  is  quantitatively  developed. 
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When  a  liquid  drop  is  subjected  to  an  aerodynamic  gas  flow,  the  stress 
produced  by  the  gas  flow  on  and  around  the  drop  initially  deforms  the 
drop.  It  is  known  from  photographic  studies  that  the  drop  initially 
flattens  on  the  leeward  side,  so  that  the  drop  diameter  increases  with 
increase  in  time,  and  then  remains  for  a  short  time  with  an  essentially 
constant  diameter  while  liquid  starts  to  be  removed  (sheared)  from  the 
drop  and  carried  in  the  gas  stream  (Reference  36)  The  drop  then 
starts  to  expand  again  this  time  very  slowly  with  time  while  breakup  »s 
occurring  it  is  known  experimentally  that  the  final  (essentially  stationary) 
drop  diameter  prior  to  breakup  is  only  a  weak  function  of  the  flow  condi¬ 
tions  and  the  properties  of  the  iiqu.d  Although  a  knowledge  of  the  exact 
shape  and  size  of  the  flattened  drop  before  shearing  is  not  essentia!  to 
the  present  studies,  it  is  of  interest  to  show  that  the  rate  and  exent  of 
flattening  is  a  function  of  the  Weber  number,  and  includes  a  mechanism 
whereby  the  lateral  flattening  or  spreading  of  the  drop  is  impeded  just 
prior  to  the  breakup  of  the  drop,  thus  allowing  the  gas-f'ow  stress  to  alter 
the  direction  of  the  iquid  flow  with  subsequent  breakup  ol  the  liquid 
This  initial  detormation  of  the  drop  a-so  appears  to  have  an  influence  on 
the  general  mechinism  (bag  or  shear'  bv  which  the  drop  is  dispersed 

The  pressure  causing  the  mitiai  drop  deformation  is  the  air  stagnation 
pressure,  which  is  opposed  by  surface  tension  and  viscous  stresses  of 
the  deforming  drop.  The  model  used  by  Gordon  to  describe  drop  breakup 
(Reference  39;  can  be  used  to  describe  the  drop  deformation  if  the  model 
is  modified  slightly  Stagnation  pressure  prod’  ces  the  lateral  flattening 
of  the  drop,  and  hence  this  pressure  acts  through  a  :aterai  shear  stress 
The  surface-tension  pressure  retarding  the  elongation  arises  from  the 
shape  of  the  drop  at  any  arbitrary  time  but  may  be  considered  to  origi¬ 
nate  mainly  from  the  head  of  the  .atera.iy  e  engating  liquid  ;n  the  latter 
phases  of  elongation.  The  shear  stress  causing  the  deformation  will  be 
roughly  one-half  of  the  norma,  stagnation  pressure 


Velocity  of  the  laterally  flowing  .;qu:d  idrop  deformation  rate,  is  hence 
given  by  Equation  23b  with  a  factor  of  4  replacing  2  in  the  A  term  and 
a  new  (unknown)  factor  K.  rep.acing  8/d  in  the  surface-tension  expres¬ 
sion  of  the  A  term  For  both  simplicity  and  purposes  of  illustration,  a 
liquid  is  considered  wi’h  a  low  viscosity,  so  that  the  exponential  term  of 
Equation  23b  may  be  exp-mded  with  the  retention  of  on.y  the  linear  term 
Integrating  Equation  23b  (assuming  a  cor.stai"  averaged  K  K't  and  con¬ 
sidering  the  drop  radius  to  be  that  of  the  or.gmat  spherical  drop  R,  when 
time  is  zero 


a 


4  f  d 


Kc 

3  ;  ci 


■> 
t  “ 


(39. 
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where  a  is  the  radius  of  the  flattening  drop  at  time  t.  Thus,  the  rate  of 
lateral  deformation  of  the  drop  is  proportional  to  the  difference  of  two 
terms  that  involve  the  stagnation  pressure  deforming  the  drop,  and  the 
surface-tension  pressure  resisting  the  deformation.  For  high  gar  flow 
or  low  surface  tension,  the  surface-tension  pressure  will  originally  be 
small  as  compared  with  the  stagnation  pressure,  in  which  case  the 
deformation  rate  is  proportional  to  the  stagnation  pressure  and  the  square 
of  the  time,  and  inversely  proportional  tr  the  liquid  density  and  the  drop 
diameter. 

The  preceding  is,  interestingly,  the  same  functional  dependence  found 
by  Burgers  (Reference  36)  for  the  deformation  rate  of  the  deforming  drop. 
In  fact,  the  pressure  term  in  Equation  39  is  identical  to  the  expression 
developed  by  Burger,  except  that  the  Burger  expression  contains  a  factor 
of  3  in  the  denominator,  whereas  Equation  39  contains  a  factor  of  2  (after 
the  other  factor  of  2  is  included  in  the  stagnation  pressure). 

While  the  liquid  drop  is  being  flattened  by  the  stagnation  pressure,  the 
curvature  of  the  front  of  the  deforming  liquid  increases.  This  change 
in  curvature  increases  the  effective  surface-tension  pressure  resisting 
the  drop  deformation.  When  the  curvature  is  sufficiently  large,  the  drop 
ceases  to  elongate,  and  the  aerodynamic  stress  then  changes  the  lateral 
liquid  flow  to  the  horizontal  direction,  leading  to  shear  and  breakup  of 
the  liquid. 

The  value  of  K  in  Equation  39  should  thus  be  replaced  by  a  function  of 
the  deformation  so  that  the  surface-tension  pressure  increases  with  an 
increase  in  deformation..  More  accurately,  this  function  should  be 
included  in  Equation  23b  before  it  is  integrated  to  give  Equation  39.  The 
nature  of  the  function  depends  on  the  geometry  that  the  drop  assumes 
during  its  deformation..  This  geometry  is  not  known  in  detail,  but  an 
inspection  of  the  available  photographs  on  the  behavior  of  liquid  drops 
in  gas  flows  suggests  that  the  drops  initially  assume  approximately  the 
shape  of  one -half  of  an  oblate  spheriod  (i.  e  ,  a  lateral  half  of  an  ellipse 
rotated  about  its  minor  axis)  before  undergoing  their  final  shearing 
leading  to  breakup  If  this  sequence  is  roughly  true,  then  the  surface- 
tension  pressure  of  the  outer  end  of  the  oblate  spheriod  must  support  the 
aerodynamic  flow  pressure  before  the  final  shearing  leading  to  breakup 
occurs.. 

Since  the  maximum  radius  of  curvature  of  an  oblate  spheriod  it  R^/a®, 
the  value  of  K  in  Equation  23b  can  be  taken  roughly  to  be  a^/P^.  The 
(essentiall-  stationary)  drop  diameter  just  prior  to  breakup  is  then  given 
by  the  solution  of  Equation  23b  for  v  =  da/dt  =  0,  i.e.  .  for  1/2  pau^  = 
o-a^/R^.  Therefore 


(40) 
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where  the  Weber  number  is  defined  by  Equation  15  This  equation 
illustrates  the  insensitive  nature  of  the  deformed  drop  size  on  the  gas- 
flow  conditions  and  the  phv  ical  properties  of  the  liquid  drop  Although 
Equation  40  gives  only  a  rough  quantitative  representation  of  the  experi¬ 
mental  data,  its  overall  general  agreement  with  experiment  is  reasonable 
in  view  of  the  assumptions  made,  and  indicates  that  a  lateral-halt  oblate 
sphenod  is  not  a  bad  approximation  to  the  true  shape  of  the  deformed 
drop,  While  breakup  of  the  drop  is  occurring,  the  drop  continues  to 
expand  slowly. 

For  bag  breakup,  the  formation  and  expansion  of  the  bag  aids  the  drop- 
deformation  process;  for  stripping  breakup,  the  shearing  i.quid  helps 
weaken  the  surface-tension  forces  that  prevent  the  drop  from  expanding 

Equation  39,  as  derived,  is  for  low -viscosity  liquids  It  is  known  that 
the  rate  of  drop  deformation  also  depends  on  viscosity  when  the  viscosity 
is  large  (References  30  and  31),  The  effect  of  viscosity  is  easily  included 
by  retaining  the  exponential  term  in  the  integration  of  Equation  23b 

Subsequent  to  the  initial  deforming  of  the  drop,  the  l.quid  is  either  sheared 
away  from  the  sides  of  the  drop  (stripping  breakups  or  the  middle  of  the 
drop  is  sheared  away  from  the  sides  (bag  breakup)  This  shearing  pro¬ 
duces  a  stretching  of  the  liquid  to  form  a  film  which  ultimately  breaks 
to  form  droplets.  The  breaking  of  the  film  is  a  very  fast  process,  and  it 
is  the  rate  of  formation  of  this  film  by  shear  that  thus  controls  the  breakup 
time  of  the  drop. 


The  rate  of  liquid  shear  is  given  by  Equation  38.  H-gh  shear  stresses 
cannot  be  built  up  in  liquids  unless  the  liquid  is  extremely  viscous,  and 
he">ce  the  hyperbob<~  sine  quentity  in  Equation  36a  will  always  bt  a  small 
number.  Since  smli  •  *  for  small  x,.  sinh  ,LP  /NakT*  •  LPxy/MakT. 
Therefore 


dv 

dr 


L2  k. 


L1  Na  kT 


P 


xy 


(41, 


The  rate  of  shear  should  be  controlled  by  self -diffusion  of  the  molecules 
of  the  liquid,,  since  a  molecular  flow  >s  involved  in  the  formation  of  the 
film  from  the  bulk  liquid  of  the  drop  P  this  assumption  is  true,  then  the 
coefficient  of  PXy  in  Equation  41  can  be-  shown  to  be  identically  equal  to 
tne  reciprocal  of  the  coefficient  of  viscosity  of  the  liquid  'References  50 
and  51).  Hence 


dv/  Ir  ^  /  ri 

x> 


(42) 
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Equation  42  may  be  applied  to  the  breakup  of  liquid  drops  as  follows: 

Consider  a  cylindrical  tube  of  liquid  of  length  d  and  radius  R  to  be 
shearing  from  a  drop  under  the  influence  of  an  average  pressure  head 
Pjj.  At  any  point  whose  distance  from  the  axis  of  the  tube  is  r,  the  shear 
stress  is  related  to  the  pressure  head  by  itr^Pjj=  2irrdPXy.  Hence 

Pxy  =  Phr/2d  (43) 

Substituting  Equation  43  in  Equation  42  and  integrating  under  the  condition 
that  v  =  0  when  r  =  R  (assuming  Uv/dr  is  negative)  gives 


v  =  Ph  (R2  -  r2)/4d  q 


(44) 


The  maximum  velocity  occurs  at  r  =  0,  and  it  is  related  to  the  pressure 
head  by 


Ph  =  4dqv/R2 


(45) 


The  pressure  head  is  the  difference  between  the  aerodynamic  stagnation 
pressure  and  the  pressure  required  to  provide  the  kinetic  energy  of  the 
liquid  flow.  Thus  (Reference  54) 


P 


h 


P  -  l/2pjV 


2 


(46) 


Substituting  Equation  46  in  Equation  45 


4d  r)  v 

- 5 — 

R  ' 


+  l/2p(v2 


(47) 


where  the  velocity  gives  the  time  rate  of  displacement  of  the  liquid  being 
sheared  from  the  drop.  To  obtain  the  proper  order  of  magnitude  for  the 
breakup  time  of  the  liquid  being  sheared,  it  is  assumed  that  K  is  the 
radius  of  ne  original  liquid  drop,  and  d  is  the  drop  diameter.  Integrating 
Equation  47  and  assuming  breakup  to  occur  when  the  length  oi  the  tube 
being  sheared  from  the  drop  is  2d  gives  the  bre- ",«,«»  time 
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d  _ 

(A2,  BP)~~A 


(48; 


where 


A  =  16  r|  /ripj  ;  B  =  2 / f.  j 


(48b 


Equation  48a  is  a  completely  general  expression  Tor  the  breakup  time' 
it  may  be  compared  to  the  general  breakup-time  expression  of 
Gordon  (Equation  24).  The  pressurt  expression  will  be  o*  the  form 


P  =  1/2  p  u2C  -  k  (r/d 

a  D 


(49) 


where 

n  =  relative  velocity  between  the  air  stream  and  the  liquid 
drop 

Pa  =  density  of  the  air  stream 

Cp  -  drag  coefficient 

a  s  surface  tension  of  the  liquid  drop 

d  -  diameter  of  the  original  liquid  drop  before  breakup 

k  -  a  constant  that  reflects  the  diop  curvature  during 

breakup  (which  determines  the  effective  surface-tension 
pressure  retarding  the  breakup) 


Equation  48a  may  be  greatly  simplified  for  certain  conditions  "lhus, 
for  liquids  or  flow  conditions  for  -.'hich  viscous  and  surface-tension 
forces  are  negligible  (and  for  convenience  assuming  Cjy  to  be  unity) 
Equation  47  becomes 


d  1/2 

7  !rl  ' 


I  50a ) 


This  expression  is  very  similar  to  that  nd  by  Gordon  (Equation  25a t, 
and  also  by  Hinze  (Equation  17}  for  sinnl  .  r  conditions  differing  essen¬ 
tially  only  in  the  value  of  the  numerical  co-.stant  preceding  the  function 
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(unity  in  Equation  50a).  For  extremely  viscous  liquids  and  negligible 
surface  tension.  Equation  48a  becomes 


=  32  rj  /  p^u 


(50b) 


This  expression  is  identical  to  that  of  Gordon  (Equation  25b),  and  differs 
from  that  of  Hinze  (Equation  18)  only  in  the  numerical  constant. 

Equation  48a,  as  derived,  describes  the  shearing  of  a  plug  of  liquid  from 
a  drop  as  a  result  of  pressure  being  exerted  on  the  middle  (thus  excluding 
the  edges)  of  the  drop.  Hence,  this  equation  will  apply  to  the  bag  break¬ 
up  of  a  drop.  The  same  general  expression,  however,  also  applies  to 
the  shear  (stripping)  breakup,  as  may  be  shown  by  assuming  dv/dr  in 
Equation  42  to  be  positive  and  integrating  under  the  condition  that  the 
velocity  is  zero  when  r  is  zero.  This  condition  thus  provides  for  the 
situation  that  the  maximum  pressure  is  exerted  on  the  edges  of  the  drop, 
as  occurs  in  stripping  breakup.  Thus,  there  remains  only  to  identify  the 
pressure  distribution  and  magnitude,  using  Equation  48a  to  delineate 
between  the  bag  and  stripping  breakup  mechanisms.  Before  considering 
these  individual  breakup  mechanisms,  however,  it  is  necessary  to  examine 
certain  aspects  of  Equations  48a  and  4<>. 

Equation  48a  is  a  generalized  expression  for  the  breakup  time  of  a  liquid 
drop  containing  given  physical  properties  and  subjected  to  an  aerodynamic 
flow  with  given  properties.  This  equation  is  expected  to  be  generally 
valid  except  whe»  the  liquid  is  extremely  viscous,  i.,  e.  ,  when  the  liquid 
is  a  solid.  Under  the  latter  condition,  it  is  not  permissible  to  replace 
the  hyperbolic  sine  function  of  Equation  36a  by  the  function  itself,  since 
very  large  shear  stresses  can  be  built  up  within  the  solid.  In  this  case. 
Equation  3fea  must  be  used  as  written  in  integrating  Equation  41  (final 
results  are  given  in  Reference  55),  and  the  kinetic  energy  term  in  Equa¬ 
tion  46  will  usually  be  negligible.  Thus,  th-'  presented  theory  for  the 
breakup  of  liquids  can  also  be  used,  with  slight  modification,  to  estimate 
the  mechanical  erosion  (breakup)  of  solid  particles  in  aerodynamic  flows. 

It  is  possible  to  use  Equation  48a  to  predict  the  breakup  time  of  a  liquid 
drop  without  regard  to  the  general  mechanism  (bag  or  stripping)  if  a 
suitable  value  is  chosen  for  k  in  Equation  49.  As  will  be  shown  later, 
Equation  48a  is  in  very  good  agreement  with  all  of  the  experimental  data 
if  the  drag  coefficient  is  considered  to  be  unity,  and  the  value  A  k  is  taken 
as  about  tw  j.  This  equation  thus  affords  a  rapid  calculation  of  the  breakup 
time  of  a  liquid  drop  to  the  accuracy  with  which  it  can  be  measured. 
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Photographs  of  the  breakup  process  (provided  in  a  subsequent  section 
of  this  report)  show  that  a  finite  time  exists  between  the  time  when  the 
liquid  starts  to  break  and  the  time  when  breakup  is  es^t  :i- ial!y  complete 
This  time  during  breakup  is  often  comparable  to  the  time  required  for 
the  breakup  to  begin  to  occur  Thus,  there  is  some  latitude  in  defining 
an  experimental  breakup  time  for  comparison  with  theory 

For  consistency  and  general  accuracy,  the  experimental  breakup  time 
has  been  taken  to  be  the  time  between  the  aerodynamic  flow  hitting  the 
drop  and  the  time  at  which  breakup  of  the  drop  just  starts  to  occur  The 
total  breakup  time,  however,  will  usually  be  slightly  larger,  and  it  is 
thus  natural  that  the  theoretical  breakup  time  may  be  slightly  larger  than 
the  experimental  breakup  time  The  latitude  allowed  in  comparing  experi¬ 
ment  and  theory  makes  it  difficult  to  check  with  great  accuracy  certain 
points  of  the  theory,  such  as  the  validity  of  the  breakup  criteria,  as  well 
as  the  laminar -flow  behavior  implied  by  Equation  42  including  *’  eglect 
of  end-effects  and  nonaccelerated  flow  (Reference  5fe>  This  latter 
assumption  may  not  be  quite  true,  but  photographs  of  the  breakup  process 
suggest  the  acceleration  to  be  very  small,  if  any  (after  the  initial  liquid - 
flow  phase)  The  assumption  of  an  accelerated  flow  is  implicit  in  the 
treatments  of  Gord">n  and  of  Hinze 

The  equations  leading  to  Equation  48a  have  used  the  maximum  liquid 
velocity,  rather  than  the  average  velocity,  produced  by  a  pressure  on  the 
liquid.  This  usage  is  believed  to  be  essentially  correct,  since  employing 
an  average  velocity  implies  that  a  parabolic  velocity  distribution  has 
been  established  in  the  flowing  liquid  However,  the  distances  required  to 
establish  this  distribution  are  very  large  as  compared  to  (he  size  of  a 
liquid  drop  (Reference  54j 

The  effects  of  using  the  average  velocity  are  easily  estimated  Average 
velocity  of  liquid  ejection  is  given  by  va  -  J  vr  dr/J  r  dr  Hence, 
va  -  R^P/8d  q  If  the  average  velocity  13  used,  then  the  total  kinetic 
energy  of  flow  is  twice  that  used  fo-  the  maximum  velocity  flow  because 
a  parabolic  velocity  distribution  is  implied  (Reference  54i  Upon  inte¬ 
grating  the  equations  it  is  found  that  if  an  average  velocity  is  involved, 
then  a  factor  of  2*^  and  2  should  respectively  be  included  in  Equations 
50a  and  50b,  and  the  factor  of  2  should  be  removed  from  the  B  term  in 
Equation  48b  These  factors  are  generally  smaller  than  the  uncertainty 
in  measuring  the  breakup  time 

Another  consideration  is  the  value  of  the  drag  coeffit,ent  to  be  used  ir. 
Equation  49  An  inspection  of  the  available  literature  on  the  expsrimental 
drag  coefficients  of  various  geometry  bodies  mciud'iig  solid  spheres  "nd 
deformed  liquids,  shows  that  a  rather  «•  range  of  values  has  been 


found  {References  36,  38,  53,  54,  and  57  through  63).-  Many  of  the 
values  reported  are  for  steady  flow,  whereas  accelerational  drag  forces 
may  also  be  of  importance.  During  the  breakup  of  a  liquid  droplet, 
however,  the  shape  of  the  drop  changes  continuously  throughout  the 
breakup;  hence,  only  some  suitable  average  drag  coefficiert  can  properly 
be  used  in  the  breakup-time  expression. 

In  a  later  section  of  this  report,  experimental  illustrations  are  provided 
shoeing  the  change  of  the  drag-coefficient  value  during  the  breakup 
process.  The  value  of  the  drag  coefficient  changes  from  less  than  unity 
to  greater  than  unity  during  the  breakup  process,  and  for  most  practical 
work  it  seems  sufficient  ,o  consider  the  drag  coefficient  to  be  unity. 

Rabin  et  al. (Reference  38)  actually  found  the  value  to  be  essentially  unity 
over  a  large  range  of  conditions.  The  preceding  shows  that  if  the  value 
k  in  Equation  49  is  chosen  properly,  then  Equation  48a  is  an  adequate 
representation  of  all  of  the  present  experimental  data  if  the  drag  coefficient 
is  taken  to  be  unity.  For  most  practical  work,  therefore,  this  representa¬ 
tion  appears  adequate. 

As  previously  discussed,  it  is  known  from  experiments  that  the  breakup 
of  a  liquid  drop  may  occur  from  either  of  two  (extreme)  mechanisms, 
viz. ,  bag  or  stripping  breakup.  Ar.  adequate  model  of  liquid-drop  breakup 
should  provide  an  explanation  of  tnese  two  mechanisms;  the  explanation 
proposed  here  is  that  bag  breakup  is  the  result  of  pressure  drag,  whereas 
stripping  breakup  is  the  result  of  friction  drag.  Thus,  two  individual 
expressions  of  the  form  of  Equation  48a  may  be  written  for  the  breakup 
time  of  a  liquid  drop,  with  one  expression  containing  the  pressure-drag 
stress  in  the  pressure  expression  (Equation  48b),  and  the  second  expres- 
containing  the  friction  drag  stress  :n  the  pressure  expression. 

Breakup  of  a  liquid  drop  with  given  physical  properties  by  an  aerodynamic 
gas  flow  with  given  properties  will  occur  by  the  mechanism  that  takes 
place  with  the  fastest  rate,  i.  e. ,  with  the  shortest  breakup  time.  When 
the  rates  are  comparable,  the  drop  should  exhibit  both  bag  and  stripping 
breakup  behaviors. 

A  quantitative  description  of  the  dual  mechanism  proposed  for  liquid-drop 
breakup  requires  a  knowledge  of  the  average  (effective)  stress  distribution 
on  the  drop  undergoing  breakup.  The  details  of  this  distribution  are  not 
fully  understood,  but  a  theoretical  study  by  Tomotika  and  Aoi  (Reference 
64)  indicated  that  for  a  sphere,  the  frictional  drag  is  twice  the  pressure 
drag  for  all  Reynolds  numbers.  It  is  assumed  that  this  explanation  is 
also  t-ue  for  liquid  drops  undergoing  breakup,  even  though  the  shape 
deviates  from  that  of  a  true  sphere.  The  total  drag  stress  on  the  diop 
during  breakup  is  given  by  1/2  pa  u**  Cp,  and  hence  fc*-  bag  and  stripping 
breakup,  the  pressure  expression  given  by  c.  ation  49  may  be  written. 
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where  k^  and  kg  are  constants  that  reflect  the  surface -tension  pres¬ 
sure  tending  to  hold  the  drop  together  during  bag  and  stripping  breakup 

If  the  fric*ional  drag  is  always  greater  than  the  pressure  drag,  as  indi¬ 
cated  by  Tomotika  and  Aoi,  then  it  may  be  imagined  that  liquid  drops 
should  always  undergo  sLnpping  breakup  Since  the  preceding  is  not  true,, 
it  is  indicated  that  the  surface-tension  pressure  holding  the  drop  together 
at  its  edges  is  greater  than  in  its  middle,  i.e  ,  that  kg  >  In  the  early 
part  of  this  section,  however,  it  has  already  been  shown  that  the  initial 
drop  deformation  produces  a  greater  curvature  at  the  outer  edge  of  the 
drop  than  m  its  middle,  and  hence  the  surface-tension  pressure  holding 
the  outer  edges  of  the  drop  is  greater  than  in  the  middle  It  is  thus 
known  semiquantitatively  that  k3  must  be  greater  than  kj, 

A  calculation  of  k^  and  ks  would  require  a  detailed  knowledge  (which  is 
unavailable)  of  the  effective  shape  of  the  liquid  drop  during  its  breakup 
However,  a  theoretical  calculation  of  these  constants  is  presently  of 
little  value,  since  the  true  stress  distribution  on  the  drop  undergoing 
breakup  is  not  known,  but  rather  assumed  to  be  that  of  a  solid  sphere 
The  values  of  these  constants  are  later  determined  by  a  best  fit  of  Equa  - 
tions  48a,  51a,  and  51b  to  the  experimental  breakup  data  In  anticipation 
of  these  results,  respective  values  of  about  4  and  2  were  found  for  kj, 
and  ks;  these  values  may  be  employed  in  using  the  preceding  equations  to 
predict  breakup  times  for  drops  undergoing  either  bag  or  stripping  breakup 

A  consequence  of  the  dual-breakup  mechanism  given  by  Equations  51a  and 
51b  is  that  at  low  air-flow  velocities,  the  drop  will  usually  (depending  on 
drop  size)  undergo  bag  breakup  because  the  drop-edge  surface-tension 
forces  prevent  stripping  breakup,  whereas  as  the  a.r  velocity  is  further 
increased,  the  stripping  breakup  becomes  faster  as  a  result  of  the  overall 
greater  friction  stress.  The  equations  =*lso  predict  the  influence  of  the 
v  rious  other  parameters  on  the  breakup  mechanism  The  predicted 
behaviors  can  be  observed  experimentally 
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4.  4  INITIAL  DROP  SIZES  PRODUCED  BY  DISPERSION  OF 
LIQUID  DROPS 

Two  considerations  are  usually  of  crime  interest  in  the  aerodynamic 
breakup  of  liquid  drops,  viz. ,  the  time  required  to  produce  the  break¬ 
up  as  a  function  of  the  air-flow  properties  and  the  physical  properties 
of  the  drop,  and  the  sizes  of  the  resultant  droplets  produced  from  the 
dispersion  of  the  original  drop.  A  quantitative  description  of  the 
breakup  process  was  developed  in  the  previous  section;  the  following 
discussion  treats  the  initial  mean  drop  size  produced  by  the  (primary) 
breakup  of  the  original  drop. 

The  drop  size  produced  by  primary  breakup  may  differ  from  the  mean 
drop  size,  which  may  be  measured  at  any  arbitrary  time  after  the 
breakup  of  the  drop.  Thus,  in  addition  to  the  drop  sizes  produced  by 
the  primary  breakup  process,  the  drop  sizes  present  at  any  arbitrary 
time  will  depend  on  the  following:' 

a.  Secondary  breakup  of  drops  produced  by  the  primary  breakup 
(which  in  turn  depends  on  the  properties  of  the  gas  flow,  e.  g. , 
its  duration). 

b.  Vaporization  of  the  primary  and  secondary  drops. 

c.  Coalescence  of  the  primary  and  secondary  drops. 

d.  Settling  or  removal  of  drops  by  winds,  etc. 

The  secondary  breakup  of  primary  drops  can  be  estimated  froi..  con¬ 
siderations  involved  in  ihe  primary  breakup  process  if  the  gas-flow 
properties  are  known,  and  the  remainder  of  the  factors  may  oft-  ,  he 
estimated  from  present  knowledge.  Therefore,  this  discussion 
confined  to  the  drop  sizes  produced  by  the  primary  breakup. 

As  discussed  in  the  previous  section,  the  aerodynamic  breakup  o’. 
drops  may  come  about  through  either  the  bag  or  stripping  mecha¬ 
nisms.  Both  mechanisms  are  believed  t«  be  the  result  of  a  shearing 
of  the  liquid  drop  by  a  nonunform  pressure  distribution  produced  by  the 
aerodynamic  flow  on  the  drop.  This  shearing  produces  a  liquid  film, 
which  subsequently  breaks  to  give  the  resultant  drops.  The  breakup 
time  of  the  drop  is  controlled  by  the  rate  of  formation  of  this  liquid 
film,  since  tha  actual  breaking  of  the  film  is  a  fast  process;  however, 
the  geometry  of  this  film,  together  with  the  manner  in  which  the  film 
breaks,  are  the  determinants  of  the  resultant  droplet  oizes. 
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Although  it  would  appear  that  the  drop  sizes  produced  by  the  two 
breakup  mechanisms  should  differ,  the  present  model  suggests  that 
they  obey  essentially  some  mathematical  size-distr  uution  law. 
Moreover,  experimental  data  obtained  in  the  present  susies  (to  be 
presented  later)  suggest  that  the  mean  diameter  of  the  resultant  drop¬ 
lets  produced  by  the  primary  breakup  of  liquid  drops  from  both  oag  and 
•stripping  mechanisms  obey  (within  experimental  error(  the  same  mathe¬ 
matical  expression  for  the  resultant  drop  size  as  a  function  of  gas  How 
and  liquid  properties.-  This  suggestion  also  appears  consistent  with  the 
observations  that  a  liquid  film  from  which  the  droplets  are  produced  is 
formed  in  both  mechanisms,  and  that  the  breakup  time  is  essentially 
monotomcally  continuous  in  the  transition  between  the  two  mechanism's 
breakup  behaviors.  For  the  present  study,  therefore  r.o  distinction  is 
made  between  the  two  breakup  mechanisms  in  treating  the  resulting 
droplet  sizes,  and  it  appears  that  more  experimental  drop-size  data 
are  required  to  warrant  distinguishing  the  two  mechanisms  m  the  future 

Another  consideration  concerns  the  manner  of  describing  the  droplet 
sizes  produced  by  the  primary  breakup  procesj  it  is  experimentally 
known  that  a  distribution  of  droplet  sizes  is  produced,  and,  ideally,  a 
theory  should  provide  a  function  that  predicts  the  droplet -size  distribu¬ 
tion  as  a  function  of  the  gas -flow  and  liquid  drop  properties  Although 
it  was  not  possible  to  attain  this  goal  during  thi  present  studies  experi¬ 
mental  examples  of  the  distribution  are  included  in  a  later  sect. on 

In  most  theories  of  liquid  breakup,  including  the  treatment  to  follow,  it 
is  assumed  that  breakup  is  brought  about  by  instabilities  'vibrations; 
that  grow  (or  decayi  with  time.  The  (optimum)  wavelength  that  grows 
the  fastest  leads  to  breakup  of  the  liquid,  with  the  subsequent  formation 
of  droplets.  This  wavelength  defines  a  mean  or  most  probable  droplet 
size,  and  it  is  this  average  that  is  computed  here  However  it  may  be 
noted  that  a  distribution  of  instability  wavelengths  is  always  present, 
and  these  values  lead  to  a  distribution  of  droplet  sizes  A  knowledge  of 
the  mean  size  allows  the  remain-ng  size  distribution  to  be  eat. mated 
(References  1  through  4  and  65).- 

The  following  discussion  presents  the  general  model  considered  to 
describe  droplet  formation  during  the  primary  breakup  of  liquid  drops. 
This  model  is  based  primarily  on  the  stripping  breakup  of  a  drop,  but 
essentially  the  same  view  may  be  taken  in  describing  the  droplet  siz'»s 
from  bag  breakup.  The  shear  stress  at  the  outer  surface  tor  middlei 
of  the  drop  produces  a  boundary  layer  of  flowing  liquid  of  average  thi<  k*. 

■  iss  6.  This  liquid  layer  (sheet)  upon  leaving  the  drop  experiences  a 
small  divergence  (expansion),  which  breaks  the  sheet  into  individual 
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strips  of  average  width  W.  Simultaneous  with  the  election  of  the  liquid 
boundary  layer  from  the  drop  and  its  divergence!  the  parallel  gas  flow 
along  the  surfaces  of  the  layer  induces  instabilities  (vibratu  r.s)  in  the 
layer  in  the  lengthwise  direction.  The  optimum  wavelength  instability, 
X ,  grows  most  rapidly  and  breaks  the  layer  into  pieces  of  average  size 
W6  X..  Surface-tension  forces  then  produce  spherical  drops  with  an 
average  diameter,  D,  from  these  pieces;  thus,  it  D ' / 6  -  W  6  X,  or 


(52) 


Photographs  of  the  breakup  process  roughly  support  the  proposed  model. 

The  value  of  X  has  been  calculated  from  the  hydrodynamic  equations  by 
Squire  (Reference  23),  who  found 
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where 


<r  =  surface  tension  of  the  liquid 

pa  =  air  density 

u  =  air  velocity 


The  value  of  6  follows  from  the  breakup  model  given  in  the  previous 
section.  The  shearing  of  the  liquid  follows  Equation  42.  Breakup  of 
the  liquid  film  is  assumed  to  occur  when  the  sheared  film  traverses  a 
distance  equal  to  the  drop  diameter,  d,  producing  a  triangular -shaped 
film  thickness  that  varies  from  zero  to  6,  or  an  average  thickness  of 
6/2.  Thus,  Equation  42  becomes 
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(54) 


For  simplicity,  liquids  are  considered  whose  viscosity  and  jurface- 
tension  forces  are  small  as  compared  with  the  aerodynamic  flow  forces 
involved.  This  consideration  includes  most  pra.  .ical  situations  of 


-38- 


interest.  The  v;tlue  of  t  is  then  given  by  Equation  50a,  and  the  shear 

stress  may  be  considered  'o  be  roughly  the  stagnation  pressure, 

1/2 p  u^.  Hence 
a 

1/2 

6  =  /u  (  papt)  (55) 


The  value  of  W  may  be  estimated  as  follows:'  Consider  an  arc  of  length 
S  in  the  plane  of  the  liquid  sheet  leaving  the  liquid  drop.  F  rom  geometry, 

S  =  rO.  where  r  is  the  distance  to  a  point  defining  the  radius  of  curva¬ 
ture  of  the  arc,  and  0  is  the  angle  subtended  by  the  arc.  Consider  the 
arc  to  undergo  expansion  in  time  t  as  a  result  of  the  divergence  of  the 
film;  thus,  ds/dt  =  0  dr/dt.-  The  derivative  dr/dt  is  the  velocity  of  the 
liquid  sheet,  v,  whereas  dS/dt  defines  a  critical  velocity,  vc,  for 
breakup  of  the  sheet  into  strips.  Thus,  vc  =  8v. 

The  distance  around  the  drop  comprises  2it  radians,  and  hence  the 
number  of  strips  N  =  ?e/0  =  2nv/vc;  the  numerical  distance  around  the 
drop  is  ird,  and  hence  W  =  rd/N  =  dvc/2v.  For  liquids,  as  in  the  pre¬ 
ceding,  having  low  viscous  and  surface-tension  forces  as  compared  v.  .th 
the  aerodynamic  flow  forces  pa  u^  =  p^v^  (^r°in  Equation  47),  an 
expression  that  defines  v  in  the  W  expression.  The  value  of  vc  may  be 
obtained  from  this  expression  and  the  conventional  equation  defining  the 
critical  aerodynamic  velocity  required  for  liquid  breakup,  i.e.  , 

1/2  Pau2CD  =  4  <t / d  (56) 


F'or  simplicity,  Cq  is  considered  to  be  unity.  Neglecting  the  small 
variation  of  air  density  with  air  velocity 


W  = 


a  v/d 


1/2 


Combining  Equations  52,  53,  55,  and  57  gives  the  average  (mass  mean) 
diameter  of  the  droplets  produced  by  the  aerodynamic  breakup  of  a  liquid 
drop 
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Equation  58  was  derived  for  conditions  where  the  aerodynamic  flow 
forces  are  much  larger  than  the  viscous  and  surface-tension  forces 
(this  situation  comprises  most  cases  of  practical  interes i,  If  the  liquid 
viscosity  is  very  great,  or  if  the  flow  velocity  is  low  and  the  surface 
tension  is  high,  then  Equation  58  as  written  may  not  be  valid,  !n  such 
cases,  the  various  approximations  indicated  in  the  derivations  should  be 
replaced  by  the  more  exact  expressions.  However,  it  is  possible  that 
Equation  58  may  also  hold  approximately  for  these  cases,  since  various 
approximations  and  errors  often  have  a  habit  of  canceling  themselves 
partially  out  in  a  multiparameter  expression.  In  any  event,  the  cube- 
root  dependence  of  Equation  58  attenuates  large  changes 

It  is  possible  that  the  factor  136  in  Equation  58  should  be  modified, 
since  various  approximations  and  assumptions  are  manifested  m  its 
value.  These  cor.siderutic-;.,  include  the  shear  length  required  for  breakup, 
the  value  of  drag  coefficient  used,  and  the  relation  of  shear  stress  to 
stagnation  pressure..  The  cube-root  dependence,,  however,  again  makes 
such  changes  seconaary. 

Equation  58  has  two  points  of  special  significance.  The  first  is  the  pre¬ 
diction  that  the  average  droplet  size  produced  varies  inversely  as  the 
4/3  power  of  the  relative  velocity  between  the  gas  stream  and  the  liquid. 
The  careful  work  of  Weiss  and  Worsham  (Reference  19)  actually  disclosed 
such  a  relationship,  which  is  further  expressed  by  the  theory  of  Mayer 
(Reference  46).-  Evidence  is  thus  accumulating  for  the  support  of  this 
velocity-droplet  diameter  relationship. 

The  second  ooint  of  interest  in  Equation  58  is  the  prediction  that  the 
droplet  diameter  should  vary  directly  as  the  1/6  power  of  the  initial 
droplet  diameter.  (This  relationship  was  found  experimentally  for  liquid 
jets  by  Weiss  and  Worsham.  )  Equation  58  gives  good  agreement  with  the 
uuia  of  Weiss  and  Worsham,  as  quoted  by  Mayer  It  will  later  be  shown 
that  this  equation  is  in  good  agreement  with  all  of  the  droplet -size  data 
obtained  under  the  present  program. 
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5 


EXPERIMENTAL  PROGRAM 


5.  1  DESCRIPTION  OK  APPARATUS 


5.1.1  Shock-Tube  Studies 

The  rapid  development  of  shock-tube  technology  in  recent  years  has  led 
to  wide  acceptance  of  the  shock  tube  as  a  means  of  inexpensively  pro¬ 
ducing  accurately  controlled,  compressible  flow  conditions  of  short 
duration  Properly  designed,  :he  shuck  tube  can  be  ustJ  .  j  generate 
flows  that  may  be  varied  over  rather  wide  ranges  without  necessitating 
the  use  of  complex  hardware  or  sophisticated  instrumentation 

Data  from  experimental  studies  made  on  the  performance  of  shock  tubes 
of  various  sizes  and  operating  ranges  have  verified  that  the  flow  conditions 
correspond  very  closely  with  values  predicted  by  theory.  For  all  but  the 
most  precise  investigations,  it  is  necessary  to  measure  only  <lj  the 
compression  and  expansion -chamber  ambient  pressures  and  temperatures., 
and  (2)  the  shock  velocity  to  determine  with  sufficient  accuracy  the  flow 
conditions  in  the  test  zone.  The  duration  of  flow  can  be  determined  from 
the  flow  parameters  and  the  lengths  of  the  various  portions  of  the  shock 
tube. 

A  brief  discussion  of  shock-tube  theory  is  presented  in  Appendix  A  of 
this  report.  Because  literature  on  shock-tube  theory  is  easily  available,, 
the  discussion  in  Appendix  A  has  been  limited  to  a  qualitative  description 
of  the  flow  in  the  shock  tube,  a  presentation  of  the  equations  describing 
the  flow  parameters,  and  the  equations  used  in  determining  the  duration 
of  flow 

The  shock  tube  used  in  this  study  (Figure  1)  employs  a  29-m  -long 
compression  chamber  made  of  steel  pipe  (4.  5-in  OD  x  4 -in  ID).  The 
end  of  the  compression  chamber  that  joins  the  diaphragm  holder  con¬ 
tains  a  transition  section  to  provide  a  smooth  transition  in  flow'  from 
the  circular  cross -section  of  the  pipe  to  the  2-in  -sq  cross-section  oi 
the  diaphragm  holder  and  expansion  chamber 

The  55-in  -long  expansion  chamber,  comprises  the  four  following  sections 

a. .  The  diaphragm  holder  (a  2- in  -thick  piece  of  rmcarta) 

b.  A  24-in  -long  shock  formation  section  with  0  3  -in. -thick  steel 
v.  alls 
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c  A  17-m  -long  test  section  containing  viewing  windows 

d  A  12-m.-'ong  section  made  of  0  5-m  -thick  plexiglass 

The  expansion  chamber  is  open-ended,  permitting  the  flow  behind  the 
shock  wave  to  exhaust  to  the  atmosphere 


5.  1  2  Diaphragm  Holder 

Efficient  operation  of  the  shock,  tube  requires  verv  rap'd  fern  oval  of 
the  diaphragm  separating  the  compressed  gas  m  the  compression  chamber 
from  the  gas  of  lower  pressure  in  the  expansion  chamber  Removal  of 
the  diaphragm  results  in  the  formation  of  a  plane  shock  front,  which 
moves  into  the  expansion  chamber  Behind  the  shock  is  a  region  ol  con¬ 
stant  flow  velocity,,  pressure,  temperature,  and  density  (the  region 
normally  used  as  the  test  zone) 

Bursting  or  releasing  the  diaphragm  can  he  accomplished  by  several 
methods  (1)  static  overpressure,  causing  mechanical  failure,  (2;: 
puncturing  the  diaphragm  material  with  a  pointed  object,  initiating  mech¬ 
anical  failure,  (3)  shock  overpressure,  which  exceeds  the  strength  of  the 
material,  and  (4)  electrical  discharge  or  exploding  wires  The  1.  tter 
method  was  chosen  for  use  in  this  study  Three -mi! -diameter  Nichrome 
wire,  sandwiched  between  layers  of  plastic,,  was  exploded  to  melt  and  cut 
the  plastic  to  initiate  opening  of  the  diaphragm 

Mylar  sheet  plastic,  of  0  5  and  1  mil  thicknesses,  was  selected  as  the 
diaphragm  materia’.,  since  it  is  easily  ruptured  and  does  not  shatter 
into  fragment that  may  be  carried  into  the  test  zone  Pressure  tests 
showed  that  a  0.  5-mil  thickness  of  Mylar  stretched  over  a  2-  x  2-in 
opening  is  capable  of  withstanding  approximately  a  13-psi  pressure 
differential  before  rupturing,  and  that  the  maximum  differential  is  pro. 
port,onal  to  diaphragm  thickness  regardless  of  whether  the  total  thickr.v -s 
is  in  one  layer  or  several 

Figure  2  shows  the  various  components  of  the  diaphragm  holder  and  the 
wire-positioning  jig  which  was  used  to  determine  the  proper  length  ol 
exploding  wire  and  to  ens<.ie  its  proper  position  v.  hi.  <  b<  .ng  fastened  to  i  he 
evelets  on  the  high-voltage  terminals  The  eyelets  can  be  seen  piotiud- 
.ng  into  the  2-in  -sq  opening  in  the  micarta  base 

F'lgure  3  shows  the  wire -positioning  jig  being  used  to  prepare  a  spr<  iP’cn 
diaphragm  and  exploding  wire  A  single  sheet  of  plastic,,  placed  over  the 
opening  in  the  micarta  base,  is  held  it.  p‘.,  ce  with  small  dabs  of  silicone 
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Figure  2.  Diaphragm  Holder  Component* 


grease  (or  similar  agent);  the  wire-positioning  jig  is  then  placed  over 
the  plastic  and  oriented  such  that  the  enlarged  holes  at  the  ends  of 
Slot  AD  (shown  in  Figure  3)  are  directly  over  the  eyelc.  .r.  the  high- 
voltage  terminals. 

One  end  of  the  exploding  wire  is  threaded  through  Hole  A  and  a  small 
perforation  in  the  sheet  of  plast’c,  and  is  attached  to  the  eyelet  on  the 
high-voltage  terminal..  The  wire  is  then  pulled  taut  and  is  passed 
around  the  2-in.-sq  opening  so  that  the  positioning  pins  at  A,  B,  and 
C  catch  the  wire  and  hold  it  in  position  The  loose  end  of  the  Nichrome 
wire  is  then  passed  down  throueh  Hole  D  and  the  plastic  sheet,  threaded 
through  the  terminal  eyelet,  and  tied  or  wrapped  snugly,  i’ne  wiring 
steps  can  be  expedited  by  running  the  Nichrome  wire  through  a  6-in 
length  of  capillary  tubing  (0  0625- in  OD  x  0  036-in  ID),  which  can  then 
be  used  as  a  needle  to  puncture  the  plastic  and  thread  the  wire  through 
the  eyelets. 

After  the  wire  is  in  position  and  fastened  to  the  high-voltage  terminals, 
the  wire-positioning  jig  is  removed  by  slipping  the  Nichrome  wire  off 
the  positioning  pins  and  lifting  the  jig  away  from  the  micarta  base  The 
loop  of  exploding  wire  slips  through  Slot  AD  and  is  left  attached  in 
position.;  A  second  sheet  of  plastic  is  then  placed  over  the  first,,  with 
th  loop  of  exploding  wire  sandwiched  between  the  sheets  of  Mylar.  The 
diaphragm-clamping  plate  is  placed  over  the  plastic  and  is  held  in  place 
w’ith  two  screws.  Figure  4  shows  the  diaphragm  holder,  complete  with 
a  specimen  diaphragm  and  simulated  exploding  wire 

The  Nichrome  wire  is  exploded  by  connecting  the  voltage  terminals 
across  a  capacitor,  which  is  charged  to  a  high  voltage  A  short  experi¬ 
mental  study  was  made  to  determine  the  optimum  energy  to  use  in 
exploding  the  wire;  too  much  energy  initiated  shock  waves  from  the 
explosion  of  the  wire,  wh.le  low  energy  caused  the  wire  to  burn  slowly, 
introducing  time-delay  problems.  It  was  determined  that  a  9-pf  capacitor 
charged  to  1.  2  kv  provided  sufficient  energy  (6  5  j)  to  explode  the 
Nichrome  wire  (35  to  40  ohms  dc  resistance)  without  time  delay  or  shock- 
wave  effects 

Appearance  of  the  diaphragm  after  a  test  run  is  shown  in  Figure  5;,  as 
illustrated,  the  explosion  of  the  wire  cuts  a  door,  hinged  along  one  side, 
that  opens  to  release  the  air  from  the  compression  chamber 

Behavior  of  the  diaphragm  during  the  opening  sequenci  was  investigated 
by  aking  high  speed  motion  pictures  of  the  process  These  tests,  the 
results  of  which  are  shown  in  Figures  6  and  7,  were  made  with  the 
expansion  section  removed  and  the  cam--  :-taced  to  view  the  diaphragm 
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Figure  4.  Diaphragm  Holder.. 


Figure  5.  Diaphragm  After  Test.' 


phragm  Opening  at  Low  Pressure 


Diaphragm  Opening  at  High  Pressure. 


from  a  position  on  the  shock-tube  axis  looking  toward  the  eompress’on 
chamber.  The  2-  x  2-in.  diaphragm  is  shown  as  it  is  cut  along  three 
sides  by  the  exploding  wire  and  begins  to  open  out  of  the  plane  of  the 
picture.  The  normally  transparent  Mylar  plastic  of  the  diaphragms 
was  made  opaque  for  these  tests  by  spraying  the  undersurface  with  fiat 
white  and  silver  paint 

Figure  6  shows  the  opening  sequence  for  a  pressure  of  3.  1  psig  in  the 
compression  chamber,  and  Figure  7  contains  the  sequences  for  pressures 
of  7.  5,  12,  and  17.  5  psig.  The  first  frame  in  each  sequence  represents 
the  instant  at  which  the  electrical  energy  was  applied  to  the  terminals  to 
explode  the  wire.  Electrical  arcing  ct  the  junction  of  the  exploding  wire 
with  the  terminals  was  visible  after  approximately  0  25  ms  or  less  in  all 
of  the  tests.  With  the  exception  of  the  test  conducted  at  the  lowest  pres¬ 
sure,  initial  rupture  and  movement  of  the  diaphragm  occurred  within  0  3 
ms,  and  complete  opening  occurred  in  less  than  1  ms. 

At  3.  1  psig,  the  sequence  was  considerably  delayed,  even  though  notice¬ 
able  arcing  and  burning  occurred  at  the  same  time  as  evidenced  in  tests 
performed  at  higher  pressures  Significant  rupture  and  movement  of 
the  diaphragm  did  not  occur  until  after  approximately  2  ms  elapsed,,  and 
the  opening  was  not  completed  until  approximately  3.4  ms  elapsed  The 
opening,  from  initial  rupture  to  completion,  required  approximately  1  5 
ms. 


Another  noticeable  feature  in  Figure  6  is  the  apparent  reclosing  of  the 
diaphragm,  an  effect  that  can  be  seen  ir.  the  last  row  of  frames  in  the 
sequence.  This  reclosing  is  apparently  due  to  reverse  flow,  which  was 
caused  by  the  slight  vacuum  created  by  the  inertia  of  the  air  rushing  out 
of  *he  compression  chamber. 

It  must  be  remembered  that  the  absence  of  the  expansion  chamber  creates 
a  different  flow  history  than  would  be  expected  with  the  chamber  in  place;, 
however,  the  data  in  Figures  6  and  7  should  be  valid  for  the  initial  open¬ 
ing  phases  of  the  diaphragm. 


5.-  1.  3  Instrumentation 

The  instrumentation  used  in  this  investigation  does  not  necessari'y  repre  - 
sent  the  combination  of  components  that  would  normally  be  selected  for 
use  in  a  shock-tube  study  of  th’s  nature  There  are  'arious  other  wei'- 
tried  methods  of  acquiring  shock-tube  operating  data  that,  oecavse  of 
economic,  space,  or  time  considerations,  might  be  more  desirable 
These  methods  are  described  in  detail  •  any  of  the  texts  and  reports  on 
shock  tubes  (References  6fc  through  70), 
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Items  of  equipment  used  in  this  study  were  all  readily  available  and 
fulfilled  the  various  requirements  for  operating  voltages,  response 
times,  and  accuracy.  It  was,  therefore,  more  expedient  to  adapt  the 
experimental  techniques  to  fit  the  available  instrumentation  rather 
than  to  acquire  instrumentation  to  fit  the  experimental  need. 

The  following  equipment  was  employed  during  the  studies' 

a. .  One  Beckman  Whitley  Model  3E6-3  Dynafax  Camera 

b. .  One  Beckman  Whitley  Model  357  Electronic  Flash  Unit 

c.  Two  Endevco  Model  2501-500  Piezoelectric  Pressure  Trans¬ 
ducers 

d.  Two  Photo  1P39  Tubes 

e.  One  Tektronix  Model  545  Oscilloscope  (with  Type  53/54C  Dual- 
Trace  Preamp  and  Polaroid  Camera) 

f  Two  Beckman  Model  7370  Counters 

g. .  Two  Shasta  Model  854R  Wide -Band  Amplifiers 

h.  One  Precision  Laboratory  Test  Gauge  (0  to  30  psig;  1/4%  full- 
scale  accuracy) 

i.  One  High-Voltage  Firing  Panel  (Model  177) 

The  Beckman  Whitley  Hvnafax  Camera  is  a  continuous -writing  fram¬ 
ing  camera,  capable  of  framing  rates  of  from  200  to  26,  000  frames/sec 
with  shutter  speeds  of  1.0.  2.  5,  or  5.  0  psec  at  maximum  framing  rate. 
This  camera  has  a  picture  capacity  of  224  frames  of  standard  16mm 
frame  size  on  a  33-7/8-in  length  of  35mm  film.  The  total  writing  time 
available  varies  from  1.  12  sec  at  200  frames/sec  to  8  62  ms  at 
26,  000  frames/ sec. 

A  Wollensak  Raptar  f/2.  8,  3-in.  Telephoto  lens  (Figure  1)  was  used  for 
all  of  the  tests,  although  various  lens -extension  lengths  were  required 
at  times  to  obtain  the  desired  magnification  and  field  of  view. 

Camer~  speed  is  controlled  by  a  variac  in  the  camera  base,,  and  the 
framing  rate  is  indicated  directly  on  a  meter  on  the  front  of  the  camera. 
More  accurate  determination  of  the  framing  rale  can  be  made  by  moni¬ 
toring  the  frequen  :y  of  an  ac-voltage  outpu.  f.  ,  the  camera.  The 
frequ?ncy  of  the  ac  --oltage  is  equal  to  tin-  speen  of  the  rotating  mirror  in 
the  camera,  and  is  directly  related  to  the  framing  rate. 
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The  Beckman  Whitley,  Model  357,  Electronic  Flash  Unit  (Figure  8)  is 
a  light  source  that  is  designed  for  use  with  the  Dynafax  camera.  This 
unit  is  capable  of  supplying  single,  square  pulses  of  high-intensity, 
cold  light  of  adjustable  duration.  Pulse  durations  available  are 
8.60,  11.  15,  14.85,  and  22.  35  ins,  with  rise  and  decay  times  of 
approximately  25  psec  (less  than  the  frame  separation  at  the  highest 
framing  rate). 

Five  methods  of  triggering  the  light  source  arc  available:  (1)  making 
contact  between  the  triggering  wires,  (2)  breaking  contact  between  the 
triggering  wires,  (3)  using  a  positive -voltage  pulse,  (4)  increasing  the 
illumination  that  falls  on  a  photocell  (1P39)  supplied  witu  tne  unit,  and 
(5)  decreasing  the  illumination  that  strikes  the  photocell.-  The  latter 
method  of  triggering  was  used  in  this  study;  the  falling  drops  of  liquid 
were  permitted  to  interrupt  a  beam  of  light  that  was  focused  on  the 
photocell. 

Figure  9  shows  the  arrangement  of  the  photocell  and  light  beam.  Light 
from  an  automobile -headlamp  type  of  bulb  m  the  light  source  is  focused 
into  a  beam  by  lenses  mounted  in  the  light-source  tub.-,  and  is  directed 
across  the  center  of  the  test  section  through  holes  in  the  test-section 
wall  to  strike  the  photocell  tube.  Interruption  of  the  light  beam  gener¬ 
ates  a  voltage  change,  which  triggers  the  flash  unit.- 

Shock -velocity  measurements  were  made  during  each  test  by  monitoring 
the  shock  front  time-of -arrival  at  two  pressure  transducers,  mounted 
14.  25  in.  apart  in  the  test-section  wall.-  These  transducers  (Figure  10) 
behave  electrically  as  capacitors  whose  charge  is  directly  proportional 
to  the  pressure  applied  to  the  transducer  face  (as  a  result,  they  respond 
only  to  pressure  changes).  With  proper  impedance  matching  and  control 
ol  system  KC  time  constants,  the  transducers  can  be  calibrated  to 
accurately  measure  both  the  magnitude  and  time  history  of  a  rapidly 
fluctuating  pressure..  For  pu>  poses  of  this  study,  it  was  only  necessary 
to  use  the  transducers  to  sense  t''e  step  change  m  pressure  that 
occurred  as  the  shock  front  passed  over  the  transducer  face. 

Because  the  transducer  signals  were  '.‘.a'  small  (4  to  200  mvolts)  to 
consistently  trigger  the  oscilloscope  and  counter,  it  was  necessary  to 
amplify  the  signals  by  passing  them  through  Shasta  Wide-Band  amplifiers 
(Figure  8)  set  at  a  gain  of  40  db.  Output  signals  from  the  Shasta  ampli¬ 
fiers  were  fed  into  both  the  Beckman  Model  7370  Counter  and  the 
Textronix  oscilloscope  (also  shown  in  Figure  8).  Th_  signal  from  the 
1  wer  transducer  was  used  to  start  the  counter,  which  had  been  set  on  the 
'  Time  B-A,'  mode  and  to  trigger  the  sweep  of  the  oscilloscope ;  the  i  -glial 
from  Jne  upper  transducer,  used  to  sto.  .*  counter,  was  fed  into  the 
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Figure  8.-  Instruni2ntatio<  . 
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Figure  10.  Side  View  of  Shock-Tube  Test  Section 


vertical-deflection  circuit  of  the  oscilloscope  When  operated  in  this 
manner,  the  counter  indicated  to  the  nearest  0  1  psec  the  time  elapsed 
between  the  start  and  stop  signals,  and  thus  represen*  the  time 
required  fcr  the  shock  front  to  travel  the  distance  between  the  trans¬ 
ducers 

Times  between  start  and  stop  signals  were  also  obtained  from  the 
oscilloscope  pictures  by  determining  the  length  of  the  sweep  trace  before 
it  was  vertically  deflected  by  the  stop  signal  Accuracy  comparable  to 
that  of  the  counter  was  obtained  from  the  oscilloscope  by  delaying  the 
start  of  the  sweep  for  an  accurately  known  interval  and  observing  only 
the  final  portion  of  the  trace  at  a  much  taster  sweep  rate  oscilloscope 
data  were  primarily  used  for  comparisons  with  the  counter  data  and  to 
estimate  possible  errors  in  time  measurement  as  affected  by  differences 
in  wave  shape  between  the  start  and  stop  signals  and  differences  in  trigger 
levels  between  start  and  stop  circuits  of  the  counter 

The  counter  shown  at  the  top  of  Figure  8  was  used  to  determine  the 
camera  speed  An  ac  signal  from  the  camera  (with  a  frequency  equal 
to  the  rotating-mirror  speed)  was  fed  into  this  counter,  which  was  set 
on  the  Events  Per  Unit  Time  (EPUT)  mode  and  thus  counted  the  number 
of  ?c  cycles  occurring  per  second.  Framing  rates  and  frame-separation 
times  were  then  calculated  lrom  the  mirror  speed 

Prior  to  the  test  runs,  the  pressure  gage  illustrated  in  Figure  8  was 
observed  to  determine  the  pressure  level  in  the  compression  chamber 
(It  was  necessary  to  measure  this  pressure  to  predict  the  magnitude  of 
shock  velocity  and  air  velocity  for  each  test  )  Shown  below  the  pressure 
gage  is  the  high-voltage  firing  panel,  which  was  used  to  explode  the  wires 
to  burst  the  diaphragms.  This  panel  ccntains  triggering  circuits,  a 
high-voltage  power  supply  (to  charge  the  9-pt  capacitance'),  and  firing 
circuits  (to  apply  the  energy  stored  in  the  capacitors  to  the  exploding 
wire)..  The  firing  and  triggering  cir<  its  are  designed  to  have  response 
times  of  less  than  1  psec;  triggering  is  accomplished  by  the  voltage 
pulse  from  a  1P39  photocell,  receiving  light  from  the  flash  unit  used 
with  the  Dynafax  camera 

At  the  bottom  of  Figu-e  8  is  the  control  panel  for  the  photocell  light 
source;  this  panel  contains  a  battery  charger,  which  was  used  to  re¬ 
charge  the  automobile  battery  that  .supplied  *he  dc  voltage  to  the  lamp 
in  the  light  source  The  lamp,  battery,  and  charger  were  connected  to 
a  three-way  switch  with  "off,  "  "operate,.  "  and  "recharge''  positions 
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5.1.4  Mic  roburet 


Formation  and  control  of  single  drops  of  liquid  are  most  easi'y  accorn 
piished  by  using  a  hypodermic  syringe  and  releasing  the  drops  .i<= 
desired  from  the  tip  of  a  section  of  capillary  tubing.  The  size  ol  the 
drops  is  determined  by  the  liquid  density,  surface  tension,  and  the  out 
side  diameter  of  tne  capillary  tubing  on  which  they  are  formed. 

If  it  is  assumed  that  the  drop  separates  fr.*m  the  tubing  at  the  instant 
when  the  weight  of  the  liquid  in  the  drop  is  sufficient  to  overcome  the 
surface-tension  force  attaching  the  drop*  to  the  tubing,  an  expression 
of  the  following  form  can  be  written  (Reterem-c  i  1) 

Drop  weight  =  kD^,cr 

where 

k  =  constant  of  proportionality  (experimentally  determined  = 

3.  8) 

D.j,  =  outside  diameter  of  capillary  tubing 
a-  =  liquid  surface  tension 

From  Equation  59,  the  expected  drop  size  can  be  written  as 


where 

D  =  drop  size 
P  =  liquid  density 

Figure  11  shows  drop  sizes  as  a  function  of  t  *bing  diameter  *or  the 
different  liquids  used  during  the  study.  Examination  of  the  curves 
indicate  that  it  is  not  practical  to  attempt  to  produce  drops  smaller 
than  about  1.  0  mm  in  this  manner  because  tubing  of  very  si* .all 
diameter  would  be  required. 
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Figure  11.  Drop  Sizes  versus  T*  !>e  Diann-ter  for  Various  liquids 
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io  produce  drops  significantly  smaller  tnan  1.0  mm,  a  microburet 
similar  to  that  used  by  Asset  (Reference  72)  was  fabricated.  A  short 
length  (approximately  0.  5  in.  )  of  0.  008-in.  -OD  x  0.  00  t~m.  -ID  tubing 
was  inserted  and  sealed  in  the  tip  of  an  8-in.,  -long  (0.  036-i.i.  -OD) 
hypodermic  needle.  A  length  of  0.  0625-in.  -OD  x  0.  051-in.  -ID 
capil’ary  tubing  was  slipped  over  the  hypodermic  needle  and  attached 
to  a  T  fitting  near  the  center  of  the  hypodermic  needle;  the  0.  008  in.  - 
OD  tubing  protruded  approximately  0.  125  in.  beyond  the  up  of  the 
0.  0625-in.,  capillary  tubing.  Figure  9  shows  the  0.  008-in.  -OD  tubing 
piotruding  from  the  capillary  of  larger  diameter  (the  T  fitting  is  some¬ 
what  obscured  by  the  ring-stand  clamp  holding  the  microburet  in  place). 

Air,  under  pressure,  is  admitted  to  the  T  fitting  and  flows  through  the 
annular  space  between  the  two  pieces  of  tubing  and  out  around  the  tip 
of  the  0.  008-in.  tubing.  As  the  droplets  begin  to  form  at  the  tip  of  the 
small  capillary,  they  are  blown  away  by  the  force  of  the  airflow.  By 
adjusting  the  rate  of  air  flow,  the  drop  sizes  can  be  varied  at  will  over 
a  considerable  range;  however,  once  a  fixed  airflow  is  established,  the 
drop  sizes  are  as  consistent  as  those  produced  with  no  airflow. 

Droplets  as  small  as  500  p  were  produced  for  the  study  with  no  difficulty. 
Droplets  smaller  than  this  size,  however,  showed  a  tendency  to  follow 
erratic  trajectories  after  being  blown  off  the  capillary  tubing  (because 
of  both  their  small  mass  and  the  air  turbulence  created  by  the  higher 
air  flow  rate).  Two  problems  resulted  from  the  preceding  behavior: 

(1)  it  was  difficult  to  get  the  drops  to  consistently  enter  the  light  beam 
and  trigger  the  flash  unit,  and  (2)  the  drops  were  frequently  outside  the 
field  of  view  of  the  camera.  As  a  result  of  these  problems,  the  test 
data  in  this  report  are  limited  to  drop  sizes  of  approximately  500  p 
and  above.. 

Figure  12  shows  the  arrangement  of  the  microburet,  the  tygon  tubing 
for  the  air  supply,  the  pressure  gage,  and  the  valve  used  to  adjust  the 
airflow  through  the  microburet.  Drops  larger  than  1. 0  mm  required 
for  the  study  were  all  created  by  selecting  various  sizes  of  capillary 
tubing  to  attach  to  the  end  of  the  hynodermic  syringe. 


5.  2  TEST  PROCEDURE 

Prior  to  each  run,  each  item  of  equipment  and  instrumentation  was 
checked  to  ensure  ils  proper  operation,  thereby  preventing  invalidation 
of  .he  test  by  failure  to  obtain  all  of  the  desired  data  After  the  preced¬ 
ing  steps  were  taken,  the  sequence  of  events  described  in  the  following 
section  occurred  during  the  test  runs.. 
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5.2.1  Photographic  Tests 


At  the  start  of  the  photographic  tests,  the  pressure  in  the  compression 
chamber  was  raised  to  the  desired  value.  Verv  slow  leakage  usually 
occurred  around  the  diaphragm  holder,  and  some  adjustment  of  the  air 
supply  was  necessary  to  stabilize  the  chamber  pressure.  As  the  adjust¬ 
ment  was  being  made,  the  capacitors  in  the  high-voltage  panel  were 
charged  to  1. 2  kv,  and  the  panel  was  adjusted  to  maintain  the  charge. 

A  second  operator  gradually  increased  the  speed  of  the  camera  until  the 
desired  number  of  revolutions  per  second  was  displayed  on  the  mirror- 
speed  counter.  The  camera  shutter  was  then  opened,  and  a  drop  of 
liquid  was  squeezed  from  the  hypodermic  needle  and  was  permitted  to 
fall  into  the  light  beam  to  initiate  the  test  sequence. 

Because  of  its  high  framing  rate  and  short  exposure  time,  the  Dynafax 
camera  requires  a  very  intense  source  of  light  for  adequate  exposure; 
therefore,  at  high  framing  rates  the  camera  shutter  may  be  left  open 
(exposing  the  film  to  normal  room  lighting)  for  as  long  as  30  sec  without 
film  fogging.  As  th<-  falling  drop  of  liquid  entered  the  light  beam,  ti.e 
electronic  flash  unit  was  triggered  to  emit  an  8.2-ms  pulse  of  light, 
which  provided  the  illumination  for  proper  film  exposure.  The  pulse  of 
light  from  the  flash  unit  also  triggered  the  high-voltage  panel,  exploding 
the  Nichrome  wire  and  releasing  the  air  in  the  compression  chamber. 

The  camera  continued  to  take  pictures  of  the  falling  drop  for  approxi¬ 
mately  2.  5  ms,  the  time  required  for  the  shockwave  to  form  and  travel 
the  distance  from  the  diaphragm  to  the  test  zone.  As  the  shockwave 
passed  through  the  test  section,  the  falling  drop  was  subjected  to  the 
onset  of  flow  behind  the  shockwave  in  the  time  required  for  the  shock 
to  traverse  the  drop  diameter  (less  than  10  (isec).  The  flow  simul¬ 
taneously  initiated  deformation  and  breakup  of  the  drop,  subjecting  it 
to  large  acceleration  forces  sufficient  to  halt  its  downward  fall  and 
carry  it  out  of  the  top  of  the  camera's  field  of  view  in  a  few  milliseconds. 

The  duration  of  constant -flow  conditions  behind  the  shockwave  was 
approximately  3  ms,  a  time  span  that,  except  for  isolated  cases,  was 
more  than  sufficient  to  cause  breakup  and  carry  the  shattered  droplet 
out  of  the  field  of  view.  Since  the  camera  had  a  total  writing  time  of 
8  ms,  and  only  2..  5  ms  of  this  time  was  used  before  arrival  of  the 
shockwave,  the  remaining  5.  5  ms  was  more  than  sufficient  to  record 
the  events  of  interest.  Passage  of  the  shockwave  over  the  pressure  trans 
ducers  at  the  bottom  and  top  ot  the  test  sect. on  generated  voltage  signals, 
which  triggered  the  counter  and  oscilloscope  to  record  the  time  required 
for  the  shockwave  to  traverse  the  listanrs  between  the  transducers 


The  test  was  terminated  by  essentially  following,  in  reverse,  the  step 
taken  to  prepare  lor  the  ICSE¬ 


S'.  ?.  Size-Distnbutic.i  Tests 

The  procedure  followed  in  conducting  the  tests  to  sample  the  droplet-size 
distribution  resulting  from  breakup  of  larger  single  drops  was  similar  to 
that  discussed  in  Paragraph  5.  2.  1.  It  varied,  in  that  the  Dynafax  camera 
was  not  utilized  and  it  was  necessary  to  place  the  shock  tube  in  a  horizontal 
position  to  gain  operating  space  beyond  the  open  end  of  the  expansion  cham¬ 
ber. 

To  remove  all  unnecessary  airflow  obstructions  from  the  test  section, 
the  droplets  were  introduced  into  the  test  ? one  by  mounting  the  hypo¬ 
dermic  syringe  and  needle  outside  the  shock  tube  and  permitting  the 
drops  to  fall  from  the  needle  lip  through  a  hole  in  the  test-section  wall. 

The  light  source  and  photocell  were  relocated  and  adapted  so  that  when 
the  falling  drop  approached  the  shock-tube  centerline,  the  light  beam 
was  interrupted  to  both  trigger  the  high-voltage  firing  pane!  and  explode 
the  Nichrome  wire  to  release  the  air  in  the  compression  chamber. 

Test  preparations  consisted  of  arranging  12  standard  microscope  slides 
in  a  plywood  holder,  creating  a  4-  x  6-in.  plane -sampling  surface.  The 
plywood  holder  was  then  mounted  on  a  moveable  stand,  with  the  center 
of  the  sampling  surface  concentric  with  the  shock-tube  axis  and  the 
plane  of  the  surface  normal  to  the  tube  axis;  the  stand  was  positioned  to 
place  the  sampling  surface  a  predetermined  distance  from  the  open  end 
of  the  expansion  chamber. 

The  test  was  conducted  following  the  same  preeeuure  used  in  the  photo¬ 
graphic  tests.-  After  the  test  was  completed,  the  microscope  slides 
were  removed,  and  photomicrographs  of  the  sampled  material  were 
taken  at  three  locations  along  the  centerline  of  the  slide.-  To  provide 
data  from  which  the  radial  distribution  of  sampled  material  could  be 
determined,  recordings  were  made  of  the  coordinates  of  each  of  the 
photomicrographs  on  the  slide  surface.  The  sampled  material  was  also 
washed  from  the  slides  and  put  into  a  solution,  permitting  assays  to  be 
made  ol  the  amount  of  material  collected  for  material-balance  calculation! 


5.3  PARAMETRIC  STUDIES 


5.  3  )  Measurement  of  Breakup  Time 

Correlation  of  the  experimental  results  »\'th  the  theoretical  model  of 
breakup  mechanism  developed  in  Section  1  of  this  report  requires  the 
following  data; 
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Time  required  for  breakup  to  occur 

b.  Physical  properties  of  the  liquid  being  shattered 

<..•  Airflow  conditions  to  which  the  liquid  is  subjected 

To  prove  useful,  th»  correlation  requires  that  the  data  represent  a  wide 
range  ol  values  ol  the  various  sigmlicant  parameters. 

Although  considerable  experimental  data  on  single -drop  breakup  are 
available  in  the  literature,  this  information  proved  to  be  of  limited  use 
in  the  present  study  for  two  rearons;  (?)  it  does  not  cove1"  tV,e  range  of 
variation  in  liquid  physical  properties  needed  for  correlation  with  theory, 
and  (2)  n  ch  of  the  work  reported  did  not  include  the  use  of  high-speed 
motion  pi  -tare  equipment  to  record  the  breakup.  As  a  result,  data  on 
the  break  o  time  are  not  readily  available.; 

A  comple'  e  parametric  study  was  therefore  undertaken  to  record,  with 
high-speeo  motion  picture  equipment,  the  breakup  of  liquid  drops  of 
different  sizes  and  different  physical  properties  while  subjected  to  a 
wide  range  of  flow  conditions. 

The  following  listing  shows  the  parameters  employed  for  the  test  and 
the  range  of  values  over  which  they  were  varied. 


Parameter 

Approximate  Range  of  Variation 

Relative  velocity 

50  to  450  ft/ sec 

Surface  tension 

18  to  487  dynes/cm 

Liquid  density 

0.  75  to  14  gm/cra^ 

Liquid  viscosity 

0.  5  to  170  centipoise 

Drop  size 

500  to  3000  ^ 

To  cover  the  preceding  range  of  variation,  several  liquids  were  chosen 
to  be  used  in  the  study..  These  liquids  and  their  physical  properties  are 
shown  in  Table  1. 

Experimental  results  obtained  during  the  study  are  presented  in  Table  2 
plates  made  from  motion-picture  films  showing  the  qualitative  behavior 
of  the  drops  during  breakup  are  provided  i  Appendix  B.  The  breakup- 
time  values  listed  in  Table  2  represent  the  -lapsed  time  from  the  first 
signs  of  deformation  until  the  drop  begins  co  break  up. 
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Table  1.  Physical  Properties  of  Liquids. 


Liquid 

Temperature 

(°C> 

Surface 
Tension 
(dynes/ cm) 

Viscosity 

(centipoise) 

Density 
(gm,/  cra^) 

Water 

25 

72.  00 

0  890 

0.  998 

Mercury 

20 

487  00 

1. 550 

13. 600 

Silicone  fluids 

GESF  (9t )  0.65 

25 

17.  50 

0.  470 

0.  758 

GE  SF  (96)  10 

25 

22.  30 

9.  520 

0.  938 

GE  SF  (96)  200 

25 

2  3.  :o 

!'”  ' 

0.  966 

Bis 

25 

27.  60 

5.  970 

0.915 
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Table  2.  Experimental  Breakup  Time  Results 


Test 

Drop 

Size 

(mm) 

At  r 

Velocity 

(ft/sec) 

Relative 
Velocity 
(ft/ sec) 

B-:-"  P 

Time 

(msec) 

WATER  | 

1  17 

2.  9 

128 

134 

1  2Q  to  1.  45 

119 

'.0 

214 

221 

0.  t>9  to  0  73 

121 

3  0 

357 

362 

0.  24  to  C.  27 

122 

3.0 

448 

455 

0  14  to  0  18 

123 

3  0 

63 

71 

3.  24  to  3.  32 

158 

2.0 

63 

69 

l  63  to  2.  83 

186 

1.4 

201 

107 

1.  37  to  1  49 

187 

1.9 

78 

84 

1.  64  to  1  30 

188 

1.9 

47 

53 

2.  35  to  2.  58 

189 

1  8 

189 

195 

0.  43  to  0.  51 

191 

1.9 

347 

350 

0  20  to  0.  3  1 

192 

■  7 

43t> 

441 

0  12  to  0  23 

214 

1  0 

57 

63 

2.  08  to  2.  24 

217 

1.2 

76 

81 

1.  45  to  1.61 

236 

0  76 

79 

85 

1.  16  to  1  27 

237 

0.  78 

68 

75 

1.07  to  1  15 

2  38 

0.78 

90 

95 

0.  89  to  0  97 

239 

0.  78 

103 

ioe 

0  77  to  0.  84 

240 

0  81 

1  17 

120 

0  69  to  0  77 

241 

0.6'* 

55 

61 

1  89  to  2  00 

261T 

0.53 

122 

127 

0.  66  to  0  77 

261B 

0.  57 

*22 

127 

0  62  to  0.  66 

262 

C  53 

197 

203 

0  5C  to  0  54 

263 

0  62 

333 

337 

0.  12  to  o.  n 

2t>4 

0  44 

425 

432 

0  04  *  >  0.  12 

267 

1  9 

79 

1.  15  to  1  27 
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Table  2.  Continued 


*V»t 

Drop 

Size 

(mm) 

At  r 

Vel  city 
(ft/ sec) 

Relative 

Velocity 

(ft/scc) 

T  line 
(m?cc) 

j  MERCURY 

*32 

2  1 

437 

4  oo 

0  ill  to  0  7  5 

135 

0  82 

351 

358 

0  "7  to  0,  85 

138 

0  69 

224 

230 

l  62  to  1  66 

|t>4 

1  2 

168 

173 

2  95  to  3  14 

|  GE  Sr<Vb)0  6** 

151 

2  3 

68 

75 

1  11  to  1  23 

152 

I.  o 

131 

137 

0  61  to  0.  b9 

153 

2  1 

120 

126 

0  43  to  0  o3 

154 

l.  2 

22S 

235 

0  16  to  0  27 

155 

1  9 

375 

381 

C.  10  to  0  14 

156 

7  ] 

467 

472 

0  08  to  0  12 

160 

1  4 

43 

49 

1  49  to  1.61 

193 

1  3 

436 

441 

0  08  to  0  12 

144 

1  5 

345 

350 

0  12  to  0  20 

195 

1  4 

195 

200 

0  20  to  0  27 

197 

1.4 

41 

48 

1  45  to  1  52 

[  CESF(96)!0 

mm 

B| 

450 

0  06  to  0  14 

mm 

228 

0  29  to  0.  41 

1 9 

n 

■ 

C»  **  u  1  00 

148 

2  0 

360 

>66 

0  10  o  0  22 

150 

2  0 

83 

90 

1  43  to  1  51 

1  4 

45 

51 

1  13  to  1  65 

iia 

1  1 

48 

1«3 

0  78  t*.  0  8»> 

199 

:  4 

188 

193 

0  35  * o  0  4  3 

200 

1  4 

338 

54  > 

t  12  too.  20 

201 

1  4 

430 

4  36 

0.  08  to  0  16 
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Table  2.  Continued. 


Test 

Drop 

Size 

(mm) 

Air 

Velocity 
(ft /sec) 

Relative 

Velocity 

(ft/sec) 

Breakup 

Time 

(msec) 

jj  CE  SF(96)200 

125 

2  1 

(7 

75 

4  00  to  4  35 

12? 

2.2 

11? 

123 

>28 

2  2 

121 

127 

1  85  to  2  01 

1 2^ 

2.2 

208 

215 

0  V2  to  1  00 

no 

2  1 

3C  / 

363 

0  51  to  0  59 

131 

2  > 

450 

456 

0.  3  3  to  0.  37 

>62 

1.4 

32 

38 

203 

i.  5 

75 

81 

204 

>  5 

292 

197 

0  .4  to  0  82 

205 

I  5 

340 

345 

0  39  to  0  43 

206 

1  5 

398 

403 

C  32  to  0  35 

208 

1  4 

63 

68 

232 

1  2 

161 

164 

243 

Si 

60 

63 

244 

o 

74 

77 

248 

Warn 

100 

106 

t  6?  trt  1  77 

249 

102 

107 

2  05  to  2  13 

251 

1  0 

114 

118 

1  54  to  1  65 

w 

1.0 

160 

165 

1  34  to  1  46 

25? 

0.9 

1?) 

1?9 

0  96  to  1  04 

258 

0  9 

181 

184 

0.  8<*  to  1  00 

68 


5.3.2  Size -Distribution  Tests 


Because  of  time  limitations  and  the  large  amount  of  work  involved  in 
determining  size-distributions,  these  tests  were  not  as  compre¬ 
hensive  as  the  breakup-time  tests.  This  investigation  was  limited  to 
the  use  of  one  liquid,  three  drop  sizes,  and  five  relative  velocities. 
Sampling  distance  was  also  varied,  but  only  for  purposes  of  determin¬ 
ing  the  optimum  distance  at  which  valid  samples  could  be  taken.- 

Bis  was  used  as  the  test  liquid  because  it  has  a  very  low  vapor  pressure, 
which  reduces  evaporation  oi  the  sampled  material  to  a  .:unimuni. 
Another  advantage  of  this  material  is  that  considerable  previous  exper- 
oce  had  been  gained  on  its  behavior  in  aerosol  form,  typical  spread 
factors  for  bis  droplets  on  glass  had  been  determined,  and  chemical- 
assay  techniques  to  determine  the  quantity  of  bis  sampled  had  been 
worked  out  and  tested. 

Drops  of  bis  (1.6  and  2.  7  mm  in  diameter)  were  subjected  to  relative 
velocities  of  approximately  100,  200,  300,  and  400  ft/sec,  and  the 
aerosol  from  the  shattered  drops  was  sampled  at  distances  of  from 
12  to  48  in.  from  the  open  end  of  the  shock -tube  expansion  chamber.. 

The  results  of  these  tests  yielded  size  distributions  for  the  shear  or 
stripping  type  of  breakup. 

Bis  drops  approximately  0.  6  mm  in  diameter  were  subjected  to  air 
velocities  of  approximately  70  ft/sec  to  determine  size  distributions 
resulting  from  bag  breakup,  which  requires  a  combination  of  small 
drop  size  and  low  relative  velocity.;  Samples  of  this  material  were 
obtained  only  after  the  1  -ft  plexiglass  section  had  been  removed  from 
ihe  expansion  chamber,  thereby  reducing  to  8  in.  the  distance  from 
the  point  at  which  breakup  started  to  the  end  of  the  expansion  chamber. 

In  these  tests,  the  sampling  surface  was  placed  approximately 
4-3/8  in.  from  the  open  end  of  the  expansion  chamber. 

Results  of  the  size-distribution  tests  are  presented  in  Appendix  C  m 
three  formsf  (1)  tabulations  of  assessment  data  taken  from  photo- 
micraphs  of  the  sampled  aerosols,  (2)  plates  made  from  photomicro¬ 
graphs  of  a  portion  of  the  sampling  tests,  and  (3)  plots  of  cumulative 
pert  ent  mass  vs  drop  size  for  all  of  the  sampling  tests. 
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'"ORRELATION  OF  RESULTS  WITH  THEORY 


6 . 


6  1  QUALITATIVE  DESCRIPTION  OF  BREAKUP 

The  following  sections  discuss  the  qualitative  behavior  of  liquid  drops 
undergoing  breakup  in  a  fast-moving  gas  stream.  The  results  were 
obtained  by  photographing  free-falling  drops  of  liquid  as  they  were 
struck  from  below  by  the  flow  behind  a  plane  shockwave  moving  upward 

As  noted  earlier  in  this  report,  substantial  evidence  by  various  authors,: 
e.g.  ,  Hanson  (Reference  30),  Priem  (Reference  32),  Sat'’  (Reference 
33),  and  Rabin  (Reference  37),  indicates  that  breakup  is  brought  about 
by  the  onset  of  the  flow  behind  the  shockwave  rather  than  from  the 
shockwave  itself.  The  time  required  for  the  shockwave  to  travel  over 
the  drop  (of  the  order  of  10  psec)  is  small  as  compared  with  the  elapsed 
time  between  shock  arrival  and  the  instant  the  drop  begins  to  shatter 
(approximately  500  psec  or  greater)..  Therefore,  it  is  highly  unlikely 
that  during  such  a  short  interval  of  time  the  shockwave  could  initiate  the 
events  observed  during  breakup. 

The  first  observable  stage  of  breakup  after  the  drop  is  subjected  to  the 
gas  flow  is  a  characteristic  flattening  of  the  downwind  side  of  the  drop. 
Such  behavior  was  also  noted  by  Engel  (Reference  36)  As  the  flattening 
of  the  leeward  side  continues,  a  noticeable  growth  in  the  lateral  diameter 
of  the  drop  takes  place  at  a  point  approximately  1  /  3  of  the  drop  diameter 
from  the  leading  face  of  the  drop,  or  at  the  point  approximately  coinciding 
with  the  aero  pressure  point  on  a  sphere  in  steady  llow. 

Flattening  of  both  the  windward  and  leeward  faces  continues  until  the 
drop  assumes  a  shape  similar  to  an  oblate  spherc.iu;  if  the  flow  velocity 
is  low  enough,  the  drop  may  be  flattened  into  a  disc  shape  with  a 
diameter  roughly  twice  that  of  the  original  drop.  At  this  point,  if  the 
relative  velocity  is  sufficiently  great,  the  flow  of  gas  over  the  front  face 
and  sides  of  the  drop  begins  to  initiate  surface  disturbances,-  which 
appear  to  be  small  wavelets  The  greater  the  relative  velocity,,  the 
greater  the  numt»~r  of  the  disturban~es.  As  these  disturbances  move 
away  from  the  stagnation  point  and  apptoach  the  edge  of  the  distorted  drop, 
they  are  stripped  off  and  carried  away  by  the  gas  flow  as  sheets  and  liga¬ 
ments  of  liquid,  which  subsequently  collapse  into  numerous  small  droplets. 

Erosion  of  the  periphery  of  the  deformed  drop  continues  until  all  of  the 
mass  of  the  drop  has  undergone  breakup  through  the  stages  of  sheet  and 
ligament  formation  and  collapse  into  droplets  that  are  carried  away  by 
the  gas  flow  Measurements  made  from  the  various  films  taken  during 
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the  parametric  study  showed  that  although  the  drop  undergoes  tre¬ 
mendous  acceleration,  the  breakup  process  is  so  rapid  that  the  bulk 
of  the  material  in  the  drop  acquires  no  more  than  -  pproximately  10  to 
20%  of  the  gr.s  velocity  before  it  is  reduced  to  aerosol  form,:  the  shower 
of  smaller  droplets  is  then  rapidly  accelerated  up  to  some  velocity 
near  the  gas  velocity. 

Figure  13  depicts  the  breakup  characteristics  described  in  the  preced-. 
ing  sections.  The  sequences  shown  are  for  3- mm  drops  of  water  sub¬ 
jected  to  various  relative  velocities  ranging  from  a  low  of  71  ft/ sec  to 
a  maximum  of  455  ft/sec.  The  relative  velocity  for  each  drop  is 
noted  above  the  first  frame  in  each  sequence,  and  the  figures  below 
each  frame  denote  the  elapsed  time  (in  milliseconds)  from  the  instant 
the  drop  was  subjected  to  the  flow. 

The  solid  diagonal  lines  between  the  sequences  are  "iso-time"  lines, 
denoting  times  of  roughly  1/2  and  1  ms  after  time  zero  By  following 
a  particular  diagonal,  it  is  possible  to  determine  the  effect  of  varying 
the  relative  velocity  on  the  rate  of  deformation  and  breakup  of  the  drop 
and  the  formation  of  surface  disturbances.  Also  of  interest  is  that  the 
liquid  being  stripped  from  the  periphery  of  the  drop  does  not  follow  a 
path  parallel  to  the  general  flow  (it  appears  to  diverge  to  form  a  cone¬ 
like  cloud  of  droplets  with  the  remaining  bulk  of  the  drop  located  at  the 
apex)  This  behavior  was  found  to  be  particularly  useful  in  establishing 
a  model  from  which  to  derive  the  expression  describing  the  expected 
droplet  sizes  resulting  from  breakup  (see  Section  4  4). 

The  effect  of  liquid  viscosity  on  breakup  is  shown  in  Figure  14,  where 
2.  1-mm  drops  of  silicone  oils  with  viscosities  of  0  47,  9.  52,;  and 
169-  5  centipoise  were  subjected  to  a  relative  velocity  of  127  ft/sec 
|<?n-time  lines  denoting  times  of  0  5,  1  0,  and  1 .  5  ms  after  time  zeio 
are  shown  for  reference. 

As  shown  in  Figure  14,  a  change  in  viscosity  of  more  than  two  orders  of 
magnitude  does  not  appreciably  alter  the  mode  of  breakup  Increasing 
the  viscosity  retards  the  deformation  process,-  however,  an  effect  that 
permits  the  drop  to  acquire  a  larger  velocity  (and  thus  reduce  the  relative 
velocity)  before  the  actual  stripping  process  sets  in  As  a  result,  an 
increase  in  the  size  of  the  droplets  stripped  off  would  be  expected  with 
increasing  viscosity. 

In  addition,  at  higher  viscosities,:  the  material  beirg  stripped  from  the 
periphery  of  the  drop  leaves  in  the  form  of  large  sheets  and  ligaments, 
rather  than  as  small  sheets  or  ligaments  or  even  d.  oplets  (as  occurs  at 
low  viscosities)  This  effect  is  shown  in  a  comparison  of  the  sequence 
for  the  0.  47  •  centipoise  viscosity  in  fc,g  /e  14  with  that  for  a  viscosity  of 
169  5  centipoise. 
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Ihe  large  sheets  and  ligaments  must  undergo  secondary  breakup  to 
produce  droplets,  and  they  appear  to  accelerate  more  rapidly  than 
the  bulk  of  the  drop;  therefore,  a  further  decay  m  rel  five  velocity 
is  likely  before  breakup  is  completed.  Again,  largei  droplets  should 
be  expected  from  the  liquids  of  higher  viscosity 

Unfortunately,  time  limitations  prevented  undertaking  the  work  neces-. 
sary  to  experimentally  investigate  the  effect  of  viscosity  change  on  the 
droplet  sizes  produced  during  breakup. 

Lane  (Reference  28),  in  conducting  one  of  the  first  comprehensive  studies 
of  breakup  in  a  steady  air  stream,  encountered  a  mode  cl  o.vakup  quite 
different  from  the  stripping  behavior.  He  noted  that  the  drop  flattened  to 
u.  disc  configuration,  and  instead  of  stripping  occurring  at  the  edges,,  the 
ctnter  of  the  disc  was  blown  out  in  the  direction  of  flow  to  form  a  large 
hohow  bag  rttached  to  a  roughly  circular  rim  (the  open  mouth  of  the  bag 
faced  the  oncoming  flow)..  Measurements  indicated  roughly  70%  of  the 
original  drop  mass  was  contained  in  the  rim..  Lane  also  studied  droplet 
breakup  using  fast  transient  blasts  of  air  and  obtained  the  stripping 
behavior..  He  thus  concluded  that  the  two  different  breakup  modes  were 
due  to  the  differences  between  steady  and  transient  flow 

Hanson  et  al.  (References  30  ar.d  31)  showed  that  both  modes  of  breakup 
could  be  produced  in  transient  blasts.  They  noted  that  bag  breakup 
occurred  when  drop  sizes  and  relative  velocities  were  close  to  the 
critical  values. 

As  a  part  of  this  study,  specific  efforts  were  undertaken  to  investigate 
both  modes  of  breakup  and  to  attempt  to  expand  the  knowledge"  of  the 
characteristics  of  the  two  types  of  behavior.  The  results  of  these  efforts 
arc  summarized  in  Figuie  13,  which  shows  the  effect  of  increasing 
relative  velocity  on  the  mode  of  breakup  Water  was  used  as  the  test 
liqt  id,  and  the  drop  sizes  were  held  nearly  constant  The  first  sequence 
for  a  relative  velocity  of  61  ft/sec.  clearly  demonstrates  the  bag  breakuo 
mode,  i.  e.:,  typical  flattening,  bag  formatj.-n.  bag  rupture  and  collapse 
of  the  rim  into  a  ring  of  large  drops.  Hanson  et  al  (References  30  and  31) 
discovered  that  after  the  formation  of  the  bag  has  taken  place,,  circular 
ripples  (similar  to  the  surface  waves  resulting  fi  >m  a  stone  thrown  into 
a  quiet  pond)  form  o.-i  the  thin  surface  of  the  bag  and  apparently  initiate 
bag  rupture. 

The  breakup  behavior  shown  in  the  second  sequence  of  Figuic  17  was 
p*‘  educed  bv  holding  the  drop  size  nearly  constant  and  increasing  the 
relative  velocity  to  85  ft/ sec.  In  this  case,  the  usual  flattening  and 
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neg  inning  of  bag  formation  occu:  red;,  however,  a  re-entrant  or  stamen 
portion  of  liquid  formed  at  the  center  of  the  flattened  drop  As  the  bag 
continued  to  form,  the  stamen  remained  attached  to  bag  and  elon¬ 
gated  to  form  a  ligament  of  liquid,  centered  within  tne  i-  ig  run  Tht 
usual  rupture  of  the  bag  a.,d  rim  occurred,  and  the  stamen  appeared  to 
collapse  into  a  series  of  larger  drops 

The  remaining  sequences  of  Figure  15  indicate  that  as  the  relative 
velocity  is  gradually  increased,  more  of  the  mass  of  the  drop  becomes 
concentrated  in  the  stamen  portion,  with  less  mass  available  to  form 
the  bag  portion.  This  effect  continues  with  increasing  velocity  until  a 
majority  of  the  drop  mass  is  containeu  in  the  center  stamen  portion  At 
this  point,  the  bag  no  longer  forms,;  instead  a  thin  rim  or  lip  of  liquid 
forms  around  the  periphery  of  tiie  center  potion  of  the  flattened  drop 
Rather  than  being  blown  out  into  a  bag,  the  lip  is  merely  stripped  off  t!u 
edge  of  the  drop  as  a  sheet,  which  subsequently  shatters  into  small 
droplets. 

A  good  example  of  the  behavior  that  might  be  expected  to  follow  after  the 
effect  shown  in  the  last  sequence  in  Figure  15  is  illustrated  in  the  first 
sequence  in  Figure  13. 

Based  on  the  results  summarized  in  Figure  15  and  other  test  evidence 
included  in  Appendix  B  of  this  report,  it  appears  that  there  is  a  smooth 
transition  in  behavior  in  going  from  the  bag  to  stripping  breakup  modes 
In  addition,  the  evidence  indicates  no  reason  why  a  different  transition 
would  occur  in  going  from  stripping  to  bag  breakup 

A  more  detailed  study  of  the  qualitative  aspects  of  breakup  can  be  made 
by  referring  to  series  of  plates  in  Appendix  B  These  plates  were  made 
it om  the  test  films  used  to  obtain  the  data  presented  in  Table  2.  All 
pertinent  information,  including  drop  sizes,  relative  velocity  type  of 
liquid,  and  time  between  frames,  is  presented  on  each  plate  for  reference 
purposes.  A  scale  factor  has  also  been  included  to  facilitate  recording 
length  measurements  directly  from  the  plates. 


6  2  AERODYNAMIC  BREAKUP  TIMES 

For  convenience  and  consistency,  the  exper’rnentally  derived  tjreuhi.ii 
values  in  Table  2  were  taken  as  the  interval  of  tune  between  the  firs; 
signs  of  deformation  and  the  time  at  which  the  drop  j  .st  began  lu  break 
V  is  evident,  however,  from  the  breakup  photographs  (Appendix  B)  th  .1 
the  entire  breakup  of  the  drop  requires  a  finite  tune  thereafter,  and  ♦  »>  .t 
this  time  is  usually  comparable  to  the  roqui'ed  to  initiate  the 

breakup.  Thus,  it  is  lo  be  expected  that  t  ere  should  be  e  small  factor 
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(of  th°  order  of  1.  5  to  3)  between  the  experimental  breakup  times  as 
measured  and  those  predicted  by  the  theory,  with  the  theoretical  times 
being  the  greater.  A  theory  for  the  aerodynamic  breakup  time  of  liquid 
drops  was  developed  in  Section  4.  3,  and  the  final  expression  without  re¬ 
gard  to  breakup  mechanism  (i.  e  ,  bag  or  stripping)  was  given  m  Equa¬ 
tion  48a.-  Figure  lb  compares  Equation  48a  to  all  of  the  experimental 
breakup  times;  the  agreement  is  seen  to  be  very  good  (except  perhaps 
for  a  small  constant  factor,  which  may  be  expected  on  the  basis  that  the 
experimental  times  represent  the  time  for  initiation  of  the  breakup  rather 
than  for  all  of  the  breakup  to  occur).: 

The  value  used  for  the  surface-tension  pre«~nre  in  the  theory  ?.&/d. 

A  factor  of  2  was  used  because  it  allowed  a  good  fit  of  the  data,  and  also 
because  if  a  larger  value  was  used  the  surface-tension  pressure  was  in 
some  cases  greater  than  the  aerodynamic  flow  pressure.  The  preceding 
would  imply  that  breakup  should  not  occur,  whereas  experimental  break¬ 
up  did  occur.  It  may  be  noted,  however,  that  if  a  drag  coefficient  greater 
than  unity  is  operative,  then  a  value  greater  than  2  could  be  used  (Equa¬ 
tion  49).  There  is  some  evidence  that  this  assumption  may  be  true,  as 
will  be  shown  later  when  the  values  of  the  drag  coefficient  throughout  the 
breakup  of  the  drop  are  computed.  However,  there  is  great  convenience 
in  not  attempting  to  use  a  drag  coefficient  in  the  breakup-time  equations.- 
Hence,  Equation  48a  may  be  considered  to  be  a  reliable  expression  in  pre¬ 
dicting  the  aerodynamic  breakup  times  of  liquid  drops  if  a  value  of  about 
’  is  used  for  k  in  Equation  49  and  the  drag  coefficient  is  considered  to 
be  unity. 

Euqation  48a  was  derived  using  the  maximum  liquid-flow  velocity,  but  it 
was  discussed  that  an  average  velocity  may  be  involved  if  a  hyperbolic 
velocity  distribution  was  established.:  Because  short  distances  were 
employed,  it  was  not  believed  that  this  involvement  co-id  _-ccor.  Figure 
17  shows  a  comparison  of  the  experimental  data  to  theory  if  an  average 
velocity  is  used.  The  agreement  with  theory  is  still  good,  except  that 
the  constant- factor  difference  is  slightly  larger  than  using  Equation  48a 
as  written.  From  the  photographs  of  the  times  involved  for  complete 
breakup  to  occur  once  breakup  is  initiated,  it  appears  that  use  of  the 
maximum  velocity  expression,  i.e.,  Equation  48a,  is  to  be  preferred. 

In  Section  4.  3,  an  explanation  was  presented  for  the  dual-mechanism 
breakup  behavior  oi  liquid  drops,  and  the  results  were  summarized  in 
Equations  51a  and  51b.  A  comparison  of  the  experimental  data  to  these 
equations  is  shown  in  Figures  18  and  19;  the  d«ag  coef-.cients  were  taken 
as  unitv  and  best  theoretical  fits  of  the  data  were  obtained  using  respec¬ 
tive  surface-tension  pressures  of  4<r/d  and  cr/d  for  the  stri>  ping  and  bag 
breakup  modes.  The  agreement  between  experiment  an-!  theory  for  strip¬ 
ping  breakup  is  very  good,  except  for  the  exp.ic  d  constant  factor. 
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The  agreement  for  bag  breakup  is  not  quite  so  good;  however,;  the 
agreement  can  oe  made  very  good  if  a  drag  coefficient  greater  than 
unity  is  used.  Since  the  factors  1/3  and  2/3  in  Equations  51a  and  b 
are  not  exact  but  merely  representative,  it  can  be  said  .hat  the  ex¬ 
perimental  data  support  the  hypothesis  presented  in  Section  4.  3  to 
explain  the  dual-mechanism  breakup  behavior  of  liquid  drops. 

A  theory  of  the  bag  breakup  times  of  liquid  drops  was  givei,  hy  Gordon 
and  was  summarized  in  Equations  24  and  26..  A  lomparison  of  Equa¬ 
tion  26  with  the  experimental  data  is  given  in  Figure  20;  the  agreement 
is  seen  to  be  fairly  good  (except  for  a  constant  factor),  although  several 
points  deviate  rather  considerably  from  the  data.  It  is  louud,  however,, 
that  the  Gordon  pressure  expression  as  written  in  Equation  23  predicts 
that  breakup  should  not  occur  for  a  number  of  the  drops  where  breakup 
was  experimentally  found  to  occur  (as  a  result  of  the  large  surface- 
tension  pressure,;  8<r/d,  which  he  used) 

The  preceding  data  were  replotted,  replacing  the  factor  8  by  4,  2,  and 
1.,  A  value  of  4  was  found  to  give  the  best  agreement  with  experiment,, 
and  this  value  is  used  in  Figure  21  (a  value  of  2  is  also  not  too  bad, 
although  not  as  good  as  4,  the  latter  value  allowed  most,  but  not  all.-  of 
the  data  to  be  used).  It  is  apparent  from  Figure  21  that  the  Gordon 
theory  with  a  modified  surface-tension  pressure  is  in  good  agreement 
with  experiment,  except  for  an  expected  constant  factor  (This  factor 
is  generally  larger  than  that  found  by  Equation  48a.  )  The  Gordon  theory 
reduces  to  a  simplified  breakup  time  expression,  Equation  25a,  for 
cases  where  viscosity  and  surface  tension  are  negligible.  This  expres¬ 
sion  is  of  interest  because  of  its  computational  ease. 

Figure  22  compares  the  theory  to  experimental  data  for  water  and  the 
iow-viscosity  silicone  fluid..  Although  the  agreement  is  generally  good, 
it  appears  that  the  effects  of  surface  tension  may  be  entering  in  slightly 
and  should  be  included.  The  breakup  theory  deiived  in  Section  4.  3  also 
reduces  to  the  Gordon  expression,  except  it  is  a  factor  of  2  smaller 
(see  Equation  50a);  thus,  it  also  is  in  satisfactory  agreement  with  the 
experimental  data.  The  limiting  expression  of  Gordon  for  large  viscosity 
and  negligible  suiface  tension  (Equation  25b)  is  identical  to  that  derived 
in  Section  4..  3  (Equation  50b). 

The  theory  is  compared  with  experimental  da’.a  for  the  high-visccsity 
silicone  fluid  in  Figure  23.:  A  comparison  of  the  results  with  figures  >6 
and  21  indicates  that  this  fluid  is  not  sufficient)!-  viscous  for  reliable  use 
of  the  simplified  equation. 

Hir.ze  developed  a  theory  of  the  breat-.j-  of  liqe.d  drops  for  the  (1  )  limit¬ 
ing  cases  of  negligible  viscosity  and  so  face  tension,  and  (2)  large  viscosity 
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,.nd  negligible  surface  tension.  His  expressions  for  the  breakup  times 
are  given  in  Equations  17  and  18.  These  expressions  are  compared 
with  the  experimental  breakup  times  of  water  and  the  Hoh-viscosity 
silicone  fluid  in  Figures  24  and  25;  the  agreement  appears  >o  be 
generally  less  satisfactory  than  for  the  Gordon  theory  and  that  devel¬ 
oped  m  this  report. 

As  summarized  in  Equation  32,  Morrell  developed  a  theory  for  the 
transverse  stripping  breakup  of  liquid  jets.  A  comparison  of  this 
theory  with  the  experimental  data  is  given  in  Figures  ’.6  and  27 
Figure  26  shows  that  there  is  a  drop- size  effect  not  apparently 
accounted  for  in  the  theory,  anti  both  illustrations  show  Ihu.  there  is 
considerable  deviation  between  the  experimental  data  and  the  theory. 

The  early  experimental  data  of  Lane  (Equation  13)  as  well  as  that  of 
Hanson  and  Domich  (Equation  19a)  showed  that  there  is  a  relationship 
between  the  critical  flow  velocity  required  to  break  the  drop  and  the 
diameter  of  the  drop.  These  relationships  are  plotted  in  Figure  28 
together  with  all  of  the  experimental  breakup  data  tabulated  as  to 
breakup  mechanism.  As  expected,  Llag  breakup  occurred  near  the 
critical  velocity  for  a  given  drop  diameter,  whereas  as  the  velocity 
was  further  increased,  the  breakup  behavior  exhibited  both  breakup - 
mechanism  characteristics.  At  higher  velocities,  only  the  shear 
breakup  behavior  was  observed. 


6.  3  DROPLET-SIZE  DISTRIBUTIONS 

Experimental  droplet-size  distributions  produced  by  the  aerodynamic 
breakup  of  bis  drops  are  given  in  Appendix  C.  Figure  29  shows  the 
mass-mean  diameters  of  the  droplets  produced  by  various  velocity 
gas  flows,  for  various  initial  drop  sizes,  and  for  both  bag  and  stripping 
breakup  behaviors. 

A  theory  for  the  average  droplet  size  produced  by  the  aerodynamic 
breakup  of  liquid  drops  was  developed  in  Section  4.4.  The  average 
droplet  size  is  given  by  Equation  58,  and  a  comparison  between 
experiment  and  theory  is  given  in  Figure  29.  Agreement  between 
experimental  data  and  theory  is  very  good,  considering  the  relative 
scatter  in  the  experimental  data.  In  view  of  this  scatter,  it  is  difficult 
to  draw  further  definite  conclusions  regarding  the  validity  of  Equation  58. 
Experimental  data  suggest  that  the  effect  of  the  initial  drop  diameter 
nr-'y  be  slightly  greater  than  the  1/6  power  dependence  predicted.  How¬ 
ever,  Weiss  and  Worsham  (Reference  19)  found  the  I  >6  power  for  the 
stripping  of  liquid  jets. 
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Figure  28.  Comparts on  of  Experimental  Results  from  Breakup  of  Water  wi'h  Critical  Size  Curves 
Obtain-. d  by  Hanson  and  Domich  (References  30,  31)  and  Lane  (Reference  28). 
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The  experimental  data  also  suggest  that  the  avc-rrge  drop  size  pro¬ 
duced  by  bag  breakup  may  not  quite  he  on  the  curves  tor  shear  breakup 
Again,  however,  it  is  difficult  to  draw  definite  concksions  but  the 
comparison  between  experiment  and  theory  is  generally  very  encour¬ 
aging  (further  experimental  data.,  however,  are  required  for  a  fin^l 
evaluation! 

.As  summarized  in  Equation  33.  Mayer  developed  a  theory  for  the 
average  droplet  size  produced  by  the  aerodynamic  breakup  o'  .iquid 
drops  A  comparison  of  the  theory  with  the  experimental  data  is  given 
in  Figure  30  A  value  of  B  ol  about  0  i  i=  required  to  b  the  theory 
to  the  same  magnitude  as  the  experimental  data,  but  the  functional  curia 
tion  still  does  not  appear  to  agree  with  the  experimental  results  The 
value  of  B  is  somewhat  arbitrary,  but  was  stated  to  be  conceptually  of 
the  order  of  unity  It  would  appear  that  the  Mayer  theory  is  less  ade  ¬ 
quate  than  the  theory  developed  in  this  report  (which  contains  no 
arbitrary  constants) 


6  4  DRAG-COEFFICIENT  MEASUREMENTS 

As  indicated  in  Section  4  3,  experimentally  determined  drag  loefficH  nts 
for  spheres  and  deformed  liqu’ds  being  accelerated  in  a  gas  stream  show 
a  wide  range  of  values 

The  results  of  drag-coefficient  measurements  made  on  3/32-in  -diameter 
nylon  spheres  and  some  of  the  liquids  used  in  this  study  are  shown  in 
Figures  31  through  40  In  these  illustrations  the  position,  celocity, 
acceleration,  lateral  diameter  (except  for  the  nylon  spheres*  and  drag 
co-fficient  of  the  sphere  or  drop  is  shown  as  a  function  of  tim»  after 
being  subjected  to  the  flow  In  each  case,,  the  sphere  or  drop  was  in 
free-fall  prior  to  onset  of  flow,  and  is  thus  depicted  as  having  a  slight 
negative  velocity  (assuming  the  air-flow  direction  to  he  positive)  before 
bong  accelerated  in  the  flow  direction 

Inspection  of  the  data  indicates  that  there  is  no  characteristic  shape  of 
the  drug-coeff ic ient -vs -time  curve  and  that  the  peak  values  of  the  curve 
may  vary  from  Cq  -  0  90  to  4  4  for  different  tests 

The  value  of  the  drag  coefficient  for  the  drop  or  sphere  undergoing 
acceleration  is  written  as  shown  in  Equation  61 
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Figure  iO.  Comparison  of  Experiments:’.  Average  Mass  Droplet 
Siae  with  the  Mayer  The<  ,  y. 
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Figure  31.-  Motion  of  a  3/32-in.  -Diameter  Nylon  Sphere  Subjected  to 
an  Airflow  of  208  ft/ sec. 
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Figure  32.  Motion  of  a  3/32 -in. -Diameter  N-,lo:i  Sphere  Subjected  to 
an  Airflow  of  386  ft/s  < c . 
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Figure  33.  Motion  of  a  3/32-in. -Diameier  Nyiin  Sphere  Subjected  to 
an  Airflow  of  434  ft/sec. 


-103- 


DISTANCE  ON) 


Figure  34.  Motion  of  a  1.  36- min.  -Distmete.  Drop  of  G.  E.  SF(96)  0.  65 
Silicone  Fluid  Subjected  to  an  Airilew  of  43  ft/ sec. 
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Figure  35.  Motion  of  a  1.25-min.  -Diar^.-Mer  Drc,-.  of  G.  E.  SF(96)  200 
Silicone  Fluid  Subjected  to  an  As  .low  of  161  ft/ sec. 
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Figure  36.  Motion  of  a  2.  2mm-Diameter  L>..  >p  of  G.  E.  SF{96)200 
Silicone  Fluid  Subjected  to  an  Airflow  •>£  208  ft/ sec. 
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Figure  37.  Motion  of  a  1.  45mm-Diameter  Drop  of  G.  E.  SFj96)Z00 
Silicone  Fluid  Subjected  to  an  Airfli  *  of  3  40  ft/ sec. 
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Figure  3f.  Motion  of  a  1.  2mm-Diamei .  x  '  rop  ot  Mercury 
Subjected  to  an  /irflow  of  168  It 'sec. 
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Figure  39..  Motion  o£  a  0.  69mn-Diamf  ter  Drop  of  Mercury 
Subjected  to  an  Airflow  of  2d*  tt/ser. 
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dv 

dt 


liquid  or  solid  density- 
air  density 
air  velocity 

velocity  of  sphere  or  drop 
original  drop  or  sphere  diameter 
lateral  diameter  or  deformed  drop 
acceleration  of  sphere  or  drop 


Equation  61  indicates  that  the  drag  coefficient  is  directly  proportional 
to  the  acceleration  and  inversely  proportion?*  to  the  square  of  the  rela¬ 
tive  velocity  and  the  square  of  the  diameter  of  the  drop  normal  to  the 
airflow  direction. 

A  comparison  of  the  drag-coefficient  curves  to  the  corresponding  accel¬ 
eration  curves  indicates  that  the  drag  coefficient  is  influenced  more  by 
acceleration  changes  (which  are  very  great)  than  by  changes  in  relative 
velocity  and  lateral  diameter.  This  effect  leads  to  the  possibility  of 
rather  large  uncertainties  in  the  calculated  values  of  the  drag  coefficient,, 
since  the  acceleration  is  determined  by  a  numerical  differentiation  of 
the  velocity  curve  (and  the  velocity  curve  is  obtained  in  a  like  manner 
from  the  dis’ance  or  position  curve)..  The  errors  arise  from  the  pro_e.» 
of  taking  the  difference  between  two  value®  that  are  close  together  (close- 
spacing  of  the  readings  is  necessary  to  retain  reasonably  accuratt  rep¬ 
resentation  of  the  curves). 

Attempts  to  eliminate  these  errors  by  utilizing  digital  computer  tech¬ 
niques  to  generate  best-fit  curves  to  the  raw  data  and  compute  di  ae 
coefficients  were  unsuccessful  bee  ••  <•  complex  expressions  were 
necessary  to  yield  a  satisfactory  fit  ov..  r  the  range  of  data..  Hand- 
calculation  techniques  were  then  resorted  to. 


-Ill- 


Although  uncertainties  in  the  drag-coefficient  values  are  present,,  it  is  not 
felt  that  they  are  of  sufficient  magnitude  to  account  for  the  -.’ide  range 
cf  peak  values,  nor  to  account  for  the  widely  differing  shapes  of  the  drag- 
coefficient  curves. 

It  thus  appears  that  no  simple  expression  can  be  found  to  predict  drag- 
coefficient  values  for  small  spheres  or  liquid  drops  being  accelerated  in 
a  moving  gas  stream.. 


7.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  studies  described  in  this  report  have  encompassed  a  considerable 
amount  of  literature  work,  original  theoretical  studies,  and  unique 
experimental  studies.  The  following  discussion  summarizes  major 
points  of  these  studies  and  provides  recommendations  for  future  work. 

The  studies  have  shown  that  the  breatcup  of  a  liquid  is  a  rate  process 
and  should  be  considered  from  a  time- dependent  point  of  view.  This 
consideration  is  particularly  important  when  high-speed  gas  flows  of 
transient  nature  are  involved  (e.  g. ,  shock  loading),  since  the  gas-flow 
acting  on  the  liquid  may  be  a  rapidly  changing  function.  It  is  also 
suggested  that  the  aerodynamic  breakup  ofbulk  liquid  with  any  geometry 
may  be  considered  to  occur  in  two  main  phases.  The  first  phase  com¬ 
prises  the  shearing  of  layers,  sheets,  films,  or  jets  from  the  liquid  by 
the  gas- flow  forces;  the  second  comprises  the  rapid  breakup  of  these 
liquid  streams  by  instabilities  that  grow  exponentially  with  time..  An 
estimate  of  the  aerodynamic  breakup  time  of  the  liquid  may  be  made  by 
quantitatively  treating  the  first  phase,  whereas  the  resultant  droplet 
sizes  may  be  estimated  by  quantitatively  treating  the  second. 

The  preceding  concept  v'as  applied  to  liquid  drops,  anu  theoretical 
expressions  were  derived  for  both  the  breakup  times  of  the  original  drop 
and  the  resultant  droplet  sires  produced  by  breakup.  Experimental 
breakup  times  of  drops  of  a  variety  of  liquids  with  a  wide  range  of 
physical  properties  were  determined  ior  several  drop  sizes  over  a  wide 
range  of  gas-flow  velocities.  These  studies  were  conducted  using  shock- 
tube  and  high-speed  photographic  techniques;  agreement  between  experi¬ 
mental  data  and  theory  was  very  good,  and  it  n-.ay  be  concluded  that  the 
breakup-time  theory  developed  in  this  report  is  valid  and  adequate  for 
general  use.  It  would  be  of  interest,,  however,  to  extend  thv  experimental 
studi'  3  to  much  higher  gas-flow  velocities  and  much  smaller  initial  crop 
sizes.. 
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The  experimental  breakup  times  were  also  compared  to  the  available 
theories  in  the  literature.  It  was  found  that  there  is  good  agreement 
between  the  Gordon  theory  and  experiment  if  the  surface-tension 
pressure  used  in  the  Gordon  theory  is  decreased  by  a  factor  of  2  to  4.. 
The  breakup  theories  of  Hinre  and  Morrell  were  found  to  be  generally 
unsatisfactory  for  use. 

Liquid  drops  are  observed  to  break  by  either  a  bag  or  stripping  mech¬ 
anism,  and  it  was  postulated  that  bag  breakup  is  due  to  pressure  arag, 
whereas  stripping  breakup  is  due  to  friction  drag  (friction  drag  was 
considered  to  be  always  greater  then  *hc-  pressure  H-a^'  A  theory  w\..s 
presented  showing  that  the  original  drop  deformation  produced  by  the 
gas  flow  increases  the  curvature  and  hence  the  surface-tension  pressure 
at  the  outer  edges  of  the  drop.  This  effect  explains  why  bag  breakup 
occurs  near  the  critical  velocity  required  for  breakup,  whereas  strip¬ 
ping  breakup  occurs  at  higher  velocities.:  The  experimental  results  are 
generally  in  agreement  witli  the  theory.  It  would  be  of  interest  to  extend 
the  theoretical  studies  so  that  a  more  rigorous  treatment  may  be 
developed. 

In  the  development  of  the  model  for  the  aerodynamic  breakup  of  liquid 
drops,  it  was  pointed  out  that  with  certain  modifications,  the  treat¬ 
ment  should  also  be  suitable  to  describe  the  breakup  (erosion)  of  solid 
particles  in  high-velocity  gas  flows.  It  would  be  desirable  to  conduct 
both  experimental  and  theoretical  studies  along  these  lines  (these 
studies  may  have  direct  applicability  to  the  problem  of  the  breakup  of 
particle  agglomerates). 

Preliminary  experimental  studies  were  conducted  to  obtain  the  droplet-- 
size  distributions  produced  by  the  aerodynamic  breakup  of  liquid  drops, 
and  the  experimentally  determined  mass-mean  diameters  of  the  result¬ 
ant  drops  were  compared  with  the  theory  that  was  developed  in  this 
report.  Agreement  was  very  good,  considering  the  relative  scatter  in 
the  experimental  data  and  the  lact  that  the  theory  contains  no  arbitrary 
constants.-  The  experimeatal  data  were  also  compared  to  the  theory  of 
Mayer;  agreement  was  less  satisfactory  than  with  the  theory  developed  in 
this  report,  and  moreover  contains  a  constant  that  must  be  evaluated 
from  the  experimental  data.-  It  is  concluded  that  the  developed  theory  is 
at  present,  most  reliable  for  use..  As  a  final  recommendation,  it  is 
suggested  that  further  studies,  both  experimental  and  theoretical,  he 
conducted  w:th  regard  to  obtaining  the  droplet-size  distributions  pro¬ 
duced  by  the  aerodynamic  breakup  of  liquid  drop3.. 
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APPENDIX  A 

SHOCK-TUBE  THEORY 


A- 1  INTRODUCTION 


The  following  discussion  begins  with  a  qualitative  desc  iption  of 
shock-tube  behavior.  For  simplicity,  the  description  is  limited 
to  the  Darticular  case  where  air  is  used  in  both  the  compression 
and  expansion  chambers  and  the  expansion  chamber  is  open-ended. 

As  shown  in  Figure  A- la,  the  gases  in  the  compression  and  expan¬ 
sion  chambers  are  at  the  same  temperature  but  at  different  pres¬ 
sures,  with  P3  being  greater  than  PQ.-  Since  the  expansion  chamber 
is  open-ended,  the  pressure  and  temperature  therein  are  •'  s-'imed 
to  be  at  local  ambient  values.  Pressure  P3  is  maintained  at  a  value 
greater  than  P0  by  an  airtight  diaphragm  separating  the  two  chambers. 

If  the  diaphragm  separating  the  two  regions  is  suddenly  removed,  a 
compression  wave  progresses  into  the  expansion  chamber,  while  a 
rarefaction  wave  simultaneously  moves  toward  the  closed  end  of  the 
compression  chamber.  After  some  distance,  the  compression  wave 
steepens  into  a  plane  shock  wave,  which  moves  through  the  expansion 
chamber  at  a  velocity,  w,  determined  by  the  pressure  ratio, 

P3/P0  =  **30-  Passage  of  the  shock  wave  imparts  a  uniform  velocity, 
u,  to  the  expansion- chamber  gas  and  raises  its  pressure,  density,- 
and  temperature  to  Pj,  pj,  and  Tj,  so  that  PQ  <  Pj  <r  P3, 
pQ  <  pj  <  p,,  and  Tj  >  TQ..  The  pressure-distance  conditions  are 
depicted  in  Figure  A- lb,  where  the  shock  wave  has  reached  the  test 
location  and  the  rarefaction  wave  is  progressing  toward  the  closed 
end  of  the  compression  chamber. 

The  contact  surface  (sometimes  called  the  "cold  front  !  represents  the 
boundary  between  the  compression  and  expansion- chamber  gases. 

It  also  moves  through  the  expansion  chamber  with  the  velocity  u,  as 
does  the  compression  chamber  gas  at  Pj  behind  the  contact  surface. 

The  region  behind  the  contact  surface  is  compression- chamber  gas 
that  is  expanded  down  to  pressure,  temperature,  and  density  Pj,  T £, 
p,  by  the  rarefaction  wave  so  that  PC1  <  P,  <  P  ,  T  <  T0,  and 

p,  <  p2  <  P3- 

In  Figure  A-lc,  the  shock  wave  is  shown  as  having  arrived  at  the 
open  end  of  the  expansion  chamber;  the  contact  surface  has  continued 
to  move  through  the  expansion  chamber  3t  the  velocity  u,  and  the 
rarefaction  wave  has  collided  with,  and  rebounded  from,  the  closed 
end  of  the  compression  chamber. 
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Figure  Al.  Shock  Tube  Pressure  Diagrams. 
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Th"  exit  of  the  shock  wave  from  the  open  end  of  the  expansion  cham¬ 
ber  creates  a  rarefaction  wave,  which  moves  back  onto  the  expansion 
chamber  as  shown  in  Figure  A -Id,  where  it  has  reached  the  test 
location. 

The  region  between  the  shock  wave  and  the  contact  surface,  being  a 
zone  of  constant  pressure,  temperature,  density,  and  particle  veloc¬ 
ity,  is  ideally  suited  as  an  aerodynamic  test  region  of  short  duration.. 
Testing  time  available  is  the  elapsed  time  from  the  passage  of  the 
shock  wave  to  the  arrival  of  one  of  the  following:' 

•  The  contact  surface 

•  The  reflected  rarefaction  from  the  compression  chamber 

•  The  rarefaction  wave  from  the  open  end  of  the  expansion 
chamber 

In  Figure  A-l,  the  testing  time  would  be  represented  by  the  time  for 
the  events  from  (b)  to  (d)  to  occur.  For  this  particular  case,  the 
testing  time  is  terminated  by  the  arrival  of  the  rarefaction  from  the 
open  end  of  the  expansion  chamber. 


A -l  SHOCK-TUBE  THEORY 


The  equations  describing  the  conditions  in  the  shock  lube  are  derived 
by  making  use  of  the  equations  of  continuity,  momentum,  motion,  and 
energy,  plus  the  equation  of  state  for  ideal  gases..  Assumptions  are 
made  that  the  flow  is  inviscid,  that  shock  compression  occurs  adia- 
batically,  and  expansion  is  isentropic.; 

One  of  the  most  important  expressions  in  shock-tube  behavior  is  that 
showing  the  relationship  between  the  pressures  in  the  various  regions. 
Adopting  the  subscript  notation  of  Figure  A-l  and  writing  the  pressure 
ratio  Pj/Pq  as  Pjq,  the  relationship  between  the  compression-chamber 
pressure  and  the  shock  pressure  can  be  written  as 
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where 


P 3q  =  ratio  of  compression-chamber  p-essure  to  expansion 
chamber  pressure 
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ratio  of  pressure  behind  shock  wave  to  expansion- 
chamber  pressure  ahead  of  sho-k  wave 


Cp/Cy  =  specific  heat  ratio  for  the  gas 


In  Equation  A-l  and  the  following  expressions,  it  has  been  assumed 
that  y  remains  constant  and  has  the  same  value  in  all  of  the  flow  regions. 
Work  reported  in  Reference  68  substantiates  this  assumption  for  the 
normal  operational  range  of  shock  tubes  similar  to  those  u-  ed  in  the 
present  study. 


The  speed  of  the  shock  wave  is  given  by 
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where 


=  speed  of  sound  in  the  gas  ahead  of  the  shock 
The  flow  velocity  behind  the  shock  wave  is  given  by 
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The  density  ratio  across  the  shock  wave  is  written  as 
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Using  the  ideal  gas  law,  the  temperature  ratio  across  the  shock  wave 
becomes 
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A  complete  derivation  of  Equations  A-l  through  -5  can  be  found  in 
numerous  sources  available  in  the  open  literature  (References  66 
through  70  and  73). 

To  illustrate  the  variation  of  the  parameters  given  in  Equations  A-l 
through  -5  with  change  in  compression-chamber  pressures,;  values  were 
calculated  for  various  pressures  used  in  this  study.  The  results  of  the 
calculations  are  shown  in  Figure  A-2,  where  the  initial  conditions  in  the 
expansion  chamber  were  assumed  as  follows: 


p  = 

o 

14.  25  psia 

T  = 
o 

535  °R 

0.  0719  lb/ft 

ao  = 

1 1  35  ft/ sec 

The  shock  tube  used  in  this  study  was  located  at  an  elevation  of  approxi¬ 
mately  1000  ft  above  sea  level  in  an  air-conditioned  room,  with  the 
temperature  maintained  at  an  average  value  of  75°F.  Thus,  the  con¬ 
ditions  used  for  the  calculations  are  typical  of  values  encountered 
during  the  tests 

Figure  A-3  shows  the  variation  of  pressure,  temperature,  and  densit\ 
behind  the  shock  wave  with  flow  velocity  for  the  range  covered  in  this 
study. 

Duration  of  the  testing  time  was  calculated  from  the  various  expressions 
describing  flow  conditions  in  the  shock  tube  (Equations  A-l  through  -5) 
and  the  lengths  of  the  various  portions  of  the  tube.  To  compute  the 
testing  time  available,  the  following  nomenclature  was  employed- 

Lv  =  distance  from  diaphragm  to  test  location 

L(.  -  distance  from  diaphragm  to  end  of  .  xpansion  chamber 
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Figure  A2.  Shock  Tub  *  Flow  Condition*  as  a  Function  of 
Compression  Chamber  Pressure 
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length  of  compression  chamber 

time  for  shock  wave  to  arrive  at  test  location 

time  for  cold  front  to  arrive  at  test  location 

time  for  foot  of  rarefaction  to  arrive  at  test  location 

time  for  head  of  rarefaction  to  catch  its  foot 

time  for  head  of  rarefaction  to  overtake  cold  front 

time  for  head  of  rarefaction  to  arrive  at  test  location 

distance  for  head  of  rarefaction  to  catch  its  foot 

distance  for  head  of  rarefaction  to  overtake  cold  front 

time  for  rarefaction  head  from  open  end  of  expansion 
chamber  to  arrive  at  test  location 


All  other  nomenclature  is  identical  to  that  used  in  Equations  A-l 
through  -5  and  Figure  A-l . 

The  time  required  for  the  shock  to  arrive  at  the  test  location  is 
given  by 


t  =  .^A/ia  4  m 

8W  aoV*P10  +  1 


(A-6) 


The  time  required  for  the  contact  surface  (cold  front)  to  arrive  at 
the  test  location  is 
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The  time  required  for  the  foot  of  the  rarefaction  to  arrive  a.,  the  test 
location  is 
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where  a,  =  a  if  T  ,  =  T  as  assumed  for  this  study. 

3  o  3  o 

The  time  required  for  the  head  of  the  rarefaction  to  overtake  its  foot 
is 
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Also,  it  can  be  shown  that 
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The  time  required  for  the  rarefaction  head  to  arrive  at  the  test  loca¬ 
tion  (if  the  cold  front  arrives  first)  is  given  by 
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The  time  required  for  the  rarefaction  head  to  arrive  at  the  test  loca¬ 
tion  (if  the  head  has  overtaken  the  cold  front  and  arrives  first)  can  be 
expressed  as 


t  =  t  + 
rw  rc 


L  -  X 


w 
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rc 
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where  a  ,  =  a 
1  o 

The  distance  required  for  the  head  of  the  rarefaction  wave  to  over¬ 
take  its  foot  is  given  by 


Xrf=  (a  -  l)a3trf 
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The  distance  required  for  the  head  of  the  rarefaction  wave  to  over¬ 
take  the  contact  surface  is 


X  =  X  ,  + 
rc  rf 
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The  time  required  for  the  head  of  the  rarefaction  from  the  open  end 
of  the  expansion  chamber  to  arrive  at  the  test  location  is  given  by 
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The  time-of-ar rival  of  the  various  phenomena  described  in 
Equations  A-6  through  -14  is  shewn  in  Figure  A-4  for  the  ra-  ge  of 
chamber  pressures  used  in  this  study  Figure  A-5,,  obtained  from 
Figure  A-4,  shows  directly  the  testing  time  available  for  the  same 
range  of  compression- chamber  pressures 


A-2  SHOCK-TUBE  PERFORMANCE 

In  all  but  the  most  precise  shock-tube  investigations,  it  i3  customary 
to  determine  the  flow  conditions  behind  the  shock  wave  eitiiei  from 
the  slate  of  the  expansion  and  compression- chamber  gases  before 
rupturing  the  diaphragm,  or  fiom  the  state  of  the  expansion- chamber 
gas  ahead  of  the  shock  wave  and  a  direct  measurement  of  shock  velocity . 
The  latter  method,  which  yields  a  more  accurate  prediction  of  flow 
conditions  than  the  former,  is  the  procedure  used  in  this  study. 

The  sources  of  disagreement  between  theory  and  actual  behavior  in 
the  shock  tube  are  primarily  the  result  of  departure  from  the  assump¬ 
tions  of  adiabatic  compression,  isentropic  expansion,  and  constant 
specific-heat  ratio,  and  to  losses  due  to  boundary-layer  effects. 

Variation  of  the  specific-heat  ratio  has  been  investigated  (Reference 
68);  the  results  indicate  that  for  shock  velocities  below  about  Mach  3, 
the  differences  between  using  constant  and  variable  specific-heat 
ratios  are  not  noticeable  The  assumptions  of  adiabatic  compression 
and  isentropic  expansion  are  also  considered  quite  good  over  this 
range  of  shock  velocities. 

Since  this  study  was  conducted  at  shock  velocities  beiow  Mach  2,:  it 
would  be  expected  that  the  results  should  agree  quite  well  with  theory 
This  assumption  is  verified  in  Figure  A-6,  where  experimental  shock 
velocities  in  the  shock  tube  used  ir  this  study  are  compared  with 
theory.  Assuming  no  error  in  shock-velocity  measurement,  the 
experimental  values  show  a  maximum  departure  of  1. 6%  from  the 
theoretical  values. 

Errors  in  shock-velocity  determinations  are  caused  by  errors  in  the 
measurement  of  the  time  required  for  the  shock  wave  to  traverse  a 
known  distance  and  in  the  measurement  of  this  distance.  The  distance 
between  the  transducers  used  in  the  shock  tube  was  1*  25  in.  (accurate 
to  .vithin  -  1  / 64  in  ),  and  the  counter  used  to  record  the  time  required 
for  the  shock  wave  to  traverse  the  distance  between  the  transducers 
was  calibrated  to  the  nearest  0.  1  psec 
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The  primary  source  of  error  in  shock-velocity  measurements  is 
associated  with  the  signal  output  from  the  transducers:  this  output 
does  not  rise  instantaneously  to  its  peak  value  with  the  at  rival  of 
the  shock  wave,  and  it  is  difficult  to  adjust  both  the  start  and  stop 
circuits  of  the  counter  to  trigger  at  exactly  the  same  voltage  level. 

At  the  start  of  the  study,  measurements  were  made  of  the  wave 
shape  of  the  output  signals  from  the  two  transducers.  The  results 
showed  the  waves  to  be  nearly  identical  in  shape  over  the  range  of 
interest.  Typical  output  signals  from  the  transducer  used  to  stop 
the  counter  are  shown  in  Figure  A- 7:  superimposed  on  the  curves 
are  the  voltages  used  to  adjust  the  start  and  stop  circuits  of  the 
counter.  Voltages  equal  to  or  less  than  the  lower  "NO  GO"  trigger 
level  failed  to  trigger  the  circuit,  while  signals  equal  to  or  greater 
than  the  upper  "GO"  trigger  level  resulted  in  triggering.  The  differ¬ 
ence  between  these  levels  then  represents  the  maximum  uncertainty 
in  trigger  voltage  between  the  start  and  stop  circuits,  and  a  corres¬ 
ponding  uncertainty  in  the  recorded  time  is  thus  determined. 

The  total  uncertainty  in  shock-velocity  measurements  is  represented 
by  the  combination  of  errors  in  distance  and  time  measurements. 
Assuming  the  worst  combination  of  these  effects  to  occur  simultane¬ 
ously,'  a  maximum  error  of  t  0.  85%  in  shock  velocity  could  be  expected 
for  the  range  of  velocities  covered  in  this  study 

Using  shock- tube  theory  to  obtain  a  relation  between  shock  velocity 
and  air  velocity  behind  the  shock  wave,  a  curve  similar  to  that  shown 
in  Figure  A-8  is  obtained.  Applying  the  possible  errors  in  shock- 
velocity  measurements  to  the  illustrated  curve  results  in  a  possible 
error  of  approximately  i  28%  for  an  air  velocity  of  60  ft/sec,,  while 
the  same  technique  >ields  an  error  of  only  -  3%  for  an  air  velocity  of 
450  ft/sec.  Inspection  of  the  curve  in  Figure  A-8  shows  that  while 
the  accuracy  of  shock-velocity  measurement  is  practically  constant, 
large  errors  in  air  velocity  can  be  expected  at  shock  velocities  near 
the  speed  of  sound,  since  the  a>r  velocity  approaches  zero  as  the  shock 
velocity  approaches  the  speed  of  sound  ahead  of  the  shock  wave 

It  should  be  mentioned  that  while  the  preceding  discussion  implies  the 
possibility  of  rather  large  errors  in  the  air  velocities,  it  is  not  likely 
that  errors  of  this  magnitude  are  present  in  the  data  With  careful 
attention  given  to  the  adjustment  of  the  trigger  levels  or.  the  counter,; 
the  operator  could  soon  develop  a  technique  of  reaching  settings 
extiemely  close  to  the  lower  trigger-level  voltage.  A.s  a  result,  it 
is  improbable  that  errors  in  air-velocity  measurements  as  large  as 
those  derived  from  Figure  A-8  actually  .  .  -red.. 
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Figure  A7.  Transducer  Output  Signals. 
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APPENDIX  B 


This  appendix  contains  plates  made  from  the  high  spt  ed 
motion  picture  films  taken  durine  the  study  The  sene-> 
of  plates  contained  herein  represent  each  of  the  test 
conditions  listed  in  Table  2  in  Section  5  i  of  this  report 
They  have  been  included  to  provide  a  source  of  com¬ 
prehensive  information  on  the  qualitative  aspects  of 
drepiet  breakup  as  well  as  detailed  quantitative  data  on 
the  behavior  of  drops  during  deformation  acce'erat.on 
and  subsequent  breakup  m  a  fast  moving  air  stream 
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APPENDIX  C 


The  results  of  the  size  distribution  tests  (summarised  in  Table  C-l) 
are  presented  in  this  appendix  in  the  three  following  forms:' 

a.  Tabulations  of  assessment  data  taken  from  photomicrographs  of 
the  sampled  aerosols. 

b  Plates  made  from  photomicrographs  of  a  selected  number 
of  the  tests,. 

c.  Plots  of  cumulative  percent  mass  vs  droplet  size  for  all  of  the 
data  in  (a)  above. 

The  tabulated  assessment  data  has  been  included  for  reference  purposes 
and  is  more  detailed  in  nature  than  that  presented  in  other  sections  of 
the  report. 

The  plates  have  been  included  to  show  visually  the  reduction  in  particle 
size  that  occurs  with  increase  in  relative  velocity. 

The  cumulative  percent  mass  vs  particle  diameter  curves  show  the  shape 
of  the  size  distribution  curves,  and  were  used  to  determine  the  MMD 
values  used  for  correlation  with  theory  in  other  sections  of  the  report 


Table  C-l.;  Summary  of  Size  Distribution  Tests. 


Shear  Breakup 


Orig 

Drop 

Size 

iramt 

V 

(ft/sec) 

Sample 
Dist 
(in.. ) 

Film 

No. 

Test 

No. 

MMD 

(pi 

2. 

,  7 

138 

15 

926 

35 

85 

1 

110 

18 

850 

6 

90 

112 

24 

849 

5 

120 

234 

27 

Q£7 

41 

60 

1 

212 

30 

865 

12 

73 

212 

30 

872 

12 

71 

1 

i 

200 

36 

845 

1 

59 

! 

315 

27 

928 

44 

32 

1 

l 

319 

30 

866 

13 

36 

319 

30 

873 

13 

41 

395 

24 

874 

18 

34 

383 

27 

929 

48 

19 

394 

30 

867 

14 

24 

1 

394 

30 

875 

14 

24 

1. 

6 

131 

12 

879 

22 

94 

1 

130 

18 

878 

21 

75 

132 

24 

880 

23 

93 

201 

24 

859 

10 

45 

206 

30 

882 

24 

74 

i 

202 

36 

858 

9 

53 

! 

316 

24 

885 

26 

21 

i 

320 

30 

886 

27 

30 

} 

j 

395 

24 

891 

30 

18 

396 

30 

912 

31 

18 

394 

36 

913 

32 

15 

Bag  Breakup 

0. 

6 

69 

4- 3/8  in 

930 

55 

100 

0. 

6 

66 

4-3/8  in. 

931 

56 

140 

0. 

6 

4-3/8  in. 

932 

57 

115 

*  Sample 

Distance  with  1-ft  T>ection  Removed 
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NOMENCLATURE 


L.  P. 

M.  F. 
U.  P- 
FO 
PFO 
SFO 
SPFO 
PMO 
SPMO 


lower  point  of  size  interval,  microns 
mid-point  of  size  interval,  microns 
upper  point  of  size  interval,  microns 
frequency  observed 
percent  frequency  observed 
cumulative  frequency  observed 
cumulative  percent  frequency  observed 
percent  mass  observed 
cumulative  percent  mass  observed 


-C3- 


TEST  NO.  ?>5  DROP  SIZE  (MM)  2.7 
RELATIVE  VELOCITY  (FT/SEC)  138 
FILM  NO.  ??.£■  MMD  (m)  85 
SAUPt  ING  DISTANCE  (IN.)  i  5 
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Cl 
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3.19 
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1  .44 
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'.4. 

1.1? 
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.'5 
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• 
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a  -  • 1  a 
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.  •  “ 

14  7?. 

7-.  51 

1.02 

*4  S  .  . 

~  t  ~  *♦ 

<“  3 

✓  i 

.  1  ' 

•  9  74. 

7-  •  74 

2.43 

c  3  «  a  A 

9  '3 

£ 

.  *: 

14’4, 

4  9*0- 

*  •:•*<• 

-  «  *  1 
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?# 

.i? 
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9  '/  .  1  - 

'  ♦ 

t  .  • 

7  ?  .  7  4 

IT).!  J 

!  • 

•  0; 
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i  «  *  * 
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• ;  • 
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DROP  SIZE  (MM)2.  7 

110 


L.  P. 

M.  P. 

0.  P. 

0.00 

3.01 

6.02 

6.32 

8.70 

11.39 

11.39 

14.08 

16.77 

16.77 

19.99 

23*22 

23.22 

25.58 

27.95 

27.95 

30.63 

33.3? 

33.32 

36.01 

36.70 

38.70 

41.49 

44.2« 

46.29 

47.08 

49.  ae 

49.88 

52.56 

55.25 

55.25 

57.94 

60.63 

60.63 

63.42 

66.2? 

66.22 

68.58 

73.95 

70.95 

73.6? 

76.32 

76.32 

78.58 

80.84 

80.84 

83.63 

86.43 

86.43 

68.58 

90.73 

90.73 

93.31 

95.89 

95.89 

93.04 

100.19 

100.19 

101.69 

IC3.20 

103.20 

108.8? 

114.5? 

1 1':-59 

117.28 

119.97 

119.97 

12  3.19 

126.42 

126.42 

128.78 

131.15 

131.15 

133.83 

136.5? 

136.52 

139.21 

141,90 

141.90 

144.69 

14  ?  .4  9 

147.49 

150.28 

153.08 

153. C8 

155.76 

158.45 

TEST  NO.  6 
RELATIVE  VELOCITY  (FT/SEC) 
FILM  NO.  350  MMO  W  90 
SAMPLING  DISTANCE  (IN.)  18 


FO 

PFO 

SFO 

0. 

0 .00 

0. 

210. 

24.59 

210. 

193. 

22.59 

403. 

166. 

19.43 

569. 

65. 

7.61 

634. 

74. 

8.66 

708. 

29. 

3.39 

73  7* 

22. 

2.57 

759. 

19. 

2.22 

778. 

19. 

2.22 

797. 

5. 

.58 

802. 

9. 

1.05 

811. 

10. 

1.17 

821. 

9. 

1.05 

830. 

3. 

.35 

833. 

2. 

.23 

835. 

3  . 

.35 

838. 

0. 

3.00 

838. 

4. 

.46 

842. 

2. 

.23 

844. 

1. 

.11 

d45. 

3. 

.35 

848. 

1. 

.11 

849. 

2. 

.23 

851. 

2. 

.23 

853. 

0. 

o.oo 

853. 

0. 

3.00 

853. 

0. 

3.00 

S53. 

1. 

.11 

854. 

SPFG 

PMw 

S8WG 

0.00 

0.00 

0.00 

24.59 

.32 

.32 

47.18 

1.06 

1.41 

66.62 

2.50 

4.02 

74.23 

2.03 

fc.05 

o2.9D 

3.96 

10.02 

86.29 

2.50 

12.53 

69.87 

2. 90 

15.43 

91.10 

3.6  5 

19.08 

93.32 

5.06 

24.15 

93.91 

1.78 

25.54 

94.96 

4.20 

30.15 

96.13 

5.90 

36.05 

97. IB 

6.57 

42.62 

97.54 

2.65 

45.28 

97.77 

2.13 

47.42 

98.12 

3.60 

7  1.23 

98.12 

0.00 

61.23 

VS.  5? 

6.6  7 

.6.11 

98.82 

3.83 

61.94 

99.94 

2.37 

t4.31 

99.29 

6.03 

I  i,\ 5 

99.41 

3.41 

/B.VO 

99.64 

7.79 

6  4.35 

99.88 

8.  74 

43.10 

99.88 

0.00 

o 

99.88 

0.00 

<3.10 

99.88 

0.00 

93.10 

H'0.00 

6.69 

100. OC 
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TEST  NO.  5  DROP  SIZE  (MM)  2.  7 
8ELATIVE  VELOCITY  <FT/SEO  112 
F1LVNO.S49  UMDOi)  120 
SAMPLING  DISTANCE  ON.)  24 


L  «  0. 

4.  o. 

U.  P. 

FO 

PFC 

SFO 

5PF0 

PMO 

■7.00 

3.01 

6.02 

0. 

0.00 

0. 

0.00 

0.00 

0.00 

6.02 

8.70 

11.39 

61. 

17.28 

61. 

17.28 

•  10 

.  10 

11.35 

14.06 

16.77 

61. 

17.28 

122. 

34,56 

.39 

.30 

16.77 

19.99 

23.22 

52. 

14.73 

174. 

49.29 

.93 

l.<*3 

23.22 

28.58 

27.95 

28. 

7.93 

202. 

57.22 

1.00 

2.4? 

27.95 

30.63 

33.32 

37. 

10.48 

239. 

67.70 

2.26 

4.70 

33.32 

36.01 

38.73 

22. 

6.23 

261. 

73.93 

2.17 

6.89 

38.70 

41.49 

44.?9 

23. 

6.31 

284. 

BO. 45 

3.46 

10*37 

44. 29 

47. oe 

49.88 

12. 

3.39 

296. 

d3.85 

2.63 

12.96 

49.98 

52.56 

55.25 

20. 

5.66 

316. 

69.51 

6.09 

55.25 

57.94 

60.63 

6. 

1.69 

322. 

91.21 

2.44 

21.33 

60.63 

63.42 

66.22 

3. 

.84 

325. 

92.06 

1.60 

*  ?*  1  3 

66.22 

68 . 58 

70.95 

3. 

.84 

328. 

92.91 

2.02 

23*16 

70.95 

73.67 

76.32 

?. 

.56 

330. 

93.48 

1.6  7 

2  6  •  H  ? 

76.  72 

78.58 

80.84 

3. 

.84 

333. 

94,33 

*.04 

29.67 

60.04 

83.63 

86.43 

4. 

1*13 

3  3  7. 

95.46 

4.89 

34.7 1> 

66.4  3 

88.58 

90.73 

4. 

1.13 

341. 

96.60 

5.60 

40.36 

90.73 

97.31 

95.89 

2. 

.56 

343. 

97.16 

3.39 

4  3.93 

96.89 

98.04 

100.19 

C  . 

9.00 

34  3. 

97,16 

0.03 

4  3.97 

130.19 

101.69 

103.20 

1. 

.26 

344* 

97.45 

2*  1 

«  J  7* 

1 33.20 

ICS. 89 

114.59 

1. 

.28 

345. 

97.73 

2.71 

**  8  .  *•  t* 

1 14.59 

117.28 

1  19.97 

0. 

•0.00 

745. 

97.  93 

o.n  0 

43.6s 

] 19. 9? 

123.19 

126.42 

1. 

.28 

046. 

98.01 

3.  VO 

;»  2 « 76 

126.42 

128. 78 

131.15 

1. 

.28 

347. 

98.30 

4*45 

>  ?  *  2 1 

131.15 

133. 6* 

136.52 

1. 

.26 

348. 

98.58 

3*00 

6  <r  *  2  1 

136.52 

1 39.21 

141.90 

2. 

.56 

350. 

99.15 

11  •  4  •_> 

/  £  «  &  7 

141.90 

144.69 

14  /.49 

1* 

.28 

351. 

99.43 

6*  l  1 

79*73 

147.49 

150.28 

153.08 

0. 

0.00 

351. 

99,43 

0*00 

7  V  #  7fy 

153.08 

155.76 

159.45 

0. 

•">.U0 

751. 

99.43 

0.00 

79*78 

1 58.45 

161.14 

163.83 

C. 

o.oc 

351. 

99.43 

0*00 

7#.  7  0 

1 63. 8  3 

166.62 

169.42 

1. 

.2d 

352. 

99.71 

t  • 

li  T  ♦  ^ 

1 69.42 

171.76 

1  74.15 

1. 

.2d 

35  3. 

100.00 

10.36 

;o  3. m 
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TEST  NO.  4 1  DROP  SIZE  (MM)  2.  7 
RELATIVE  VELOCITY  (FT/SEC)  234 
FILM  NO.  927  MMD  (p)  60 
SAMPLING  DISTANCE  (IN.)  27 


4  '  • 

"«  f  • 

a. 

?o 

PFO 

SFQ 

SPFO 

PMO 

SPMO 

• 

•  •  .'  I 

7  .'  > 

o. 

o..;o 

r. 

C.00 

0.00 

0.00 

4 

3.7 

li.J9 
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.75 
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16.  1£ 

16.  •  1 

182. 

26.09 

379. 

fS.?9 

2.59 

3.29 

Itj  •  :  t 

1«.95 

4  3.2? 

100. 

15.90 

9  79. 

73.80 

3.  ’’ 

7.16 

23.22 

8?.  66 

2  7.96 

36. 

3.69 

515. 

79.35 

2.7  8 

9.9  7 

>7.95 

3 '  .63 

37.72 

29. 

9.96 

599. 

83.82 

3.89 

13.82 

76.31 

16.70 

22. 

7.36 

6f6. 

87.21 

9.7’ 

16.53 

3  8.  77, 

9  •  .69 

99.29 

16. 

2.96 

782. 

89.67 

5.22 

23.75 

99  .33 

9  7.08 

99.se 

17. 

2.61 

593. 

'2.29 

8.05 

11.69 

).f>i 

64.66 

66.26 

12. 

1.89 

611. 

99.19 

7.92 

39.7  7 

66.25 

67.99 

60.63 

16. 

2.31 

626. 

>6.95 

13.29 

•53.01 

m. f; 

63.94 

66.22 

6  . 

.92 

632. 

97.38 

6.99 

69.9o 

66  .  22 

68.6? 

'0.36 

7. 

.96 

635. 

97.89 

9.78 

69.  *9 

>  '.96 

77.63 

'6.32 

-  . 

.  " 

69 

'?.fl 

9.0'* 

7  7.72 

'6.34 

7P.58 

60.89 

3. 

.96 

69  3. 

99.07 

6.56 

79.56 

6 '■•.89 

93.63 

66.9  3 

3. 

.96 

696. 

99.53 

7.99 

c  7.90 

86.-3 

99.69 

90.73 

0. 

0.00 

696  « 

99.53 

0.00 

o  7  .  0 

9"'.73 

93.31 

96.89 

2. 

.30 

6.8. 

95.89 

7.39 

9  0.25 

95.Rc 

96.09 

100.19 

0. 

0.00 

698. 

99.89 

0,00 

95.25 

‘  7  .  i  9 

101.69 

1.03.10 

1. 

.15 

699. 

100.00 

9.79 

100.00 
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TEST  NO.  1  i  DROP  SIZE  <MM)2 .  7 
RELATIVE  VELOCITY  (FT/SEC)  2 12 

FILMN0.865  MMO((i)7S 
SAMPLING  DISTANCE  !IN.)  30 


L.  P» 

M.  P. 

U.  P* 

"1.00 

3.01 

6»n2 

e."'2 

8.70 

11.39 

11.3? 

14.  CS 

16.77 

16.77 

19.9° 

2?.  22 

:  1.22 

26.69 

27.95 

27.9 

30.63 

33.32 

33.32 

36.01 

39.70 

38.73 

41.49 

44.20 

44.39 

47.09 

49.86 

49.98 

62.66 

55.25 

96.26 

57.94 

60.63 

60.6  3 

63.42 

66.22 

66.22 

68.59 

70.95 

70.96 

73.6? 

76.32 

76.32 

76. 63 

80.84 

83.34 

63.6  3 

66.4  3 

86.43 

88. 5a 

90.73 

90.  7  3 

93.31 

95.8? 

96.99 

98.04 

100.1° 

130.19 

101.6? 

103.20 

1 03.23 

ioe.69 

1 14,5? 

1 14.69 

117.29 

119.97 

119.97 

123, IP 

126.4? 

FO 

PFO 

SFO 

0. 

0.00 

0. 

135. 

22.13 

135. 

98. 

16.06 

233. 

99. 

14.59 

322. 

44. 

7.21 

366. 

5  7. 

9.34 

423. 

30. 

6.39 

462. 

25. 

4.09 

487. 

29. 

4.75 

516. 

16. 

2.45 

531. 

20. 

3.27 

551. 

13. 

1.63 

561. 

U. 

1  •  90 

572. 

1  ?. 

2.13 

585. 

3. 

.49 

589. 

9. 

1.31 

596. 

5. 

.91 

6C1. 

2. 

.32 

603. 

3. 

.49 

6C6. 

1. 

.16 

607. 

1. 

.16 

609. 

1. 

.16 

609. 

1. 

.16 

613* 

SPFO 

PMO 

SP«0 

0.00 

0.00 

0.00 

22.13 

.24 

.28 

38.19 

.63 

.67 

52.78 

1.60 

2.48 

60.00 

1.58 

4.06 

69.34 

3.51 

7.58 

75.73 

3.87 

11.46 

79.8? 

3.79 

15.25 

84.59 

6.40 

21.66 

87.04 

4.60 

26.26 

90.3  2 

8.20 

34.46 

91.96 

5.37 

39.84 

93.77 

7.46 

47.30 

95.90 

10.91 

56.22 

96.39 

3.05 

61.26 

97.70 

9.8? 

71.11 

98.52 

7.29 

78.41 

98.85 

3.41 

61.82 

99.34 

5.93 

67.75 

99.50 

2.20 

69.96 

99.67 

2.72 

92.68 

99.83 

3.38 

96.07 

100.00 

3.92 

100.00 
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TEST  NO.  12  DROP  SIZE  (MM)  2.  7 
RELATIVE  VELOCITY  (FT  SEC)  212 
FILM  NO.  872  MMD(r)  71 
SAMPLING  DISTANCE  (IN.)  30 
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L  •  P. 


o.oo 
6.02 
1 J .  39 
16.77 

23.22 

27.95 
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TEST  NO.  1  DROP  SIZE  (MM)  2.  7 
RELATIVE  VELOCITY  (FT/SEC)  209 
FILM  NO.  845  MMD  (|<)  59 
SAMPLING  DISTANCE  (IN.)  36 


M.  P. 

U.  P. 

FO 

PFO 

3.01 

6.02 

0. 

0.00 

8.70 

11.39 

147. 

23.00 

14.08 
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0.00 

117.23 

119.97 

0. 
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TEST  NO.  44  DROP  SIZE  (MM)  2.7 

RELATIVE  VELOCITY  (FT/SEC)  31^ 

FILM  NO.  928  MMO  ip)  32 

SAMPLING  DISTANCE  (IN.)  27 


I 

I 

I 


r>.03 
6.02 
5  I  .39 
.6.77 
23.22 
27,95 
33  •  32 
36.77 
7.4.  26 
A9.3r 
‘5.25 
6''. 63 
ft-  .2? 
73.9b 
76.  "<2 
S'1. 86 


v. 


o. 


3.71 

9.72 
16.08 
1  =  .79 
26.83 

<0.63 
'<6.7  1 
<*  1  ,«9 
<■7.08 
c2  «  66 
6  I  ,tU 

53. 63 

68.88 

73.63 
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TESTNO.  13  DROP  SIZE  (MM)2.  7 
RELATIVE  VELOCITY  (FT/SEC)  3  19 
FILM  NO.  866  MMD  (p)  36 
SAMPLING  DISTANCE  (IN.)  30 
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TEST  NO.  13  CHOP  SIZE  (MM)  2.  7 
RELATIVE  VELOGTY  (FT/SEC)  SI' 

FILM  NO.  873  MMO  (ft)  41 

SAMPLING  DISTANCE  (IN.)  30 

L.  P. 

.1.  P. 

U.  P. 

FO 

PFO 

SFO 

SPFO 

PMO 
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TEST  NO.  18  DROP  SIZE  (MM)  2.7 
RELATIVE  VELOCITY  (FT/SEC)  395 
FILMNO.  874  MMD(ft)  34 

SAMPLING  DISTANCE  (IN.)  24 
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TEST  NO.  48  DROP  SIZE  (MM)  *.7 
RELATIVE  VELOCITY  (FT-SFO  383 

FILM  NO.  929  MMD  (,')  19 

SAMPLING  DISTANCE  (IN.)  27 
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TEST  NO.  14  DROP  SIZE  (MM) 
RELATIVE  VELOCITY  (FT.  SEC) 
FILM  NO.  867  MMD((<)  24 
SAMPLING  DISTANCE  (IN.)  30 
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TEST  NO  14  DROP  SIZE  (MM)  2*  7 

RELATIVE  VELOCITY  (FT  SEC  394 
FILM  NO.  875  MMD  W  24 
SAMPL'NG  DISTANCE  (IN.)  30 
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TEST  HO  22  DROP  SIZE  (MM)  1.6 
RELATIVE  VELOCITY  (FT/SEC)  13  1 
FILM  HO.  879  MMD(|i>  94 
SAMPLING  DISTANCE  (IN.)  12 


L.  P. 

M.  P. 

U.  p. 

FO 

o.oo 
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0. 
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8.70 

11.39 
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11.39 

14.08 
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19.99 

23.22 
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25.58 

27.95 

26. 

2  7.95 

30.63 

33.32 

33. 

33  32 

36.01 

38.70 

12. 

36.70 

41.49 

44.2? 

9. 

44.29 

47.08 

49.88 

6. 

49.08 

52.56 

55.25 

4. 

55.25 
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60.63 

5. 

60.63 

63.42 

66.22 

7. 

66.22 
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3. 
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1. 
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0. 

66.43 

88.58 

90.73 

3. 

90.73 

93.31 
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0* 
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97.91 

3.42 

.48 
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0.00 
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o.or 
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0.00 
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0.00 

0.00 

622. 
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0.00 

0.00 

622. 
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0.00 

0.00 

622. 

99.03 

0.00 

.16 

623. 

1 jO.00 

22.5? 
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0.00 
.96 
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0.52 
12.95 
13.  ‘s 
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3  3 . 3  3 
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77.  - 

77.40 

77.40 

100. 00 


C14- 


TEST  NO.  21  DROP  SIZE  (MM)  1.  C 
RELATIVE  VELOCITY  (FT  SEC  130 
FILM  NO.  878  MMOW  73 
SAMPLING  DISTANCE  (IN.)  18 


L.  P. 

M.  P. 

U.  P. 

FO 

PFO 

SFO 

$PFO 

PMC 

' 

0.  CO 

3.01 

6.02 

0  . 

■\,ro 

A 

J  • 

0.00 

0.00 

o.oo 

6.02 

8.70 

11.39 

61. 

23.82 

6 1  • 

23.82 

.26 

.?o 

11.39 

14.09 

16.77 

5C  . 

19.53 

in. 

43.35 

.78 

1 .04 

16.77 

19.99 

23.22 

41. 

16.01 

152. 

i?."»7 

1.78 

2.82 

23.22 

25.58 

27.95 

13. 

5.C7 

165. 

fc-*.45 

1.12 

3.95 

77.96 

30.63 

3*. 32 

19. 

7.03 

183. 

71.46 

2.67 

6.62 

33.32 

36.01 

38.70 

18. 

7.03 

201. 

78.51 

4.31 

10.93 

38.70 

41.49 

44.29 

8. 

3.12 

209. 

61.64 

2.92 

13.85 

44.29 

47.08 

49.88 

13. 

5.07 

222. 

06.71 

6.91 

20.7b 

49.88 

52.56 

55.25 

10. 

3.90 

232. 

90.62 

7.38 

6  8.15 

88.25 

57.94 

60.63 

2. 

-.78 

234. 

91.40 

1.97 

30.12 

60.63 

63.42 

66.22 

5. 

1.95 

239. 

93.35 

6.47 

36.59 

66.22 

68.58 

70.95 

5. 

1.95 

244. 

95.31 

8.16 

44  •  76 

70.95 

73.63 

76.32 

3. 

1.17 

247. 

96.49 

6.06 

50.32 

76.32 

78.58 

80.84 

2. 

.76 

249. 

97.26 

4.90 

8  5.73 

SO.  64 

83.63 

86.43 

0. 

0.00 

249. 

97.26 

0.00 

65. 73 

86.43 

88.58 

90.73 

2. 

.78 

251. 

98.04 

7.02 

62.75 

90.73 

93.31 

95.89 

1. 

.39 

252. 

98.44 

4.10 

b6.36 

98.89 

98.04 

100.19 

1. 

.39 

25  ■>. 

93  .02 

4.76 

7  1  .62 

100.19 

101 .69 

103.70 

0. 

0.00 

253. 

98.62 

0.00 

7  1.62 

103.20 

108.89 

114,59 

0. 

9.00 

253. 

98.82 

0.00 

71.62 

1 14.89 

117.29 

119.97 

1. 

.39 

254. 

99.21 

6.14 

79.7o 

119.97 

123.19 

126.42 

1. 

.39 

255. 

99.60 

9.44 

6  9.2  1 

126.42 

128.78 

131.15 

1. 

,39 

256. 

1C-0.00 

10.76 

100.00 

-C15- 


TEST  NO.  23  DROP  SIZE  (MM  1.6 
RELATIVE  VELOCITY  (FT/SEC)  1 3  6 
FILMNO.880  «*»DM  93 
SAMPLING  DISTANCE  (IN.)  24 


L.  0. 

M,  P. 

U.  P. 

FO 

PFO 

SFO 

SPFO 

PMO 

SPMO 

0.00 

3.01 

6.02 

0. 

0.00 

0. 

0.00 

.  0.00 

0.00 

6.02 

8.70 

11.39 

21. 

19.09 

21. 

19.09 

.08 

.08 

'  1.39 

14.08 

16.77 

15. 

13.63 

36. 

32.72 

.23 

.31 

16.77 

19.99 

23.2? 

24. 

21.81 

60. 

54.54 

1.02 

1.34 

23.22 

25.58 

27.95 

4. 

3.63 

64. 

58.18 

.34 

1.66 

27.  95 

30.63 

'  33.32 

10. 

9.09 

74. 

67.27 

1.45 

3.14 

33.32 

36.01 

38.70 

4. 

3.63 

78. 

70.90 

.94 

4.06 

38.70 

41.49 

44.29 

1. 

.90 

79. 

71.81 

.35 

4.44 

44.29 

47.08 

49.88 

7. 

6.36 

86. 

78.19 

3.66 

0.10 

49.88 

52.56 

58.25 

2. 

1.81 

88. 

80.00 

1.45 

9.55 

85.25 

57.94 

60.63 

1. 

.90 

89. 

60.90 

.97 

10.52 

60.63 

63.42 

66.22 

1. 

.90 

90. 

61.81 

1*27 

11.79 

66.22 

66.58 

70.95 

3. 

2.72 

93. 

04.54 

4 . 0  1 

16.61 

70.95 

73.63 

76.32 

2. 

1*81 

95. 

66.36 

3.97 

20. 5t> 

76.32 

78.58 

80.84 

4. 

3.63 

99. 

90.00 

9.04 

30.23 

80.84 

83.63 

86.43 

2. 

1.81 

101. 

91.81 

5.61 

36.05 

86.43 

88.58 

90.73 

3. 

2.72 

104. 

94.54 

10.35 

46.40 

90.73 

93.31 

95.89 

0. 

0.00 

104. 

94.54 

0.00 

46. 4C 

95.69 

98.04 

100.19 

2. 

1.81 

106. 

96.36 

9.35 

50.76 

100.19 

101.69 

103.20 

1. 

.90 

107. 

97.27 

6.21 

60.97 

103.20 

108.89 

114.69 

0. 

0.00 

107. 

97.27 

0.00 

60.97 

1 14.59 

117.28 

119.97 

0. 

0.00 

107. 

97.27 

0.00 

60.97 

1 19.97 

123.19 

126.42 

0. 

0.00 

107. 

97.27 

0.00 

60.97 

126.42 

128.7.3 

131.15 

1. 

.90 

108. 

98.18 

10.59 

71.57 

131.15 

133.83 

136.52 

-» • 

0.00 

108. 

98.19 

C.00 

71.57 

136.52 

139.21 

141.90 

i  • 

.90 

109. 

99.09 

13.38 

•  04.96 

141.90 

144.69 

147.49 

1. 

.91 

no. 

IOC. CO 

15.03 

100.00 

-C16- 


TEST  NO.  10  DROP  SIZE  (MM)  1.  6 
RELATIVE  VELOCITY  (FT'SEC)  201 
4**  FILM  NO.  859  MMD  Qi)  45 

SAMPLING  DISTANCE  (IN.)  24 


l.  p. 

«.  ?. 

c 

• 

V 

• 

►  A 

opo 

ilT  ^ 

SPf 

°M0 

SPVG 

i.on 

3.0  1 

#•- 

6. "2 

0. 

**r.?C3- 

*  ...Of 

0.00 

O.QO 

0.00 

6.02 

8.70 

1  1.3? 

2  39". 

39.89 

239." 

39*89,*. 

2.02 

2.02 

i:  i’1’ 

14.08 

18.77 

137. 

22.87 

376. 

62.77 

4.22 

6.25 

15.77 

19.99 

23.22 

90. 

15.02 

466. 

77.79 

7.73 

13.98 

27,  2? 

25.58 

27.95 

48. 

8.31 

514. 

85.80 

S.21 

22.19 

27.95 

30.83 

33.32 

36. 

6.01 

550. 

91.81 

10.55 

32.75 

53.32 

36.01 

38. 7C 

12. 

2.00 

562. 

93.82 

5.67 

38.43 

^ 

41  .49 

•*4.29 

16. 

2.50 

577. 

96.32 

A  •  0 

‘.9.2o 

9 

•W.CS 

49.8)'- 

7. 

1.16 

584, 

97.49 

7.^ 

56.62 

49#  6* 

5  2.56 

5-3.2  6 

4. 

.66 

588. 

93.16 

5.83 

62.45 

55.25 

57.94 

60.63 

4. 

.66 

592. 

96.83 

latiO 

tO. 2b 

si. 63 

63.42 

66.22 

2. 

.33 

594. 

99.16 

5.11 

75.37 

55.22 

68.68 

7 O'. 9 5  '■ 

f.  " 

.  16 

595. 

99.33 

3.22 

78.60 

73.95 

7’. 6? 

7  6  .  2 

0. 

0.00 

•395. 

99.33 

C.00 

78.60 

75. 52 

78.58 

80.84 

2. 

.33 

597. 

99,66 

9.69 

o?.-»0 

89.84 

83.6? 

?6.4  3 

2. 

.33 

599. 

loO.OO 

11.70 

100. CO 

TEST  NO.  24  DROP  SIZE  (MM)  1.  6 
RELATIVE  VELOCITY  (FT/SEC)  206 
FILM  NO.  882  MMD  (?)  74 
SAMPLING  DISTANCE  (IN.)  30 


P. 

V. 

W.  7. 

FO 

PFO 

SFC 

£»FU 

f>W0 

j.:o 

3.01 

6.02 

0. 

0.00 

0. 

e.co 

0.00 

o.oc 

6.92 

3.70  ‘ 

11.79 

47. 

25.82 

4?* 

25.82 

.4  1 

.4) 

11  .  *v 

•4.64 

16.  "7^ 

41. 

22.32 

88. 

-5.76 

1.33 

1.75 

16.77 

19.99 

2  7.22 

22. 

12.08 

1 10. 

60.43 

1.99 

3.74 

2  1  •  •'  - 

25.56 

27.95 

17. 

9.34 

127. 

69. 78 

3.06 

6.80 

27. 

30. t>"5- 

.33.32 

17. 

9.34 

1  44  • 

79.12 

5.25 

1  2  .  C  5 

’1,7? 

76.91 

1  ?  *  '  u 

11. 

6.04 

155. 

o5.  1  6 

5.43 

17.54 

7  «  .  ’ 

4  J  .  -  9 

44.79 

4. 

4.94 

164. 

90.10 

£-.84 

<.  4 . 3  6 

-4.2'' 

4  7.95 

,..»3 

4. 

2.14 

168. 

•,*?.  '0 

4.47 

^  b  •  *  I 

4  9  ,  “  r 

52.56 

65.25 

3. 

1.64 

171. 

43.96 

“.4 ; 

*’  *  • 

55.2 

».  »  ^  ~  ^ 

;  ~  ' 

•  • 

.54 

1  72. 

-"*•50 

_  • . 

* 

,  •  • 

66 . 2  2 

2. 

1  .09 

1  7  4. 

4  6.60 

6.  Jr 

*  .  • 

-  •  - 

0. 

0.00 

1  74. 

*5,60 

C.00 

• 

7 *  #  6  * 

3. 

1.64 

1  77. 

47.25 

1 2 . 6  2 

'*  2 

»6.  <; 

,  3  • 

.83.84 

0. 

0.00 

177. 

4  7.76 

0.00 

::  ^ 

6  1  •  6  ? 

*36#**^ 

V  • 

0.00 

177. 

47.25 

o.oo 

”  ;  •  <* w 

5  A  .  “  < 

‘  r,  -  .  * 

40.  1  2 

i  • 

.34 

178. 

4  7.80 

7.30 

V  -  •  J 

-  '  ^  *  J 

' » *a  * 

‘  -  .>  .2 

2. 

1 .  C  9 

80. 

98.90 

17.09 

/  '  •  V  / 

- 

•  I  - 

0. 

0  .  '*0 

.60. 

46.90 

0.00 

•  * 

'  • 

'  • 

: .  4 

'2. 

1 jO.OC 

22.09 

■  \  • 

C17- 


TEST  NO.  9  DROP  SIZE  (MM)1  ■  6 
RELATIVE  VELOCITY  (FT/SEC)2C- 
FILM  N0.858  MM)  53 
SAMPLING  DISTANCE  (IN.)  36 


L.  P. 

V.  °. 

y.  P. 

*0 

PFO 

SFO 

S°FO 

PMO 

SPMQ 

0.00 

3.01 

6.02 

0. 

0.00 

0. 

0.00 

0.00 

o.oc 

6. Of 

8.70 

11.39 

38. 

27.94 

38. 

27.94 

.50 

.50 

11.39 

14.08 

16.77 

22. 

16.17 

60. 

44.11 

1.06 

1.56 

16.77 

19.99 

23.22 

20. 

14.70 

80. 

58.82 

2.69 

4.25 

23.22 

25.58 

27.95 

6. 

3.67 

85. 

62.50 

1.34 

5.60 

27.95 

30.63 

33.32 

10. 

7.35 

95. 

69.85 

4.59 

10.19 

33.32 

36.01 

38.70 

7. 

5.14 

102. 

75.00 

5.19 

15.38 

38.70 

41.49 

44.29 

5. 

3.67 

107. 

73.67 

5.65 

21.04 

44.29 

47.08 

49.88 

8  • 

5.88 

115. 

84.55 

13.18 

34.22 

49.88 

52.56 

55.25 

ii. 

9.08 

126. 

92.64 

25.15 

59.37 

66.25 

57.94 

60.63 

4. 

2.94 

130. 

95.58 

12.23 

71.60 

60.63 

63.42 

66.22 

3. 

2.20 

133. 

97.79 

12.02 

63.62 

66 .22 

68.58 

70.95 

2. 

1.47 

1  35. 

99.26 

10*11 

93.73 

70.95 

73.63 

76.32 

1. 

.73 

136. 

100.00 

6.26 

100.00 

TEST  NO.  26 

DROP  SIZE  (MM)  1.6 

RELATIVE  VELOCITY  (FT/SEC)  316 

FILM  N0.885 

MMD0J21 

24 

SAMPLING  DISTANCE  (IN.) 

L.  ?. 

v.  P. 

U.  P* 

FO 

PFO 

SFO 

S=FO 

PKO 

SPMO 

0.00 

3.01 

6.02 

C. 

0.30 

A 

c  • 

0.00 

0.00 

O.OC 

6.02 

8.70 

11.39 

1139. 

52.53 

1138. 

52.53 

10.0b 

10.08 

11.39 

14.08 

16.77 

528. 

24.37 

1666. 

76.91 

17.08 

27.17 

17  .77 

19.99 

23.22 

335. 

15.46 

2001. 

92. 

30.14 

?7.32 

23.22 

25.58 

27.95 

88. 

4.06 

2089. 

96.44 

15. 7d 

'3.10 

27.95 

30.63 

y*.32 

64. 

2.95 

2153. 

99.39 

19.66 

92.77 

35.32 

36.01 

38.70 

10. 

.46 

2163. 

99.86 

4.95 

47.73 

38.70 

41.49 

44.79 

3. 

.13 

2166. 

100.00 

2.26 

130.00 

-CIS- 


TEST  NO.  27  drop  SIZE  (mm)  1  6 

RELATIVEVELOQTY (PT/SECi  72.0 

FILM  NO.  886  MMO  (p)  30 

SAMPLING  DISTANCE  (IN.) 

30 

L  .  0  . 

M.  P. 

U.  P. 

FQ 

PFO 

SFO 

SPKO 

PMO 

SPmo 

o.co 

3.91 

6.02 

0. 

0.00 

0. 

0.00 

C#G3 

o.CG 

6.  02 

6.70 

11.79 

170. 

37.44 

170. 

37.44 

3*19 

i*  19 

11.39 

14.08 

16.77 

95. 

20.92 

265. 

78.37 

6*5  1 

9.71 

16. 77 

19.99 

23.22 

89. 

19.60 

374. 

77.97 

16*97 

2  0  •  &  3 

27.22 

25.78 

27.95 

74. 

7.48 

388. 

7;.  <• 

39.61 

27.97 

39.63 

33.72 

39. 

8.59 

427. 

94.07 

2 '#40 

66.0c 

73.32 

36.01 

38.70 

20. 

4.40 

447. 

48.47 

*3 1;  •  0‘6 

33.  70 

41.49 

44.29 

6  • 

1.32 

473. 

99.77 

9*61 

96*66 

94.29 

47.08 

49.88 

0# 

0.00 

453. 

99.77 

0*00 

95*66 

49.88 

72.56 

77.25 

0. 

0.00 

453. 

99.77 

0*00 

95*66 

75.23 

7  7.99 

60.63 

1. 

.22 

474. 

100.00 

4*33 

100.00 

TEST  NO.  30 

DROP  SIZE  (MM)  16 

RELATIVE  VELOCITY  (FT/SEC)  395 

FILM  NO.  891 

HMD  (4) 

18 

SAMPLING  DISTANCE  (IN.)  24 

!_.  ». 

M.  P. 

U.  P. 

FO 

PFO 

SFO 

Sr  F  0 

PMO 

SPMO 

7.70 

3.01 

6.02 

0. 

0.00 

0. 

o.oc 

0.00 

0.00 

6*^2 

8.70 

11.3  7 

B95. 

45.92 

895. 

47.92 

10.  Oi 

1  0  •  0  c 

1!  •  1 9 

14.06 

16.7  7 

659. 

33.76 

1553. 

79.68 

27.05 

37.1*. 

16*77 

19.99 

23.22 

311. 

15.95 

1864. 

95.63 

34.57 

72.71 

23*22 

27.58 

27.95 

52. 

2.66 

1916. 

98.30 

11.87 

84.5  7 

27*95 

30.63 

33.32 

27. 

1.38 

1943. 

99.69 

10.54 

95.11 

^  •  32 

36.01 

38.70 

4. 

.20 

1947. 

99.69 

2.5? 

97.63 

38*  70 

41.49 

44.29 

1. 

.05 

1948. 

99.94 

.96 

96.60 

44*29 

47.08 

4P.88 

1. 

.05 

1949. 

lJO.OO 

1.40 

100. oc 

TEST  NO. 

DROP  SIZE 

(MM) 

RELATIVE  VELOCITY  (FT/SEC) 

FILM  NO. 

SAMPLING  DISTANCE  (IN.) 

L  * 

M.  P. 

U.  °. 

FO 

“FO 

7f  ? 

i.3 

=  -  - 

5  -  V  , 

o*°“ 

?.l'l 

6.72 

7  . 

0.00 

0. 

0.07 

C.0& 

5*A2 

P.7n 

11.39 

387. 

79.17 

787. 

7.9.  •  ’ 

16.20 

1  6  •  ’  0 

11*39 

14.74 

16.77 

155. 

23.70 

742. 

H?.»7 

23.68 

jt),  g- 

16*7? 

19,99 

23.22 

91. 

13.95 

•  33. 

96.78 

36.68 

78.ro 

2  3*22 

2:  •  58 

27.95 

15. 

2.29 

.48. 

99.06 

12.71 

9 1 .2  3 

27.95 

30.63 

33.32 

6  * 

.91 

1.74. 

IrO.OC 

.1 . 7  0 

lOO-'V 

-C19- 


TEST  NO.  32  DROP  SIZE  (MM)  1.6 
RELATIVE  VELOCITY  (FT/SEC)  394 
FILMN0.9B  MMO^)  IS 
SAMPLING  DISTANCE  (IN.)  36 


L  * 

P. 

Mr 

,  p. 

|J.  P. 

FO 

PfO 

SFO 

»PFO 

PMO 

SPMO 

3. 

34 

.01 

6.92 

0. 

0.00 

• 

0.00 

0 

.00 

0.00 

B  4 

i  70 

11.39 

?99. 

60.77 

299. 

60.77 

2  j 

.95 

21.95 

i  *  • 

3'; 

14, 

.08 

16.77 

141. 

28.65 

440. 

89.43 

37 

.79 

59.74 

26m 

77 

19, 

.99 

23.22 

50. 

10.16 

490. 

99.59 

37 

.27 

97.02 

23. 

22 

26« 

.98 

27.99 

2. 

.40 

492. 

1 jC . 00 

2 

.97 

100.00 

-C20- 


TEST  NO.  55  DROP  SIZE  (MU)  0.  f, 

RELATIVE  VELOCITY  (FT/SEC)  69 
FILM  NO.  930MMOU)  100 
SAMPLING  DISTANCE  (IN.)  4-3/8 

.  P.  P.  FO  PFO  SFO  SPFO  PMO  h™ 


3. 

3.01 

A. 02 

6.  ">? 

8.73 

11.79 

U.^9 

14.38 

16.77 

if>.  n 

19.  « 

23. 72 

23.22 

25.58 

27.95 

27.  96 

30.67 

3  3.7.2 

13.12 

36.01 

3°. 70 

3fl. 7'. 

-1.49 

64.2° 

'*<..29 

47.08 

43.66 

4?. 86 

52.56 

56.75 

93.  1~ 

57.94 

60.6  3 

63.  *■! 

67.47 

66.72 

66.72 

66.68 

70.95 

7". 06 

73.63 

76.3  2 

76.3? 

7s. 68 

60,86 

80.84 

S3. 63 

bA.43 

66.6  1 

66.5s 

90.73 

95.71 

?3.71 

96. 69 

06.67 

96.06 

100.19 

1 30.1? 

101 .69 

103.70 

133.70 

1  Is .69 

116.69 

1  16.6  4 

117.76 

••  13.9  7 

113. 77 

123.1? 

1  /<■»(»? 

1 26.6? 

12-.73 

131.16 

131.16 

1  ’  ’  •  P  ’ 

i 36.-7 

1 ’6.6? 

179. ?1 

141.00 

161. 9  0, 

146*6? 

147.6'- 

1  6  7  .  /.  o 

150.78 

153. On 

1  *-3.0? 

1  5  *  •  7  r 

1 5s .4  5 

l 93.66 

161.1- 

163.8  3 

167.67 

1 Af .87 

169.47 

169.62 

171 .78 

1 74. 1  - 

1  ?  u  ,  1  6 

17-. 03 

17  J.s2 

1  70. 92 

1  H  !  .  78 

1 8- •  0  — 

IjA.O- 

18-.83 

!«  9.6  ' 

1  69.6  1 

1-31 .78 

1  3 .  4  1 

A 

«  • 

0.00 

0. 

11. 

31.80 

111. 

9-. 

26.93 

205. 

57. 

16.33 

762. 

21. 

6.01 

783. 

26. 

7.44 

309. 

70. 

5.73 

070. 

9. 

’.57 

338. 

3. 

.85 

341. 

?. 

.57 

343. 

0. 

0.00 

343. 

0. 

O.00 

343. 

7. 

.57 

345. 

0. 

o-.OO 

345. 

0. 

O.OO 

345. 

l. 

.28 

346. 

0. 

0.00 

346. 

1. 

.28 

347. 

0. 

0.00 

347. 

1. 

.28 

->48. 

0. 

0.00 

368. 

0. 

o.OO 

368. 

0. 

o.oo 

348. 

0. 

0.00 

3-8. 

r. 

0.30 

348. 

9. 

n 

’66. 

o. 

'I.CO 

366. 

o. 

O.OO 

’68. 

• 

0.00 

346. 

o 

/  • 

0.00 

34e. 

0. 

0.00 

36  P. 

0. 

4.00 

’4P. 

0. 

7.30 

348. 

0. 

0.00 

346. 

c. 

0.30 

348. 

1. 

.28 

349. 

0.00 

0.00 

0.00 

31.80 

.65 

.65 

58.73 

2.02 

2  .66 

75.07 

3»4l 

6.10 

61.08 

2.51 

o.6l 

38.53 

5.32 

1  '3 . 4  3 

V4. 26 

6.60 

20.4- 

96.84 

4.53 

25.07 

97.70 

2.20 

2  t.  2b 

98.28 

2.03 

29.31 

98.28 

0.00 

29.31 

99.28 

0.00 

29. ’1 

98.85 

4.50 

33.82 

98.85 

0.00 

33.32 

98.H5 

0.00 

33.82 

•<9.14 

4.08 

37.90 

99.14 

0.00 

3  7.90 

99.42 

5.66 

-3.57 

99.42 

0.00 

43.67 

99.7! 

7.32 

5  0  «  «  9 

99.71 

0.00 

50.64 

59.71 

0.00 

50.87 

99.71 

C.OC 

50.69 

99.71 

0.0  — 

5  0  •  8  9 

99.71 

0.00 

50.37 

99.71 

0.00 

50.89 

99.71 

0.00 

50*69 

99.71 

o.oo 

5  0  •  O  9 

99.71 

0.00 

5  0.55 

9° ,  7  1 

C.00 

"  ^ 

99.71 

0.00 

1/  0  •  **  * 

59. 71 

0.00 

99.7} 

0.00 

>  w  •  3  * 

7  9.71 

0.00 

79.71 

C.Op 

9  0  •  *  '  * 

1-0.00 

44.10 

1 ^0.^0 
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TEST  NO.  56  DROP  SIZE  m)  0.6 
RELATIVE  VELOCITY  (FT/SEC)  66 
FILM  NO.  931  MMOtf)  140 
SAMPLING  OISTANCE  (IN.)  4-3/8  . 

*.  a.  It.  P.  FO  JFO  SFO  SPFO  pvO  $PVj 


•  • 

7.71 

A*’'? 

■"1 

• 

6*0? 

6.7; 

11.79 

106. 

u .  i  -• 

14.38 

If. 77 

74. 

IS.  77 

1  7.9-s 

?7.  ?? 

33. 

7’.  7/ 

?6 ,4a 

27.36 

6. 

77. S- 

7'. 6  • 

77.7? 

1?. 

17.7? 

74. o; 

79.7/1 

6. 

18.7'. 

4  1.47 

44.79 

?. 

*4.  29 

4  7  ,5r 

49.98 

6. 

41.0  5 

5  ;  •  56 

56.78 

5. 

5a»  7 c' 

-7,3„ 

w>f* 

7.- 

47.47 

4  7.4? 

t*  *  ?  ' 

1. 

Sf.,7  7 

4C  .4  j 

1. 

71.4=) 

77.C.7 

74.7? 

1. 

76.7? 

1 7. .  78 

60.94 

0. 

4  -  . 

-7.61 

66.43 

0. 

86.4  * 

a  s  .  6  5 

99.73 

*  • 

9'.  77 

57.7  1 

56.89 

c. 

96.  F~, 

54 .14 

105.19 

0. 

l  on  .I-* 

171.65 

1  7 . 7  c 

1. 

1 33. 2  0 

114.65 

c. 

t  14.4 

117.76 

119.57 

0. 

1)9.5/ 

177.19 

t  2  .  4  2 

0* 

176.4/ 

139.7  = 

111.16 

0. 

’71.'6 

i 77.47 

1  74 . 5  . 

0. 

1 76. =7 

179. ?! 

141.30 

?. 

141,9 

154.63 

1  4  7.47 

0. 

\  U  7  •  u  ? 

150. ?" 

i  3  -.18 

9. 

* 

1—5.76 

If --.4  4 

0. 

1  »*  ?  *  C  - 

161,1- 

145.9  » 

7. 

I  A  "1  #  A  3 

1  6,4 ,4  7 

169.4? 

n 

•/  * 

1 VS.-7 

171.74 

174.16 

0. 

•  -*4. 1  * 

174.57 

i  7*5*  ^ ? 

l. 

9.30 

0.05 

0.00 

0.00 

40.00 

1C5. 

40.69 

.46 

.46 

28.06 

179. 

6V.37 

1.13 

1.85 

1  ?  .  79 

212. 

52.1  / 

1.48 

3.14 

2.32 

719. 

64.43 

.53 

3.57 

4,45 

7’3. 

63.14 

1.8? 

5.11 

7.3? 

?76.. 

91.47 

1.48 

7.00 

.77 

273. 

92.2** 

.  75 

7.75 

2.3? 

244, 

.4.57 

7.r'3 

1 1.3- 

1.9? 

?49. 

76*51 

1.81 

14.87 

.77 

261. 

77,78 

.7.04 

16.91 

.78 

76?. 

9  7,67 

1.71 

16.2“ 

.38 

253. 

99.06 

1 .66 

1  7.93 

.38 

254. 

76.44 

2.08 

22.02 

0.00 

264. 

96.44 

C.OC 

22.02 

3.00 

254. 

78.44 

0.00 

22.04 

5.00 

264. 

93.44 

0.00 

22.02 

3.00 

254. 

98.44 

0.00 

22.02 

7  .  ',‘0 

264. 

98.44 

0.00 

22.02 

.  36 

263. 

78.63 

5.4n 

/7.3l 

3. CO 

265. 

9S.63 

0.00 

27.51 

3.  CO 

285. 

98.93 

0.00 

27.51 

0."3 

258. 

33.83 

0.00 

27.51 

3.00 

765. 

>3.»7 

0.00 

*7.51 

9.CC 

255. 

58.87 

0.00 

27.61 

9  ,  C  5 

265. 

9«.p7 

0.00 

27.51 

0.00 

266. 

96.83 

0.00 

2  7.51 

'  •  .  0 

265. 

76.63 

0.00 

<  7.-1 

3. CO 

265. 

93.87 

0.00 

27.61 

.77 

?47. 

9  a  .  6 1 

43.45 

?  1.18 

3.00 

2C7« 

99.61 

0.00 

7  1.16 

9 . 79 

747. 

99.61 

0.00 

71.16 

.  36 

268. 

I  >)•(',> 

/  7  .  r  4 

100.00 
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TEST  NO.  57  MOP  SZE  (MM)  0.6 
RELATIVE  VELOCITY  (FT/SEC)  66  approx. 
FILM  NO.  932  MMD  (p)  115 

SAMPLING  DISTANCE  «N.)4-3/S 


L.  P. 

M .  P. 

U.  P. 

FO 

PFO 

SFO 

SPFO 

PMO 

SPMO 

*1.00 

3.01 

6.  *02 

0. 

0.00 

0. 

0.00 

0.00 

0.00 

6.02 

8.70 

11.09 

26. 

30.95 

26. 

30.95 

.15 

.15 

11.09 

16. OF 

16.77 

11. 

13.09 

37. 

66.06 

.26 

.60 

16.77 

19.99 

23.7? 

10. 

11.90 

67. 

55.95 

.61 

1.02 

20.22 

26. 5C 

27.96 

7. 

9.33 

56. 

66.28 

.86 

1.88 

27.96 

30.63 

33.02 

8. 

9.52 

62. 

73.80 

1.68 

3.57 

33.02 

36.01 

38.70 

1. 

1.19 

63. 

75.00 

.36 

3.91 

38.70 

61.69 

66.79 

2. 

2.38 

65. 

77.38 

1.03 

6.95 

66.29 

67.08 

69.86 

1. 

1.19 

66. 

78.57 

.75 

5.70 

69.86 

52.56 

55.75 

2. 

2*38 

68. 

80.95 

2.10 

7.80 

65.26 

57.96 

60.63 

3. 

0.57 

71. 

86.52 

6.21 

12.02 

60.60 

60.62 

66.72 

7. 

2.38 

73. 

66.90 

3.68 

15.70 

66.22 

68.58 

70.95 

r% 

• 

6.00 

73. 

86.90 

0.00 

15. 7C 

70.95 

73.63 

76.02 

3. 

3.57 

76. 

90.67 

d.62 

26.33 

76.32 

78.58 

8C.86 

2. 

2.38 

78. 

92.85 

6.97 

31.31 

80.86 

83.63 

86.6? 

1. 

1.19 

79. 

96.06 

6.21 

35.52 

86.63 

88.58 

90.7? 

1. 

1.19 

80. 

95.23 

6.99 

60.51 

?•".  70 

90.01 

95.89 

•9. 

0.00 

80. 

95.23 

0.00 

60.51 

95.89 

98.06 

100.19 

0. 

0.00 

80. 

95.23 

0.00 

•«0.51 

100.19 

101.69 

103.70 

0. 

O.CO 

80. 

95.23  . 

0.00 

60.51 

103.20 

138.89 

116.5? 

c. 

0.00 

80. 

95.23 

0.00 

60.51 

116.59 

117.29 

119.97 

1. 

1.19 

81. 

96.62 

11.59 

52. 1C 

129.97 

123.19 

126.67 

0, 

0.00 

81. 

96.62 

0.00 

52*10 

126.62 

128.78 

131.15 

2. 

2.38 

83. 

98.80  • 

30.67 

82.78 

101.15 

130.80 

136.5? 

1. 

1.19 

86. 

ICO. 00  ■ 

17.21 

100.00 
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TEST  NO.  35  DROP  SIZE  (MM)  2.7  DIAMETER  <|i) 

RELATIVE  VELOCITY  (FT/SEC)  138 
FILM  NO.  9  lb  MMDI*)  35 
SAMPLING  DISTANCE  (IN  )  15 


SSVH  lN3083d  3AUVnnwfr» 


-C25- 


1 


TEST  NO  41  DROP  SIZE  (I 
RELATIVE  VELOCITY  (FT/SI 
FILM  NO  0?  7  MMDOii  60 
SAMPLING  DISTANCE  (IN.) 


SSVHiN3DS3d  3AliVinWD 


TESTNO.  12  DROP  SIZE  (MM)  2.  V 
RELATIVE  VELOCITY  (FT/SEC)  212 
FILM  NO.  £12  MMO(ji)  71 
SAMPLINO  DISTANCE  (IN  )  10 


SSVW  lN3D83d  3Aliy-)r<*r  - 


-C30- 


RELATH  E  VELOCITY  (FT  'SfcO  i  I  S 
FILM  NO  'W  MMDW  3? 
SAMPLING  DISTANCE  (IN  )  7 


-C  >2  - 


5 

f 


TEST  Nl  1?  DROP  SIZE  (MM)  2.  7  DIAMETER  04 

RELATIX  E  VELOCITY  (FT/SEC)  -  19 
FILM  NC.$’?3  MMDfci  'rl 
SSMPLIMG  DISTANCE  (IN.)  30 


ssvw  iN3DU3d  3Aiivnnwn: 


TEST  NO.  -4-5  DROP  SIZE  (MM;  _\  7  DIAMETER  <„) 

RELATIVE  VELOCITY  (FT 'SEC'  V-i  5 
FILM  NO  9>l>  MMD  (»*)  19 
SAMPLING  DISTANCE  0 N  i  ?  7 


TEST  HO.  14  DROP  SI2E  (MM)  2.  7  DIAMETER  (u) 

RELATIVE  VELOCITY  (FT/SEC)  394 
FILM  NO.  367  MMOW24 
SAMP -INI  DISTANCE  fIN.)  30 


2g8S*S8$882o 


SSywlN3JJJ3d  3AllV"!nH(X) 


I 


\ 

1 


-C38- 


I 


TEST  NO.  23  DROP  SIZE  (MM)  j.  6  DIAMETER  (u) 

RELATIVE  VELOCITY  (FT 'SEC)  i36 
FILM  NO  MMD(n)  93 
SAMPLING  DISTANCE  UN.)  24 


•C4 


TEST  NO,  1)  DROP  SIZE  (MM)  1,6  DIAMETER  (,.) 

RELATIVE  '  El.OCITY  (FW  C)  ZO  1 
FILM  NO  66  »  MMD<jj)45 
SAMPLING  D'L'^NCE  (IN.)  24 


TEST  NO.  24  OROP  SIZE  (MM)  1.  6  DIAME  IBR  I* 

RELATIVE  VELOCITY  (FT  SEC>~,)6 
FILM  NO  882  MMD(|i)74 
•iAMpUNJ  DISTANCE  (IN.)  30 


Ml* 


RELATIVE  'ELOCITY  (FT  SEC)  LG2 
"IwMNO  *  P  MMD(f()  53 
SAMPLING  DISTANCE  (IN  ,  3f 


i  i-  i 

QQVM  I  3a  i  l - 


8 


8 


5 
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»- 

Ui 

3 
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10 
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- 
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— i  vC 
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< 

—  VL.  r- 

00  O 
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-C45- 


-C46- 


TEST  NO.  30  DROP  SIZE  (MM)  1.  6  DIAMETER  W 

RELATIVE  VELOCITY  (FT  'SEC)  395 
FILM  NO.  8’ 1  MMD(m)  18 
SAMPL.NO  DISTANCE  (IN.)  24 


SAMPLING  .DISTANCE  'IN  .  TO 


TFST  NO.  55  DROP  SIZE  (MM)O.  6  DIAMETER  (#,) 

PELATIVt  VELOCITY  (FT  SEC)  69 
FILMNO.  >90mMD{/<)  100 
SAMPUNC  CI5TANCE  ilK  )  1  1/ 8 


TEST  NO.  56  DROP  SIZE  (MM)  0.  6  DIAMETER  (p) 

RELATIVE  VELOCITY  (FT  SEC)  66 
FILM  N  >.  931  MMD  f/i)  MO 
SAMPLUG  DISTANCE  (IN.)  4-3 /8 


SSVH  lN3DM3d  3AUVinwm 

-C51- 


RELATIVE  VELOCITY  (FT'SEC;  66  .'.pprox. 
FILM  NO.  S  MMD  Ip)  113 
SAMPLING  DISTANCE  ON.)  4-3/8 


o 


•  I 


/o° 


.o 


«  < 


-  O 
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•  v 


i  0  •  cO 


af-  v 

'  <.  >.*.  *<*• 
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.  f  .*  •  «t'  ?  ».r 
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ft-  • 
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*4 


.  f  st  n',  *>  i>o*‘  b  ;f  *am  ? »’ 

M  1  AT  .'f  vfi  Of  »7f  ♦  T  »K  1  ? 
M  ft4'j  MVi  -  *0, 


IsfeKJS.-  •« 


-C52- 


*  4 


TEST  NO  I#  DROP  SIZE  (MADE 7 
RELATIVE  VClOOT"  -  |  SEC!  S9S 
Cl'  u  NO  SIR  MM!  iA 
SAMPLING  DISTANCE  i IN  '■  ZA 


TtST.-JOre  PROP  SIZE  <MMM6 
RELATIVE  VELOCITY  iFT  SEOlSl 
r  I_»1  wn  A  70  UMO  ,h\  94 

SAwri.  no  Distance  *n  '  12 


TEST  NO  DROP  SITE  WM  16 
RtL*.Ti  .£  VELOC'TY  FT  SfC  ?  1 

r  NO  959  wD  ,  4  5 

'AV“- '».G  D'cTA»jrf  •,  ?4 


TEST  nO  35  DROP  SIZE  **V)  0  * 
RELATIVE  v£lOC*T*  .fT  SEC: €9 
f.lMn''  *30  v*0  V  tOO 
:\*'L  *.G  DISTANCE  <tN  )4  3/4 


